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Introduction 


This Is the third annual report of the research activities con- 
ducted at the Materials Processing Center under the Materials Pro- 
cessing Research Base of NASA Grant #7645 during the period October 1 , 
1981 through September 30, 1982. Research projects conducted under 
the auspices of this grant Include work In the areas of nucleatlon 
and rapid solidification, fluid flow and crystallization processes, 
and adaptive materials processing. 

NASA's Initial and continuing funding of the Materials Processing 
Research Base has contributed slgnlflcately to the research and Indus- 
trial activities of the Center to develop and mature In an expeditious 
manner. The actual dollar support provided by NASA has been effec- 
tively leveraged to promote research and Industrial Interactions In 
areas of national Interest and of particular relevance to the Space 
Program. 

Since Its Inception In 1979, the Materials Processing Center has 
sustained very rapid growth. The Center's staff has grown to some 15 
faculty members, 30 research staff, and 56 graduate students, all 
drawn from a number of MIT Departments, Including Materials Science 
and Engineering, Chemical Engineering, Mechanical Engineering, and 
Electrical Engineering and Computer Science. The total funding level 
has grown to S3.425M, of which $1.267M Is derived from Industrial 
sources. 

A very Important activity which came to fruition In 1982 was the 
creation of the Materials Processing Center/ Industry Collegium, which 
has provided an excellent focal point for close Interaction with 
Industrial organizations. Founded on November 1, 1981, the Collegium 
currently has more than forty Industrial members and has a goal of 


one hundred members for 1983. The Collegium promotes Information and 
technology exchange through its Collegium Reports and Workshops. The 
Collegium also facilitates personnel exchange between industry and 
government agencies and MIT. In addition, the Collegium provides 
fellowships and scholarships to focus attention upon and attract top 
students to materials processing research. There are currently three 
fellowships for the 1982-1983 academic year, with commitments for four 

1983 summer scholarships and ten full-year fellowships for the 1983- 

1984 academic year. The Collegium also provides seed research funds, 
assists in acquiring necessary equipment, and supports the development 
of new urriculum and text books. 

Research efforts within the Materials Processing Center strive to 
control the structure, shape, and properties of materials at acceptable 
economic and social costs, through the effective utilization of the 
scientific and engineering base. The research activities of the Center 
may be divided into the following two broad categories: 

1) the Materials Processing Research Base, that is, work which 
is directly supported by NASA; and 

2) the overall materials processing research effort, which is 
supported by other governmental or industrial sources. 

The present report describes individual projects completed under the 
Materials Processing Research Base. Information pertaining to the 
general research effort will be contained in the 1982 Materials Pro- 
cessing Center Annual Report. While very substantial in its own right, 
the Materials Processing Research Base has further served as an 
extremely important nucleus for facilitating the growth of research 
efforts in related areas. 

The work described in this report, while involving research in 
the broad field of materials processing, has two common features: 
the problems are closely related to space processing of materials and 
have both practical and fundamental significance. An interesting and 
important feature of many of the projects is that the interdisciplinary 
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nature of the problem mandates complementary analytical modeling/experi- 
mental approaches which are best undertaken by several Investigators. 

The joint projects of Professors Brown and Witt, involving crystal 
growth, and Flemings and Szekely, dealing with rapid solidification 
and grain refining, are excellent examples of such cooperative efforts. 
Very close collaboration L 'tween investigators occurs in many other 
instances, although this may not appear formally in the project des- 
cription. Inaeed, the facilitation of such close cooperation is a 
major function of the Center. 

Another important aspect of many of the projects ftinded under 
the NASA grant is the increasing use of mathematical modeling tech- 
niques as one of the research tools. The predictive capability of 
these models, when tested against measurements, plays a very impor- 
tant role in both the planning of experimental programs and in the 
rational interpretation of the results. 

Finally, and perhaps most important, many of the projects des- 
cribed in this report will have a space experiment as their ultimate 
objective. Allowing for the realities of the 1980’s, thesr space 
experiments cannot be undertaken lightly, but require extensive 
ground-based preparations. Mathematical models are proving to be 
extremely valuable in projecting the findings of ground-based experi- 
ments to microgravity conditions. 

Of the projects funded by NASA, the Grain Refining/Rapid Solidi- 
fication program directed jointly by Professors Flemings and Szekely 
is the closest to realization as a space experiment. The space 
apparatus has been designed (on the basis of very extensive ground- 
base measurements and sophisticated calculations) and a formal pro- 
posal for the space experiments is being submitted. A formal proposal 
for space experiments on electromagnetically-driven flows in trans- 
parent media (Professor Szekely) will be submitted to NASA during the 
early part of 1983. Other space experiments on crystal growth under 
conditions of microgravity (Professors Witt and Brown) and electro- 
lytic processing (Professor Sadoway) are also being contemplated. 


Nucleation and Rapid Solidification 


i 
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Basic Study of the Role of Convection in Grain Refining 

Principal Investigators: Marton C. Flaming# 

and 

Julian Szakaly 


RESEARCH SUMMARY j 

.4 

■\ 

The scientific aim of this program is to obtain a better understanding 
of fluid flow, nucleation, grain refinement and structure in undercooled 
melts. Work under Professor 0. Szekely has the following principal components: 

i 

(i) Calculation of the electromagnetic force field, the fluid flow field 
and the temperature field in levitation melted (positioned) specimens both 
under earthbound and zero gravity conditions, and the verification of these 
through a comparison with experimental results. 

(ii) The use of the above techniques for the rational design of in-fliqht 

experiments, fully utilizing the unique features available in a zero gravity 
environment. | 

(iii) Fundamental studies of electromagnetic force field-turbulence inter- ] 

actions. < 

(iv) The interfacing of these results with the undercooling studies 

conducted under Professor Flemings' direction. j 

(v) The broad dissemination of the fundamental findings through potential ■< 

applications in metals processing. j 

•4 

I 

Work under the direction of Professor M.C. Flemings and Dr. Y. Shiohara j 

is in the following four areas: j 


Ik... 
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(1) Emulsification studies on undercooling and levitation melting Sn-Pb 
alloys. 

(2) Fine droplet studies on iron and nickel base alloys. 

(3) Thermal and structural measurements on levitation melted droplets 
of iron and nickel base alloys. 

(4) Emulsification studies on iron and nickel base alloys. 

To date, a total of 3 publications have resulted from this work; these are 

(7-9) 

listed at the end of this section of the report 


The Mathematical and Physical Modelling of Electromaqneticall.y Driven Fluid 
Flow and Associated Transport in Contained and in Containerless Melts 

(i) Electromagnetic Force Field and Fluid Flow Field Calculations 
Concerning Levitated Specimens 

A formulation has been developed to represent turbulent, recirculating 
flow as driven by the combination of buoyancy and electromagnetic forces in 
levitated specimens. The general vectorial form of the governing equation 
is given as: 


V*u = 0 


0 ) 


equation of continuity 


-V-pu u - Vo - Vt^ + pF b (2) 

where V*r is the turbulent stress tensor and is the body force field, which 
in general will include both the buoyancy and the electromagnetic forces. The 
electromagnetic component of the body force was calculated using the technique 
of mutual inductances^). The buoyancy component, which is of relevance omy 
for ground based applications was obtained through the simulataneous solution 
of the thermal energy balance equation. 

The thermal energy balance equation takes the following form: 

UVT = VkVT (3) 
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where k g is the effective thermal conductivity. POOR QUALITY 

The system of equations (1-3) was solved numerically, utilizing the k-e 
model to represent the components of the turbulent stress tensor. An extensive 
set of the computed results pertaining to ground based conditions has been 
published^ Here we S H a y| present a selection of the computed results 
concerning zero gravity conditions. In the presentation of the computed 
remits emphasis will be placed on the application of these to the interpretation 
of previous experimental measurements, notably the SPAR 1 series and on their 
use in the design of the Space Shuttle experiments, which involve scale- jp. 

SPAR 1 Experiments 

-5 

In the SPAR 1 experiments zero gravity (actually 10 g) conditions were 
produced, for brief time periods, say 5-7 mins, by using a rocket. In the 
studies beryllium spheres, 9 mm in diameter were levitated, using a spherical 
coil, sketched in Fig. (1). The coil was made of copper tubing, 2.5 mm in 
diameter. The coil consisted of two turns, each hemisphere wound on a 16 mm 
diameter sphere. The aximuthal angles formed by the two turns of the coil are 
46° and 74° respectively. The power supply was a 107 kHz radio frequency 
generator, which drew 1.2 kW, DC from a battery. The coil current was 300A 
and it was estimated that a power absorption rate of about 46 watts produced a 
mean sample temperature of 1400°C. 

Fig. 2 shows f Me computed electromagnetic force field and the heat 
generation pattern, while Fig. 3 shows the computed velocity and temperature 
fields. It is of interest to compare these results with the corresponding 
values computed for earthbound conditions using a coil assembly sketched 
in Fig. 4, for a steel sphere, 6 mm in diameter, given in Fig. 5 and 6. 

It Is seen that vigorous circulation exists in both cases, but for 
the zero gravity case we have four, rather than two recirculating loops 
in each half sphere. This significant difference between the two flow patterns 
Is not due to the absence of gravity, but rather attributable to the difference 
in coil configuration and the actual position of the sphere within the field 
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and hence the actual force field induced within the sphere which was depicted 
in Fig. 3 and 5. This point will be discussed subsequently. 

Fig^ 7 and 8 show the computed values of the turbulent kinetic energy 
and the ratio: effective viscosity/molecular viscosity for the zero gravity 
beryllium run and the earthbound steel run respectively. 

It is seen that in both cases the turbulent kinetic energy and the 
turbulent viscosity are higher in the center than at the outer surface* 
however there are marked differences in the details of the behavior. In the 
zero gravity case the isopleths are symmetrical, furthermore, the ratio of 
the turbulent viscosity to the molecular value is much hiqher and varies 
much less than in the earthbound case. The behavior of the turbulent kinetic 
energy parallels that seen for the turbulent viscosity. 

Space Shuttle Experiments 

In the proposed space shuttle experiments a variety of materials will 
be examined, hut for the purpose of comparison beryllium will be used in 
the present example. The same spherical coil design will be used as for 
SPAR 1, but the linear scale will be increased by a factor of two. 



In scaling these systems the following principal factors have to be 
considered. 

(i) A twofold increase in the linear dimension will provide a fourfold 
increase in the surface area of the specimen over which heat is dissipated. 

(ii) Approximate circuit analysis indicates that this scale-up wil! 
also increase the energy absorption due to the induced current by a factor 
of four. 

(iii) The above would indicate that the system could be operated at the 
same current, however due to impendance losses associated with the larger 
coil a higher coil current would be required. 


Calculation have shown that a coil current of 500A~ /2x300A would give 
an energy input of 255W, which in turn would give the desired 1400°C specimen 
temperature. 
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Fig. 9 shows the computed electromagnetic force field and heat generation 
pattern for this larger coil, which is quite similar to that computed for the 
smaller spherical coll. 

Fig. 10 shows plots of the computed velocity and temperature field' « 
four recirculating loops, characteristic of the spherical coll desigr «re 
readily apparent, but the numerical value of the absolute velocity smaller, 
by about a factor of four than that found for the larger coil. 

Fig. 11 shows the computed turbulent kinetic energy and the ratio: 
effective viscosity/molecular viscosity. It is seen that the flow field is 
highly turbulent but there is very little spatial variation of these properties. 
The important findings assocated with these results are that there are 
significant differences in the fluid flow fields in levitation melted droplets 
under earthbound and zero gravity conditions. 

More particularly for a zero gravity condition and the coil configuration 
employed we found a symmetrical circulation pattern and four recirculating 
loops in contrast to a asymmetrical circulation pattern and two circulation 
loops predicted for the earthbound conditions. 

The actual dimension of the drop was found to play a major role in 
determining the structure of the turbulence. For the small, earthbound 
system there was « marked damping of the turbulence on approaching the free 
surface, due to the changes in the angular momentum (the centrifugal and the 
Coriolis forces). In contrast for the larger, zero gravity systems the spatial 
variation in the turbulent kinetic energy was greatly reduced. 

Examination of the larger system (i.e. the planned Space Shuttle 
experiment) has shown that notwithstanding the increase in the absolute value 
of the coil current the actual linear velocity would be significantly reduced 
(by a factor of four compared to an identical geometrical arrangement but 
having half the linear dimension). 
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Our ability to quantify these effects is thought to represent significant 
contribution to the fundamentals of fluid mechanics and at the same time is 
an important asset in the rational planning of the in-flight experiments. 

(ii) The Design of Space Shuttle Experiments 

The principal scientific trust of the proposed space experiments is the 
extension of the ground based experiments to a zero gravity environment The 
actual concept of these experiments has been defined as a result of the prior 
ground based studies; the feasibility of accomplishing the measurements during 
in-flight conditions has been as certained as a result of discussions with 
Dr. T. Frost of the GE Space Science Center, in Valley Forge, Pennsylvania. 

While essentially the same equipment will be used for both the fluid 
flow and the undercooling studies, the details of the experimental objectives 
will necessarily differ. 

The purpose of the experiments is to determine the n It velocities on 
the free surface of levitated specimens at various current levels and thus 
enable a direct comparison of these experimental measurements with the theore- 
tical predictions that have been developed as part of this program. 

In essence the experimental procedure will involve the levitation melting 
of tin spheres, 1.84 cm in diameter, at various current input levels. The 
point values of the velocity at the surface of the tin specimens will be 
determined using tracer particles, the location of which would be followed 
using a cine camera, with a frame speed of 48 s’"*. 

The specific experiments to be performed on tin would include: 

(a) the determination of the steady velocity, corresponding to a 
particular current input level at several current inputs. 

(b) the rate of velocity decay, once the driving force has been drastical- 
ly reduced. 


Significance of the Experiments 

Fluid flow In levitation melted and positioned specimens Is an interesting 
class of problems in electromagnetically driven flows, because it Involves 
a curvi linear coordinate system. 

The proposed measurements would enable us to critically test the theore- 
tical representations that have been developed, through verifying - 

(a) the predictions made regard! the flow patterns and the velocity 
fields. 

(b) the dependence of the velocity field on coil current. 

(c) the turbulent - laminar transition. 

(d) the decay times, once the driving force has been discontinued. 

These experiments could not be carried out in an earthbound environment, 
because under the conditions of the earth's gravitational field there is a 
very narrow limit of coil currents, which allows one to attain levitation; 
almost invariably this corresponds to turbulent flow. Furthermore, one is 
limited in the size of the specimens that may be levitated. Under earthbound 
conditions it would be impossible to: 

(i) measure the current dependence of the velocity. 

(ii) measure the decay times. 

(iii) determine the laminar -turbulent transition. 

(iv) ascertain the fundamental difference ir the circulation patterns 
that have been predicted for earthbound and zero gravity conditions. 

In addition to testing the theoretical predictions pertaining to zero 
gravity levitation and positioning, the proposed experimental program would 
make an important contribution to the fundamentals of electromagnetically 
driven flows - a domain where accurate measurements are difficult to make. 

Thus the verification of the modelling equations should find rather wide- 
spread application in a number of related fields. 

Last but not least, the experimentally *nd theoretically predicted 
velocity fields will represent an important contribution to the characteriza- 
tion of the temperature and velocity fields needed in the modelling and quanti- 
tative description of the undercooling experiments. 
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(iii) Fundamental Studies of Electromagnetic Force Field - Turbulence 
Interactions 


The containerless processing program inherently involves the interaction 
of electromagnetic fields and fluid flow fields. One level of these interac- 
tions, namely electromagnetically driven flows has been discussed in preceeding 
sections, in connection with levitation melting. However, there exists 
another, more subtle level of interactions namely the possible damping of 
turbulence by an electromagnetic field. 

In order to examine this problem an apparatus has been constricted to 
study heat and fluid flow phenomena in a Woods Metal melt, generated by the 
passage o' an electric current between two electrodes. The apparatus is 
sketched in Fig. 12. An important component in this project is the use of 
a hot film anemoneter for the measurement of the velocity and of the turbulent 
fluctuations. A hot film anemometer assembly usable in molten Woods Metal 
has been developed for this particular application^. 

The calibration system for the hot film anemometer is given in Fig. 13; 
this is seen to consist of a shallow boat, containing the molten Woods Metal, 
through which the probe is being dragged at a predetermined, steady velocity. 
Quite sophisticated signal processing is required, the arrangements for which 
are shown in Fig. 14. A typical calibration curve is given in Fig. 15, 
indicating the drift, which is often associated with hot film anemometry in 
molten metal systems. 

Finally Fig. 16 shows the calibration curve with drift correction, which 
is quite satisfactory. 

Having established the feasibility of hot film anemometry in molten 
Woods Metal, the work is now continuing with systematic data gathering regarding 
velocity and turbulence paramenter measurements. The results of this work 
should be quite important in the containerless processing program, but should 
also find rather more widespread application in the general field of metals 
processing. 
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(iv) The Interfacing of the Heat and Fluid Flow Work with the Under- 
cooling Studies 

An important component of the modelling program is its association with 
the undercooling studies, conducted in Professor M.C. Flemings' laboratory. 
The detailed description of this work is available in Sections of this 
report. This effort is primarily experimental with the major thrust of 
examining undercooling, the subsequent recalescence of various alloy systems 
under a variety of conditions. Particularly noteworthy features of this 
work are that very significant undercoolings have been achieved, furthermore 
the thermal history of the recalescence process has been recorded, and very 
interesting metastable phases have been produced. 

Work is in progress regarding the modelling of the heat flow and the 
recalescence and metastable phase formation in these systems. The fluid 
flow and heat flow calculations, reported in the preceeding sections are of 
course closely related to this effort; more specifically work is in progress 
to define the fluid flow conditions both earthbocnd and for zero gravity 
in these undercooling experiments. Furthermore joint work has been initiated 
regarding the mathematical modelling of the recalescence process. 

( v) Industrial Applications of the Heat and Fluid Flow Calculations 

The tacxling of the heat and fluid flow problems in containerless 
processing operations has lead to the development of very powerful computa- 
tional tools capable of predicting the electromagnetic force fields, fluid 
flow fields, heat and mass transfer in a variety of metals processing opera- 
tions. An important, but of course not dominant component of our research 
program is to disseminate these results broadly to the metals industry, 
thus providing direct industrial applications of the NASA sponsored basic 
research program. 

Examples of this work include the study of fluid flow in induction 
durnaces, and other steel processing operations, described in the preceeding 
report of the Center for Materials Processing. 
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A more recent industrial application, which is illustrated in the 
following has been the prediction of the velocity fields and mass transfer 
rates in a 4 ton induction furnace holding molten steel. 

Fig. 17 shows a schematic sketch of the furnace, while the computed 
electromagnetic force field and the velocity fields are given in Fig. 18, 

19 respectively. Fig. 20 depicts the experimental arrangement used for 
measuring the velocities in the melt, which essentially relied on the 
deflection of an immersed rod. Finally, Fig. 21 shows a comparison of the 
experimentally measured and the theoretically predicted dependence of the 
linear melt velocity on the coil current. The good agreement, for pilot 
scale type operations, is readily apparent. 

Studies were also made to predict the rate at which graphite rods, 
immersed in the melt would dissolve, and Table 1 shows a comparison 
between the predictions and the experimental measurement - here again the 
agreement is excellent. 

UNDERCOOLING EXPERIMENTS ON TIN-LEAD, AND IRON AND NICKEL BASE ALLOYS 

Emulsification Studies on Sn-Pb Alloys 

Emulsifications are prepared using a high-speed stirring device patterned 
after Perepezko et al.. Fig. 22- Some of the emulsifications were reheated 
and the metal droplets then resolidified at controlled cooling rate in a DTA. 
With the droplet sizes and experimental technique explained, undercoolings 
of about 95 to 105°C were typically obtained. Other samples were rapidly 
solidified in the device shown in Fig. 23. Thus in these cases, undercoolings 
of not less than about 110°C were combined with rapid cooling rates (calculated 
to be in the range of 10 6 °C/sec.). 

The nucleation temperatures of the different alloys studied, when cooled 
at 10°C/min. are shown in Fig. 24. T p represents nucleation of a tin-rich 
phase, and T ' of a lead-rich phase. 
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Sn-5 wtX Pb dropelts, undercooled and rapidly solidified in the agitated 
bath, show "featureless" structures. Fig. 25, which are supersaturated solid 
solution. A lead-rich pahse precipitates in the solid state after extensive 
holding at room temperature or heating above room temperature. Fig. 26. 

Having a cold bath temoperature for quenching is essential to obtain the super- 
saturated structure. Fig. 27. 

In Sn-Pb alloys of compositions where two phases nucleate, structures 
are obtained that are strongly sensitive to cooling rate, such as that of 
the Sn-45 wti Pb alloy in Figure28. When rapidly quenched, the solidification 
structure is composed of a dendritic, lead-rich phase (white, with dark solid- 
state precipitates) and a dark tin-rich pahse. When slowly cooled (10°C/min.), 
the lead- rich phase coarsens as shown in Figure 28 on the far right. 

Analytical studies are now underway as part of this research to under- 
stand and model the nucleation, solidification and heat flow behavior of 
these undercooled droplets, and there will be reported at a later date. 

I 


i 

I 
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FINE DROPLET STUDIES ON IRON AND NICKEL BASE ALLOYS i 

We have developed a method for obtaining high undercooling in metal 
alloys by mixing fine metal powder with glass powder, and then heating and 
cooling the composite. To date, we have achieved 475°C undercooling in 
316 stainless steel in this way - a substantially higher undercooling than 
has been obtained by any other worker for an iron base alloy. The former j 

maximum reported for Fe base alloys was 295°C, Table 2. - 

i 

As undercooling increases and particles size decreases, structures 
obtained in 316 tend to go from dendritic FCC; to "featureless FCC"; to 
mixed FCC/BCC; to full BCC. Figure 29 (top) shows the dendritic FCC micro- 
structure characteristic of very low undercoolings. Figure 29 (bottom) 
shows the "featureless" microstructures characteristic of higher undercoolings. 

Dark particles are BCC; lighter ones are FCC. Mixed FCC/BCC structures are 
also obtained. Figure 30. (BCC structure surrounds the smaller FCC area). 

These results are in general agreement with expectations, and calculations 
confirm that in this alloy, the BCC structure can be formed (metastably) 
at undercoolings greater than about 120°C. 
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TH ERMAL AND STRUCTURAL MEASUREMENTS ON LEVITATION MELTED DROPLETS OF IRON 
AND NICKEL BASE ALLOYS 


Extended experiments have also been carried out on levitation melting, 
jn solidification of levitation melted droplets, and on structures of 
droplets produced in this way. 

In current work, iron and nickel base alloy droplets are floated in 
Molten glasses to achieve higher undercoolings. Surface temperatures of 
■coalescing, and of subsequently cooling droplets are being measured to 
provide high-speed data for modelling studies. Figure 31 is a schematic 
illustration of the equipment, showing the sample, the high speed sensing 
device, and the digital oscilloscope. 

Figure 32 is an example of short time thermal prof i els obtained in a 
Ni -32. 5 wt% Sn (eutectic) alloy. It is an interesting and startling 
observation that the total recalescence time (at least of the surface of 
the sample) decreases sharply with increasing undercoolings. Figure 33 
shews longer time thermal profiles (up to 10 sec.) for these droplets; 
this figure ,/'Ovides an indication of the time required to complete recales- 
ciince. Note that the rapid surface cooling rates are obtained immediately 
after solidification (> 10^ °C/sec.) even in these very large droplets 
(5 mm = 5,00C ym). This rapid cooling is due to extraction of heat from 
the recalesced droplet into the much cooler glass. Much more of these data 
are being obtained and being used as foundation for developing a mathematical 
model both o. solidification of these large droplets, and of the smaller 
droplet- 3-50 ym) studied in other portions of this work. 
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,MP 1 1 CAT IONS FOR SPACE EXPERIMENTS OF POOR QUALITY 

(1) Fine Metal Droplets Suspended in Glass 

Tne fine droplet approach of mixing metal powder with glass and heating 
is an effective way of obtaining higher undercooling, in spite of one experi- 
mental probelm. That problem is coarsening or agglomeration of the particles. 
For example, starting with an initial average particle size of 15 ym, the 
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typical average particle size after melting and resolidification is approxi- 
mately 90 pm, an increase in droplet volume of over a factor of 200. This 
size increase, we believe, is much too large to be accounted for by diffusional 
coarsening. It most probably results largely from agglomeration due to gravi- 
tatinal settling in the glass, both before and after the glass is molten. 

Further experiments are perhaps necessary to confirm this probable 
gravitational cause of agglomeration. Assuming this is the correct mechanism, 
an important space experiment becomes self-evident; that is, the carrying 
out of similar experiments in space. The absence of a gravitational field 
should permit retaining very fine particle sizes inside the glass, thereby 
obtaining much higner undercoolings than we have yet been able to achieve. 

( 2) Levitation Melted Droplet in Glass Blanket 

Our current large diameter dorplet work (5 mm dia.) is on nickel base 
alloys levitated with a glass slag. Thus far, undercoolings obtained have 
been appreciably less than undercoolings obtained in the fine droplets, but 
the technique has great interest for us because we can measure the surface 
temperature of the droplet during recalescence using high speed optical 
methods. 

Performing this experiment in space would permit eliminating both 
gravitational convection and induction stirring (since the power would not 
be needed to maintain levitation). Thus, we would have to be able to achieve 
increased undercoolings in this space experiment and could also measure 
the rates of recalescence in samples at these undercoolings, recording on 
telemetering the digitized information. 
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Fig. 1. Schematic sketch of the levitation coil 
used in SPAR 1 experiments. 
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Computed electromagnetic force field and 
heat generation distribution for an 18 
mm beryllium drop at zero-g. 
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Fig. 12. Schematic sketch j of the apparatus designed 
for measuring the electromagnetically driven 
velocity field In the Woods Metal system. 
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Fig. 14. Sketch of Data Processing Scheme. 
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Fig. 15. Calibration curves showing drift. 
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Fig, 16. Calibration curves with drift correction. 
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Fig. 19. The computed velocity field for a 4 ton 
Induction furnace and for a 2800 A coll 
current. 
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Fig. 20. A sketch of the velocity measurement system 


1 - probe 

2 - ann 

3 - spring 

4 - strain 


5 - support 

6 - instruments 

plate 7 - x-x-T recorder 

gauge 8 - RMS converter 
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Fig. 22. Schematic diagram of the experimental apparatus 
for emulsifying droplets of mclten low-melting- 
point alloys. 
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Composition of Itsd (st%Pb) 


20 40 60 80 Pb 



Composition of Itsd (wt*/»Pb) 


Fig. 24. Composition dependence of nucleatlon temperatures 
for Sn-Pb alloy droplets observed In DSC measure- 
ment with 10°C/m1n cooling rate. 
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Figure 25 Microstructures of droplets of Sn-5wt%Pb alloy 
ropidly quenched in a bath at -25°C, 1000X. 
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Figure n 6 The same droplets shown in Figure 25. after one 
month at room temperature, note the substantial 
precipitation of lead-rich phase, 1000X. 
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(a) Bath Temperature: 50 U C 


(b) Bath Temperature: 25 C 


(d) Bath Temperature 


(c) Bath Temperature: -4 C 


Microstructures of droplets of ?n-5wt%Pb alloy 
superheated 5(TC and rapidly quenched in bath 
at temperatures shown, 5000X. 
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Table 2 

Previous Undercooling Results for Iron and Nickel Base Alloy 


Alloy 

Undercool inq ( °C) 

Investigators 

Sample Size 

Ni 

319 

Cech, Turnbull 

50 - 100 pm 

Ni 

365 

Ward 

1-5 gm 

Fe 

295 

Cech, Turnbull 

30 - 100 ym 

Fe 

295 

Scheil 

Bulk melt 

Fe-25%Ni 

300 

Kattamis, Flemings 

50 - 500 gm 

Fe-25%Ni-0. 3%S 

225 

II 

fl 

AISI 4330 Steel 

200 

ft 

II 

440 C. Steel 

150 

II 

II 

Ni 

200-290 

Lux, Haour, Mollard 

12mm0 x 10cm 

Ni-15%Co 

200-240 

I* 

ft 

Ni-12%Co-15%W 

130-145 

11 

II 

Ni-1 3%A1 (Ni 3 Al) 

10-15 

II 

II 

Ni-42^50%Sn 

<138 

Kobayashi, Nakao 

3 gm 

Ni-52%Sb 

229 

II 

tl 
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Fig. 29. la) SEM micrograph of cross-sectioned and 
etched particles of the as-received 
316 stainless steel powder. 

(b) Optical micrograph of cross-sectioned 
and etched 316 BCC and FCC particles 
after undercooling experiment. (FCC 
particles are lighter in color). 
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Fig. 30. SEM micrograph of a cross-sectioned and 
etched particle of 316 after high under- 
cooling. Both the BCC and FCC phases 
co-exist in the same particle. 
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Time (sec) 


Fig. 33. Eutectic solidification profiles for 
three different undercoolings. 
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Tht Structure and Propartiaa of Rapidly Solidified, 
Highly Alloyed Aluminum Materials 

Principal Investigator: Nicholas J. Qrant 


I. THE STRUCTURE AND PROPERTIES OF RS-PM 2024 A1 + LITHIUM ALLOYS 
RESEARCH ABSTRACT 


This study has been completed and a series of publications are in 
preparation. One phase of the study was concerned with Alloy 2024 modi- 
fied with additions of 1 to 2 % lithium. Further alloying modifications 
permitted study of the effect of the Cu:Li ratio on the resultant me- 
chanical oroperties. Alloy 1, incorrectly treated at 500°C for hydrate 
decomposition, was the high Cu, high Cu:Li ratio alloy, and was nega- 
tively influenced by formation of coarse CuAl^ precipitates (in excess 
of solid solubility at all temperatures). Alloys 2 and 3 showed excel- 
lent combinations of both nominal and specific strength and modulus of 
elasticity. Fatigue properties, notch toughness and crack propagation 
rates were all significantly improved over I-M 2024 and RS 2024 without 
lithium additions. 

A publication covering this study has been accepted Tor the Second 
International Conference on Al-Li Alloys, scheduled for April 1983. 


RESEARCH SUMMARY 


Table 1 lists the Cu and Li compositional variations, Young's Modulus 
and the shear modulus values, which in turn are transposed into percent 
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increases in both modulus values corrected for specific gravity, 
increases are of significant value and interest. 


These 


The major strengthening phases in these three alloys are summarized 
below for aging conditions of 20°C (T-4) and 19G°C (T-6): 


Aged at 2Q°C 

Alloy 1 GPB and GP zones 

Alloy 2 GPB zone 

Alloy 3 GPB zone 


Aged at 190°C 

GPB zone, s' , y' , 8' and e" 
GPB zone, s' and y' 

GPB zone, s' and y' 


Alloys 2 and 3 are similar in structure but otherwise differ signifi- 
cantly in properties. One would expect higher levels of E and G, and of 
specific modulus values in favor of Alloy 2 with the higher lithium con- 
tent, which is confirmed in Table 1. Yield strength values for all three 
alloys are closely matched at about 68 to 69 ksi; however, the UTS is sig- 
nificantly better for Alloy 3 in spite of having the lowest lithium con- 
tent. Perhaps even more importantly, the ductility at fracture *>s sharply 
better for Alloy No. 3. The combination of excellent YS, UTS and ductility 
shown by Alloy 3 is related to its lower content of oxide derived from the 
RS powders. The fracture surfaces showed a cleaner, finer, dimpled struc- 
ture. Alloy 1, which had very coarse CuAl^ precipitates due to initial 
overheating of the RS powders gave pocr strength and ductility values and 
is not further discussed. 


The density corrected values for Alloys 2 and 3 are excellent; the 
best combination is that for Alloy 3 for the T-6 condition. See Table 2. 


The notch tensile properties and fracture toughness values are shown 
in Table 3. Again Alloy 3 shows he best values. 


Figure 1 is a plot of stress amplitude (reversed bending) versus 
cycles to failure at room temperature in an air atmosphere, and compares 
data for I-M 2024, RS-splat 2024 (no Li), and for Alloys 2 and 3. There 
is a significant improvement for Alloys 2 and 3 (uncorrected for density) 
compared to both I-M 2024 and RS-splat 2024. In the case of Alloy 3, the 


55 


QRiGtNML '.'d 

OF POOR QUALir-f 



•» 

UJ 

£C 

3 


12 


o 

h- 

co 

UJ 

-I 

o 

> 

o 


(!S>j) 3anilldlMV SS3H1S 


56 


Fig. 1. Reversed Bend Fatigue Data, Smooth Bar, for RS-PM Alloys 2 and 3, 
for Splat Quenched 2024 and for IM 2Q24-T4. 




improvement at 10 7 cycles, over I-M 2024, Is more than 35%. 

Figure 2 plots fatigue crack growth rate versus stress intensity 
range for Alloys 2 and 3 versus that for I-M 2024. Again Alloy 3 shows 
excellent values! 

Alloy No. 3, with only about 1% Li, and apparently with lower oxide 
content than is frequently the case with RS-PM and RS-splat based alloys, 
shows excellent properties in all tests performed, including the sharp, 
precracked notch tests. Corrected for density, this alloy is clearly 
better than I-M 2024 in terms of its mechanical behavior. 


Table 1. Modulus of Elasticity of RS-PM Aluminum Alloy 2024 + Li 


Alloy No. 

Cu 

Li 

E 

G 

Percent Increase in Values 
vs. I-M 2024 


% 

% 

psi x 10® 

psi x 10 7 

AE 

AE 

AG 





P 

x 10 7 i 

JL £L 


I-M 2024 

4.5 

- - 

10.4 

3.9 


-- 


-• 

1 

5.92 

1.29 

11.3 

4.3 

8.7 

10.8 

10.3 

12.9 

2 

3.78 

1.63 

11.7 

4.5 

12.5 

18.3 

15.4 

20.5 

3 

4.16 

0.96 

11.3 

4.3 

8.7 

12.5 

10.3 

12.8 


°®ginai 
0F poor 


page 

Quality 
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CRACK GROWTH RATE, (10 in/cycie) 



CRACK GROWTH 








Table 2. Room-Temperature Tensile Properties of RS-PM 2024 
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II. THE STRUCTURE AND PROPERTIES OF RS-PM Al-Mg-Li ALLOYS 


RESEARCH ABSTRACT 


In terms of anticipated gains in specific strength and elastic modu- 
lus, the Al-Mg-Li system is particularly attractive. Even though the lit- 
erature shows that Al-Cu-Li alloys have better strength values (nominal 
and specific), the substitution of 2 to 4% Mg for 4 to 6% Cu in the A1 - 
1 to 3 % Li alloy is expected to be particularly beneficial to specific 
modulus. Rapid solidification, which has not been researched to an impor- 
tant extent for Al-Mg-Li alloys, was applied in this study in order to pro- 
duce extruded bar stock for test purposes. The first alloy and test results 
for an RS-PM alloy look quite promising. 


RESEARCH SUMMARY 


The composition of the hot extruded (30:1 ratio) Al-Mg-Li alloy was 
(in weight percent): 6.74Mg - 1.64Li - 0.19Cu - .27Cr - .21Ti - .13Fe and 
.15Si. Solution heat treating temperatures of 300° to 540°C were studied, 
along with aging from 20° (T-4) to 190°C (T-6). 

Dendrite arm spacing values for USGA powders (helium quench) varied 
from about 0.5 to about 2pmover a powder size range which varied from 

5 

about 1 or 2 ym to 250 pm, indicating a quench rate of about 10 K/s. 

Aging in Al-Mg-Li alloys is essentially based on formation of the 
6 ( A 1 3 Li ) phase; the magnesium serves primarily to enhance solid solution 
strengthening. As a secondary effect, the magnesium decreases the solu- 
bility of Li in A1 , thereby increasing the amount of 5' phase on aging, 
for a given lithium content. 

The results of exploratory aging studies for the Al-Mg-Li alloy are 
plotted in Fig. 1. Briefly, solution temperatures below about 300°C do 
not induce important aging benefits. Natural aging results in no change 
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in hardness. Best results are for a solution temperature of 510°C plus 
aging at 135°C (curve 7, Fig. 1). Of some interest but not pursued, solu- 
tion temperatures from 400 to 510°C, with aging at 165°C, result in a nar- 
row spread of hardness values. In this respect, solution temperatures 
from 400 to 510°C are apparently sufficient to achieve suitable supersat- 
uration; however, aging at 165°C leads to overaging after about 10 hours. 
Aging at 135° to 150°C seems preferable. 

In all cases the 5 ' (Al^Li) phase was identified after aging, although 
some A^Mg phase is expected for this 6.7% Mg alloy. 

With an addition of only 1.64% Li, the density of this alloy is 
3 

2.523 g/cm , compared to 2.768 for I-M 2024; this represents a density de- 
crease of 8.85% for the Al-Mg-Li alloy compared to I-M 2024. In turn. 
Young's modulus for this alloy is 11.2 x 10^ psi, yielding a value for the 
specific Young's modulus of 12.2 x 10^ in. These are quite attractive val- 
ues. 

Table 1 lists room- temperature tensile data including specific 
strength values based on the noted decrease in density. Measured strength 
and ductility values, in general, are better than those for I-M 2024 - 
T-4; corrected for density the values are excellent and are as good as 
those reported for Alloy 3 in Section I above. 

Figure 1 indicated optimum hardness on aging at 135"C; this is borne 
out in Table 1 in terms of yield strength; however the ultimate strength 
is not the highest among the tests which were made, largely due to consid- 
erably lower ductility values for this aging condition. Aging at 150° or 
165° results in lower yield values, equally high ultimate strength values 
and progressively improved ductiliy. In fact * tion at 450° plus aging 

at 165°C yields the best combination of stren and ductility while pre- 

serving a Notch Yield ratio of 1.0. Crack growth rates were the poorest 
among the group of alloys discussed in Sections I and II, for reasons 
which are not now apparent. 
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Nevertheless, it is clear that further studies of RS-PM Al-Mg-Li 
alloys are warranted, in particular in compositions which will provide an 
additional hardening effect to that provided by 4' precipitation in this 
alloy. 


Table 1. Tensile Properties of A1 - 6.74Mg - 1.64Li Alloy at 20°C 


Solution and Aging 
Conditions 

0.2% YS 

UTS 

Elong 

RA 

Spec. YS 

Spec. UTS 

O 

O 

Ksi 

Ksi 

% 

% 

4 

p x 10 in. 

4 

-p x TO in. 

510°, 1/2 h, WQ + 
135%, 170h 

64.0 

75.7 

5.5 

5.0 

70.4 

82.3 

510°, 1/2 h, WQ + 
150°, 15h, WQ 

63.4 

74.6 

7.0 

10.0 

68.9 

81.2 

510°, 1/2 h, WQ + 
165°, lOh, WQ 

60.2 

75.1 

8.9 

10.0 

65.4 

81.6 

450°, 1/2 h, WQ + 
165°, 1 Oh , WQ 

61.2 

76.7 

11.0 

15.0 

66.5 

83.4 

400°, 1/2 h, WQ + 
165°, lOh, WQ 

61.6 


8.0 

10.0 

67.0 

83.7 


III. MODIFICATION OF STRUCTURE AND PROPERTIES OF RS-PM 7075 ALLOY 
CONTAINING 1% ADDITIONS OF Zr AND Ni 


RESEARCH ABSTRACT 


In general, a great deal of the alloy studies using RS particulates 
leads to small compacts weighing between 1.5 and 5 kg, measuring 4 to 15 
cm diameter to produce hot extrusions typically 12 mm diameter to 12 x 35 
mm cross section, and representing extrusion ratios (area) of 16-30:1. 

This is often an inadequate amount of hot working to disperse oxide films 
on the particulates, or to break down occasional trapped impurities. In- 
dustrial R & D powder billets, in contrast, will weigh 50 to 100 kg and 
will receive considerably greater hot reduction. Using USGA powders, a 
35 kg batch of powders will be hot extruded to achieve from 60 to 100:1 
extrusion ratios (area). Various subsequent cold and hot work cycles will 
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be tried In order to achieve finer and more uniform dispersions of oxide 
films and dirt which may accidentally be Included in the powder batch. 


RESEARCH SUMMARY 


Excellent structures and mechanical test data were obtained from small 
extrusions (5 kg) of an RS 7075 alloy containing 0.7% Zr and 1% N1 . Yield 
strength values in excess of 90 Ksi, UTS values in excess of 100 Ksl and 
ductility values in excess of 10%, combined with excellent fatigue, frac- 
ture toughness, and notched tensile data qualified this alloy for further 
studies to evaluate the role of total reduction on structure and subsequent 
properties. The USGA alloy powders have been produced and are currently 
being canned in preparation for extrusion from an 18 cm round to a 5 cm 
round (13:1 area reduction). The 5 cm extrusion will be cut to selected 
lengths and variously processed to smaller cross-sections, for example: 

a) Hot re-extrude to 1.2 cm: total final reduction equals 225:1 

b) Hot upset the 5 cm round to obtain an alternate grain, oxide, and 

inclusion pattern at very high reduction 

c) Other TMT's including Intermediate cold work, hot rolling and cross 

rolling, etc. 

Based on expected beneficial effects, an alloy much more prone to de- 
velop heavier oxides, such as a 2 or 3% Li alloy, will be similarly proc- 
essed. 
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Development of Rapidly Solidified, Oxidation Rasistant Alloy? 

Principal investigator: Gregory Yurek 


ABSTRACT 


A fine-grained rapidly solidified 18-8 stainless steel (AISI 3-3) 
exhibited enhanced resistance to oxidation under isothermal conditions 
compared with a conventional wrought alloy of the same composition. The 
enhanced resistance to oxidation is attributed to three factors: (1) T he 
formation and growth of protective scales is promoted by the fine al y 
grain size (^10ym) and probably by the presence of a dispersion of fine 
MnS precipitates in the steel matrix, (2) scale adherence is increased by 
the formation of intrusions of oxide from the external scale into the alloy, 
which form along closely-spaced alloy grain boundaries and probably around 
the dispersed phase, and which act to key the scale mechanically to the al- 
loy, and (3) growth stresses in the scales are reduced by the formation of 
scales of uniform composition, which are promoted by the uniform composi- 
tion of the alloy. Modification of the microstructure of 18-8 stainless 
steels through rapid solidification processing appears to allow an increase 
in operating temperature of about 75°C under isothermal conditions. 

INTRODUCTION 


Very fine grained alloys are typically obtained by compaction of alloy 
powders produced by rapid solidification processing. In a number of cases, 
the fine grain sizes (<10ym) remain stable at elevated temperatures. The 
grains in these alloys are apparently pinned by very fine nonmetallic in- 
clusions that are formed in the alloys during rapid solidification and/or 
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during subsequent compaction of rapidly solidified products. The com- 
bination of a stable, fine alloy grain size, a fine, uniform dispersion 
oi nonmetallic inclusions and a homogeneous alloy composition should 
provide an excellent alloy microstructure for oxidation resistance, pro 
vided the alloy composition is adjusted properly. 

The objective of this program of research was to determine the ef- 
fects of rapid solidification processing on the high-temperature oxi- 
dation behavior of alloys. It was demonstrated in a previous in'esti- 
gation [1] that a f .e-grained, rapidly solidified 18Cr-8Ni stainless 
steel exhibited superior resistance to oxidation at 9G0°C under cyclic 
conditions compared with a conventional, wrought allr- of the same com- 
position. The present study focussed on the oxidot on behavior of the 
same rapidly solidified alloy under isothermal conditions over the tem- 
perature range of 800-1 100°C. 


EXPERIMENTAL METHODS AND MATERIAL S 

The material selected for this investigation was a fine-grained rap- 
idly solidified 303 stairless steel (designated herein as 30CRS), which 
contained a fine dispersion of MnS precipitates. Rapid solidification 
of this alloy was carried out by forced convective cooling (in helium) 

of centnfugally atomized metal droplets to achieve cooling rates on the 
5 

order of iO °C/sec. The resultant powder, which was 5 to 200 microns in 
diameter, was sieved to -140 mesh, enclosed in a stainless steul can,, and 
consolidated to a fully dense one-inch diameter rod by hot extrusion at 
900°C and a 10:1 extrusion ratio. According to Kelley and Vander Sande 
[2], the 303RS alloy comprised a duplex grain structure of 5-10 micron 
grains surrounded by heavily dislocated grains that were less than one 
micron in diameter. The alloy contained about 1 vol % o° MnS in the form 
of a uni form dispersion of fine precipitates (0. Cl -0.2 microns in dia- 
meter) . 

The oxidation behavior of the 303RS was compand to that of a con- 
ventional wrought 304 stainless steel (designated herein as 304W). The 
nominal composition of both steels, which are given in Table 1, are vir- 
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tually the same, except that the 303RS has a higher content of sulphur 
and silicon. A wrought 303 stainless steel was not used as the stand- 
ard for comparison because that material contains very large Mn3 preci- 
pitates which would disrupt the oxide scale and cause rapid oxidation. 

TABLE 1: Compositions of Steels 
Composition, Wt. Pet. 

Steel C Mn P S SI Ni Cr 

30 4 W 0.04 1.62 0.033 0.015 0.13 8.78 18.28 

303RS 0.059 1.60 0.028 0.34 0.62 8.68 17.21 


The oxidation specimens were disks approximately 1mm thick and 16 
and 11.5mm in diameter fer the 303RS and 304W alloys, respectively. These 
disks were sealed in an evacuated quartz capsule and anneal led for 2h at 
1100°C. The specimens were then ground through 600 grit SiC paper, pol- 
ished with 6 micron diamond paste, ultrasonically cleaned and then fin- 
ally rinsed with distilled water and methanol. After the stabilization 
anneal, the av rage grain sizes of the 303RS and the 304W were 6.5 and 
85 microns, respectively. The MnS precipitates were 0.2 to 0.4 microns 
in diameter, indicating that particle coarsening occurred during anneal- 
ing. 


Oxidation experiments were conducted using a thermogravi metric 
apparatus that employed an Ainsworth Semi -Microbalance as the weighing 
unit. Experiments were conducted at 800. 900, 1000 and 1100°C in an 
atmosphere of dry oxygen at 1 atm pressure. Temperature control was 
-5°C as measured by a Pt/Pt-10%RH thermocouple placed approximately 5mm 
from the specimens. The microbalance was interfaced with a strip-chart 
recorder that produced a contin;. s record of the weignt gained by the 
specimens as they oxidized. >i? a '.curacy of the balance was -.02 gm. 

The weight gains (£m) measured by the balance were then normalized by 
dividing by the specimen area (A). 

Scanning electron microscopy (SEM) was used to examine the surfaces 

and polished cross-sections of the t..idized specimens in order to deter- 
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mine the oxide morphology. The relative proportions of the elements in 
the oxides were determined from the ratios of the 1^, peaks in EDX spectra. 
An electron microprobe was used to measure concentration profiles on the 
polished cross-sections of oxidized alloys. 


RESULTS 

Values of (Am/A) are plotted as a function of time in Fig. 1 for the 
isothermal oxidation of 304W and 303RS steels at 800-1100°C. The total 
times of exposure varied, but the tendency was to terminate the tests after 
breakaway oxidation occurred. Two tests were conducted on the 303RS mat- 
erial at each temperature. The data discussed throughout this report are 
from whichever of the two tests resulted in the largest weigh, t grins. 
Figures 2 and 3 show the complete sets of data for both tests at each tem- 
perature. Reproducibility of results was poor at 1100°C, but at all lower 
temperatures the results of the two tests were within - 20» of each other. 
Data for the 304W steel were produced by a single test at each temperature. 

The kinetic data show that the 303RS alloy was more resistant to oxi- 
dation than the 304W alloy over the entire range of temperatures tested. 

The general trends may be summarized as follows. The 304W formed protec- 
tive scales at 800 and 900°C; at 1000°C a protective scale formed on 304W 
but soon broke down as thick iron oxides began to grow on the surface and 
pitting occurred. At 1100°C, breakaway oxidation began almost immediately. 

By comparison, the 303RS alloy formed protective scales at 000, 300, 
and at 1000°C. Even at 1000°C there was no evidence of pitting o r scale 
breakdown. At 1100°C, the oxidizing conditions were too severe and break- 
away oxidation set in quickly. The most significant difference between 
the two alloys is that the J03RS was resistant to oxidation at temperatures 
up to 1000°C, while the 30 A W was only resistant to oxidation up to 900°C. 
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DISCUSSION OF RESULTS 


Oxidation of 304W Alloys 

The maximum recommended temperature for the use of a wrought 304W 
stainless steel under isothermal oxidation conditions is 925 0 C [3]. 

Therefore, it was expected that the 304W would show adequate oxidation 
resistance at 900°C and below, but not at 1000°C and above. 

The scales observed on specimens oxidized at 800 and 900°C consist- 
ed of thick, relatively iron-rich oxides over the grain interiors and thin, 
relatively chromium-rich oxides over the alloy grain boundaries. A protec- 
tive chromium- rich oxide also formed beneath the iron- rich oxides over the 
grain interiors and was subsequently exposed when many of the oxide mounds 
spalled during cooling. When the stainless steel was initially exposed to 
the high- temperature, oxidizing atmosphere, oxides of every metal in the 
alloy probably formed, including the rapidly growing wustite, "FeO". Wustite 
can grow so quickly that a thick, Fe-rich scale formed first, followed by the 
formation of the underlying protective chromium-rich scale at the alloy sur- 
face. The formation of wustite in the early or transient stage of oxidation 
was apparently responsible for the rapid initial weight gains c' the 304W 
specimens (Fig. 1). The chromium- rich oxides at the grain boundaries were 
able to form rapidly enough to prevent the formation of thick iron-rich oxides 
at these sites. This phenomenon is attributed to the preferential nucleation 
of CrgOj at the boundaries and to the rapid transport of Cr to the alloy sur- 
face by grain boundary diffusion. 

Oxide-filled pits that occasionally developed on the 304W/900°C specimen 
were able to reform a protective Cr-rich layer at the base of the pits. The 
formation of such pits is generally attributed to local cracking or spalla- 
tion of the protective scale [4,5]. Once formed, the pits can continue to 
grow until they reach a depth in the alloy where the chromium concentration 
is above the critical level for the formation of a protective film [5]. Be- 
cause there were relatively few pits and a protective scale reformed at 
the base of the pits, this steel is sufficiently oxidation resistant to be 
useful up to 900°C. 
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At 1000°C the 304W began to passivate, but then breakaway oxida- 
tion kinetics set in as iron-rich oxide mounds and pits began forming 
over a significant portion of the surface. These pits usually did not 
reform a protective scale at their Dase, indicating that the underlying 
alloy no longer had a high enough chromium concentration to reform a con- 
tinuous layer of Cr^O^. In addition, the very steep chromium concentra- 
tion profile in the alloy near the interface (Fig. 4) and the high propor- 
tion of iron in the scale (30-70 wt%) indicate that the supply of chromium 
to the surface was probably insufficient to maintain a protective scale. 
That is, the oxidation conditions were so severe that the scale degraded 
to a less protective structure or composition even in the absence of crack- 
ing or spallation. Therefore, catastrophic degradation of the alloy is ex- 
pected and observed at this temperature and at higher temperatures. 


Development of Protective Scales on 303RS Alloys 
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The 303RS alloy exhibited much better oxidation resistance than the 
304W, although both had the same nominal composition. The 303RS specimens 
had a relatively smooth Cr-rich protective scale at temperatures up to 1000°C, 
although at temperatures of 900°C and above the scale was dotted with oxide 
mounds that were rich in iron and chromium. However, even the scales with 
the oxide mounds were distinctly more uniform than the scales found on the 
304W alloys. The oxide mounds appeared to be the '>mnants of rapidly grow- 
ing iron-rich oxides that formed during the initi cages of oxidation. 
Polished cross-sections of the specimens revealed that all mounds that 
formed on 3Q3RS at or below 1000°C were undercut by a protective chromium- 
rich scale. The formation of pits was commonplace on 304W at 10Q0°C, but 
they were not observed at all on the 303RS specimens. 


The preferential formation of a uniform chromium-rich scale on the 
303RS alloy is attributed to its fine grain size and the presence of a fine 
dispersion of MnS in the alloy matrix. Grain boundaries and second-phase 
particles that intersect the surfaces of alleys may be sites for the prefer- 
ential nucleation of C^O^. Since the 303RS exhibits both a fine grain size 
and a uniform dispersion of MnS p. ecipitates , it is expected that there would 
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be a high density o* closely spaced Cr^ nuclei on the alloy surfaces 
during the initial stages of oxidation. Therefore, relatively little 
sideways growth of these nuclei would be required to form a continuous, 
protective chromium oxide scale. The grain boundaries are also paths for 
rapid diffusion of chromium to the surface of the alloy. This increases 
the availability of chromium in the vicinity of the grain boundaries and 
promotes the formation of protective chromium oxide films in these areas. 

The protective chromium-rich oxides formed over grain boundaries in 
the 304W stainless steel at 800 and 900°C were approximately 10 microns 
wide. Since the grain size of the 303RS alloy is only 6-8 microns, the 
presence of a 10 micron wide protective film along the grain boundaries 
would cover the entire surface of the alloy. 


Maintenance of Protective Scales on 3Q3RS Alloys 

The 303RS alloys were more effective both at forming and maintaining 
protective Cr-rich oxide scales on their surfaces. The breakdown of Cr- 
rich scales may occur by spalling or cracking or by chemical destabiliza- 
tion caused by chromium depletion of the base metal and the simultaneous 
incorporation of iron or manganese into the scale, which causes a trans- 
formation to a less protective scale. 

The uniform composition of the 303RS alloy appears to promote the for- 
mation of external scales of uniform composition and thickness, which would 
tend to decrease stresses and stress gradients in the scales, thereby en- 
hancing the ability of the scales to resist spallation. The 304W alloys 
tended to oxidize nonuni formly, which probably contributed to breakdown 
of protective scales and promotion of the growth of iron-rich oxides. 

Although the tendency for cracking of a chromium-rich scale may be 
reduced on 303RS, it would be expected that occasional breaks in the scale 
would jccur. However, no pitting of the alloy surface or rapid growth of 
Fe-rich oxide mounds was observed at 1000°C and below. Therefore, it is 
assumed that any breaks in the scale that may have occurred were rapidly 
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healed by the formation of a new protective scale. Apparently, the inter- 
face of Cr content of the 303RS alloy (13.7 wt%, see Fig. 5) was not be- 
low the critical level for reformation of a protective scale at 1000°C. 
Rapid transport of Cr by grain boundary diffusion in the 303RS alloy, 
which had retained its fine grain size, apparently accelerated reforma- 
tion of the protective scale relative to the rate of reformation of protec- 
tive scales on 304W. 


Development Si0 2 Subscales and Spinel Layers on 303RS 


The protective scales on the 303RS alloy we^e enriched in Si at the 
alloy/scale interface and in Mn at the scale/gas interface. Both of these 
effects are commonly observed on commerc'al stainless steels. 

The kinetic data and EDAX analyses of polished cross-sections in- 
dicate that the formation of Si0 2 subscales was more prevalent on 303RS 
than on 304W and occurred more readily at higher temperatures. These ef- 
fects were expected since 303RS has both a significantly higher concentra- 
tion of silicon (0.62 vs. 0.13 wt%) and a finer grain size (6.5 vs. 85 
microns) to enhance Si transport to the outer surface. 

The silica subscales were highly convoluted and formed intrusions -into 
the alloy along the alloy grain boundaries where the transport of oxygen 
and silicon is believedto be more rapid [1]. The silica intrusions have 
been observed to increase the adherence of oxide scales to tr.e substrate 
[1,6]. Both the subscale and the intrusions decreased the oxidation rate 
since further oxidation was dependent upon the transport of reactants 
through the silica scale as well as through the chromium-rich outer scale. 

Analysis of the kinetic data for the 303RS specimens oxidized at 900 
and 1000°C demonstrated that the rate of oxidation durinq approximately the 
first 40h was controlled by diffusion in Cr^, followed by a transition 
to an even slower rate of oxidation. This transition is believed to be 
caused by the deve 1 opment of the subscale on these alloys. 
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The manganese-enriched outer scales frequently had the approximate 
composition of the manganese-chromium spinel, MnCr^. The formation of 
this spinel probably consumed some of the protective Cr 2 0 3 and produced 
the faceted crystals that were observed on the 303RS alloy surfaces. The 
spinel also appears to have spalled frequently upon cooling, particularly 
at 1000°C. 

The formation of a manganese enriched outer scale layer on the 303RS 
is due to the high affinity of manganese for oxygen and the rapid diffusion 
of Mn through C^Og [7]. The preferential transport of Mn through the chro- 
mium-rich outer scale may be assisted by the presence of an Si0 2 subscale 
which serves as a differential diffusion barrier to inhibit the transport 
of chromium [1] while allowing the flux of manganese and iron [5] to con- 
tinue unabated. 


Long-Term Oxidation Resistance of the 303RS Alloy 

The protective oxide scales which form on the 303RS alloy are expected 
to exhibit very little cracking or spallation during long-term oxidation. 
The uniformity of thickness and composition apparently minimizes internal 
stresses in the oxide, thereby limiting the incidence of cracking of the 
scale. At the same time, the oxide intrusions which key the scale to the 
alloy help prevent spallation. 

In the event of scale breakdown after the alloy has been in service 
for a long period of time, continued resistance to oxidation requires that 
the alloy be able to reform a protective surface scale. Therefore, i* is 
essentia' that the alloy retain the microstructural characteristics which 
were responsible for its original oxidation resistance. In the case of 
303RS, these characteristics include a fine dispersion of MnS precipitates 
to pin the alloy grain boundaries and to serve as heterogeneous nucleation 
sites for Cr 2 0 3 , and a fine grain size to promote both heterogeneous nucle- 
ation of CrgOj and Cr transport from the bulk alloy to the surface. 

Optical microscopy was performed on the 303RS specimens to determine 
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the grain size after oxidation at 900, 1000 and 1100°C (Table 2). The 
measured grain sizes indicate that the alloy strongly resists grain 
coarsening at 1000°C and below, and less strongly at 1100°C. Transmis- 
sion electron microscopy , studies of the alloys after oxidation were not 
performed, so the size and distribution of MnS particles is not know. 

However, Yurek et al . [1] observed that the MnS particles in a 303RS alloy 
did not noticeably coarsen during 3Q0h at 900°C, hut remained 0.2 to 0.5 
microns in diameter. 

Olson et al . [8] presented a model for the growth of grains in alloys in 
which the grain boundaries are pinned by second phase particles. This 
model predicts that the grain coarsening process will obey a cubic rate 
law of the form: 


( r - r Q ) = k t 1/3 

where r is the size of the grains, r Q is the original grain size, k is a 
constant, and t is time. The grain growth which was measured during the 
oxidation tests (Table 2) was used to calculate k for 900, 1000, and 1100°C. 
The model was then used to predict the alloy grain size expected after 1000 
and 10,000 hours at each of these temperatures. The results a^e also pre- 
sented in Table 2. The model predicts that a relatively fine grain size 
(13.9 microns) will be maintained in the 303RS alloy even after 10,000 
hours at 900°C, so the enhanced oxidation resistance of this alloy should 
be maintained for relatively long periods of time. At 1000°C grain growth 
proceeds more rapidly, and after 10,000 hours a grain size of approximately 
28 microns is predicted. This large degree of grain growth will almost 
certainly have an adverse effect on the long-term oxidation resistance of 
the alloy, but the magnitude of the effect cannot be predicted. It should 
also be noted that the silica intrusions at the surface of the alloy will 
still serve as mechanical pegging points even after grain growth, and will 
actually limit the severity of grain growth that is possible in the near- 
surface region. 
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TABLE 2: 

Predicted Grain 

Growth 

of 303RS 

Grain Size 
10,000 hrs. 

Thermal History 

Measured 
Grain Size 

k 

Predicted 
1000 hrs. 

900°C/200 hr 

8.5 micron 

0.34 

9.9 

13.9 

1000°C/135 hr 

11.7 

1.01 

14.9 

28.3 

1100°C/50 hr 

22 

4.48 

46 

92 


(The initial grain size was 6.5 microns) 


Improvements in the long-term oxidation resistance of rapidly solid- 
ified alloys may be achieved if the fine grain sizes inherent in these alloys 
can be stabilized. Recent work at M.I.T. [9] demonstrates that the control 
of small quantities of impurities in these alloys, which eventually yield 
very stable, complex oxide precipitates that are very effective in pinning 
grain boundaries, may be of paramount importance in the processing sequence. 
More studies in this area would be beneficial because additional knowledge 
of how to control alloy microstructure, and hence alloy properties, through 
rapid solidification processing would help significantly in the develop- 
ment of new alloys with improved properties. 
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Figure 1. Kinetics of the isothermal oxidation of 303RS 
and 304W steels. 
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Figure 2. Kinetics of the oxidation of 303RS at 800 to 
1000°C. 
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Figure 3. Kinetics of the oxidation of 303RS at 1100°C. 
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Concentration profiles in a 304W alloy that 
had been oxidized for 126h at 1000°C. The 
abscissa is the distance from the alloy/scale 
interface. 
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Figure 5 


Concentration profiles in a 3C1 3 F;c alloy 
that had been oxidized for I ' • n 1000°C 

The abscinsa is the distance from the 
alloy/sca! '^terface. 
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Processing of Sendust -Typ« Soft Ferromagnetic Alloys 

Principal Investigator: Samuel A. Allen 


RESEARCH ABSTRACT 


Super Sendust is a quaternary iron-base with 5.3 w/o Si, 4.2 w/o AT 
and 3.1 w/o Ni. It was developed from a ternary iron-silicon-aluminum 
alloy, Sendust, with the aim of improvements both in soft ferromagnetic 
properties and in mechanical ductility. The present study focussed on 
processing of Super Sendust by a combination of melt-spinning rapid- 
solidification, deformation, and thermal annealing techniques. The 
objective was to obtain material in the form of thin ribbons, with 
modest ductility, and excellent soft ferromagnetic properties. 

The melt spinning of Super Sendust was shown to give ribbon 
geometries that was strongly dependent on nozzle geometry of the melt 
crucible. The most rapidly cool ribbons were ductile in the as-cast 
edition. As cast, the ribbons had a coercive force H c of about 1.9 
Oe. Two-stage heat treatments were devised that enabled H c to be reduced 
to 0.5 Oe, but in addition embrittled the material. The optimum 
annealing schedule consisted of one hour at 1200 * 0 , cooling at 
-10°C/minute to 690°C, holding at 690*C for three hours, and air 
quenching. Transmission electron microscopy was used to studv micro- 
structures of the as-cast and annealed ribbons. A phase transition 
from the B2 to DO 3 superlattice structures was found to occur 00 cooling 
at about 680%. Magnetic properties of samples with the 82 structure 
were superior to those with the DO 3 structure. 


INTRODUCTION 


The Sendust composition, Fe-9.6 w/o SI-5.4 w/o Al, has long been 
recognized as an all'jy with exceptional characteristics as a soft ferro- 
magnetic material. These are: High permeability, low coercivity, high 
electrical resistance and reasonably good saturation magnetization. 

These excellent properties result primarily from the fact that the com- 
position has near-zero magnetostriction as well as near-zero magnetic ani- 
sotropy. Unfortunately, the material is also very brittle and applica- 
tions in the past have been limited mainly to cast shapes and to forms 
made from compressed powder. 

Two recent studies give promising results with respect to 
processing-related improvements for these alloys. The first is a report 
by Tsuya et al. [1] of the use of a melt-spinning rapid-solidification 
technique to produce Sendust alloy ribbons of fine grain size. As cast, 
these ribbons were ductile enough to be rolled to 25% of their original 
thickness. After annealing treatments to optimize magnetic properties, 
however, the material was embrittled. The second study, by Yamamoto and 
Utsushikawa [2] focussed or alloy modification through 3 w/o additions of 
nickel to an Fe-6 w/o Si-*r w/o Al alloy. Although the mechanical proper- 
ties to the Super Sendust alloy were not studied, there is reason to 
expect that it will have better ductility than Sendust: Nickel additions 

to soft ferromagnets are generally believed to contribute to both magnetic 
and mechanical softness. Also, Super Sendust has approximately one-third 
less aluminum and silicon, each of which embrittles iron. 

The aim of the present study was to employ rapid solidification 
techniques to produce Super Sendust alloy ribbons, and to study 
processing-stnr -re-property relations of the resulting material in order 
to optimize soft ferromagnetic properties. The results of the study are 
reported in detail in the Master's thesis of Eric P. Kvam [3]. 
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Alloy Preparation and Melt Spinning 

Super Sendust alloy buttons were prepared by vacuum arc melting. 

The ingots were then remelted and cast in a melt spinning apparatus onto 
either the inside or the outside of a rapidly rotating drum. Ribbons 
rangeed in thickness from 30 to 100 pm, and in width from 2-5 ron. Ribbon 
quality proved difficult to control. Frequently, ragged edges were 
obtained and ribbon surfaces were rough. A ribbon of typical quality is 
shown in Figure 1. It was generally difficult to produce continuous 
long straight segments of ribbon, and although process parameters like 
wheel speed, distance of the nozzle from the wheel, and degree of 
superheat of the melt were varied, no particular combination of parame- 
ters were found that consistently produced smooth uniform ribbons. 


Figure 1. Macrograph of 
ribbon showing 
typical appearance 
of surfaces and 
edges. Roughness 
of as-spun Super 
Sendust ribbons 
could not be 
el iminated. 


Effect of Processing on Magnetic Properties 



The magnetic properties of the ribbon material were studied. Our 
efforts centered on developing processing schedules to m1nT~ J ze the coer- 
cive force of the material. As a target value, the lowest value of coer- 
cive force reported for the ternary Sendust composition is about If 1 mOe [4J. 

As melt-spun, the coercive force of the Supe - Sendust ribbons was 
approximately 1.9 Oe, which is comparable to the val*ie reported for as- 



spun Sendust ribbons prior to thermal annealing treatments. [1] Two types 
of processing treatments have been used in an attempt to dramatically 
reduct coercive force in the Super Sendust ribbons: thermal annealing 

treatments, and combined deformation and annealing treatments. 

A variety of thermal annealing treatments were used. Annealing tem- 
peratures, times and quench rates were varied. The lowest values of H c 
were obtained from ribbons annealed for 1 hour at 1200°C, (T^ temper- 
ature), cooling to 690°C at -10°C/min, holding at 690°C (T2 temperature) 
for 3 hours, and finally air cooling to ambient temperature. This treat- 
ment resulted in a reduction of H c to 0.5 Oe. The influence of temper- 
ature Ti and T2 on measured values of H c is shown in Figure 2. 



Figure 2. Coercive force cf Super Sendust ribbons after two-stage 

anneal ling treatments at temperatures Ti and T2. H c is a 
minimum for T2 slightly above the critical temperature T c 
for the B2+DO3 phase transition. 

Rolling of as -melt-spun ribbons was done to reduce suface roughness 
of the ribbons in an attempt to provide further reductions in H c . After 
thermal annealing of the ro.led ribbons, H c was still at the 0.5 Oe 
level. This compares unfavorably with the 0.01 Oe value obtained for 
Sendust ribbons [4]. 

Effect of Processing on Microstructure 


An order-disorder ohase transformation from the B2 ^CsCl) nigh tem- 
perature phase to the 003 (BIF3) low temperature phase was fc. J to occur 



at 680°C by electron microscopy examination of Super Sendust samples 
annealed over a range of temperatures from 570°C to 800°C. The primary 
defects in these phases are antiphase boundaries (AP3s). Because 
antiphase domain sizes which are of similar scale to magentlc domain wall 
widths (~ 0 . 1 ym) can interfere with magnetic domain wall motion, we 
studied the kinetics of antiphase domain coarsening. The measured kine- 
tics for thee Super Sendust alloy fit a coarsening model developed by 
Allen and Cahn [5]. A 4-hour annealing treatment at 670°C produced an 
antiphase domain size in the DO 3 phaie of about 1.4ym, which should be 
coarse enough to minimize this particular source of resistance to magne- 
tic domain wall motion in specimens with the DO 3 structures. 

The minimum coercivity measured in our alloys after heat treatment 
at 690°C thus corresponds to equilibration at a temperature just slightly 
above the B 2 >D 03 phase transformation, in the single-phase B2 phase 
field. Our results show that sample with the DO 3 structure have higher 
coercivity by comparison. 

Grain size of the as-spun ribbons was found to increase with 
increasing ribbon thickness. A 30 ym thick ribbon had a 1.5 ym grain 
size, while a 100 ym thick ribbon had a 10 ym grain size. Annealing at 
temperatures over 900°C produceu nearly complete grain growth (so that 
the grain size was limited by ribbon thickness) in one hour. Only those 
ribbons with the smallest grain sizes were found to have significant duc- 
tility. 
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Fluid Flow and Crvsta!iization Processes 

Heat Flow and Segregation in Directional Solidification 

Principal Investigator: August F. Witt 


RESEARCH ABSTRACT 


The thrust of this research program continues to be the optimization 
of heat transfer control In melt and solution growth of elemental and 
compound semiconductor systems. Focus Is placed on the development of an 
optimized crystal growth facility which will permit a meaningful explora- 
tion of the potential of reduced gravity for growth of electronic 
materials. Significant developments during the past year Include the 
fol lowing. 


(1) A computer controlled, heat pipe based Bridgman growth 
facility has been comp let ;d and Is currently under- 
going extensive thermal characterization. 

(2) A two-dimensional heat transfer analysis for Bridgman 
growth, cased on boundary conditions which are 
realizable in the constructed facility, has been 
completed and Is used for the optimization of critical 
gradient zone conf Iguraf Ion tn this fact I Ity; It Is 
also used as a basic criterion for crucible and 
ampoule design. 

(3) The heat transfer control concepts developed In this 
research program had a major Impact on Army and NRL 
sponsored research and have led to the submission ot « 
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patent disclosure (wl+h joint NASA-Ar my sponsorship) 
on a "heat exchange system for optimized control of 
Interdependent thermal parameters in crystals grown by 
LEC and the conventional Czcchralski -technique". The 
extension of the heat transfer control concepts devel- 
oped for Bridgman growth to Industrial Czcchralski 
growth is considered by many as a major breakthrough. 

(4) NRL sponsorship for the preparation of a flight exper- 
iment Involving CdTe has been ascertained. 

(5) It has been determined that, contrary to assertions In 
the STAMPS report, transverse magnetic fields of up to 
j kG applied to vertical Bridgman growth will not 
establish diffusion controlled segregation during 
growth; the effective distribution coefficient for 
growth In the presence of a magnetic field Is not 
significantly larger than that in the absence of a 
magnetic field. For growth systems with a diameter of 
more than 1 cm, magnetic fields are no substitute for 
reduced gravity environment-. 

(6) During the past year this research program enjoyed two 
supporting graduate fellowships from IBM and BTL. 

(7) The research conducted during the past year Is docu- 
mented In seven publications (one published, two In 
print and jr being prepared for publication), one 
Ph.D. thet. and one Masters th6sl$. 


This research program has now reached a phase where preparation of 
flight experiments for Implementation In the Shuttle program can and 
shou I d be I n 1 1 1 ated. In this context It Is also suggested to conduct a 
detailed analysis of the developed growth facility with regard fo Its 
suitability for flight operation. 
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(Basic elements of submitted patent disclosure 
on heat transfer control in Czochralski growth.) 

Heat Exchange System for Optimized Control 
of the Following Interdependent Parameters: 

• Crystal temperature (range 500 - 1100 °C) 

• Growth interface morphology 

• Temperature gradients about the growth interface 

• Thickness of encapsulating B 2 0 3 layer 

• Radial thermal symmetry in the vicinity of the 

growth interface 

• Radial thermal gradient in growing crystal 



® Inconel Heat Pipe v ith Operating Range 
from 500 °C - IlfO 3 C 
® Peltier Heating / Cooling System 
(D Thermal Insulat on 
(D Auxiliary Heat ng Element. 

© Axial Positior of Heat Exchange Sys'em (adjustable) 
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VERTICAL BRIDGMAN TYPE CRYSTAL GROWTH; 

A- MEAT. TRANSFER ANALYSIS 


INTRODUC TI ON 

Semiconductor crystal growth experiments carried out in a reduced 
gravity environment using a conventional Bridgman configuration revealed 
deficiencies In crystalline and chemical perfection which have been 
attributed primarily to inadequate heat transfer control during growth. 
Since this growth technique is of Importance for the growth of electronic 
materials on earth and, moreover. Is a viable technique for the explora- 
tion of the potential of growth In a reduced gravity environment, a 
comprehensive, NASA sponsored effort has been made for Its optimization. 
This work addresses the thermal characteristics of crystal growth by the 
vertical Bridgman method to (1) evaluate Its Inherent limitations and (2) 
to develop an understanding of the primary parametric behavior necessary 
for a cognizant appraisal of thermal design alternatives. The first part 
of the thermal analyses, concerned with the axial temperature distribu- 
tion of the charge, has been discussed in the 1981 Materials Processing 
Center Annual Report; a brief summary of these results is Included here. 
Presently, focus is placed on the radial temperature variations within 
the charge. A complwte account of both thermal analyses Is contained in 
the Ph.D. thesis of T. Jaslnski, Department of Mechanical Engineering, 
1983. 


The basis for the thermal models Is an experimental growth system In 
which well defined thermal boundary conditions are achieved through the 
use of coaxial heat pipes, A schematic of this configuration is shown In 
Fig. 1.1. The gradient control region located between the hot and cold 
heat pipes I s env I si oned as prov 1 d i ng a degree of flexibility for the 
thermal boundary conditions of the furnace near th6 growth interfaco. 
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The one-dimensional thermal model considers the charge and crucible 
to be a fin moving through the boundary condl rlons Imposed by the fur- 
nace. This model Is hereafter referred to as the moving fin model. The 
axial temperature distribution Is described by a second order, ordinary 
differential equation within each axial region of uniform furnace 
boundary conditions and charge thermal properties. Given the assumptions 
employed by the one-dlmenslonal model, there are four such regions. The 
solution for the axial temperature distribution is readily obtained by 
computer. Carrying out the analytical solution, however, yields rela- 
tively simple expressions for (1) the axial temperature gradient In the 
liquid at the growth Interface, G^, and (2) the axial position of the 
growth Interface, $j. These expressions Include the effects of the 
thermal coupling between the hot and cold zones and the charge, gradient 
zone length, thermal conductivity change and generation of latent heat at 
the growth Interface, charge lowering rate, and charge diameter. Cri- 
teria are also developed which Indicate under which conditions the 
affects of the lowering rate of the charge, generation of latent heat and 
finite length of the charge may be neglected. Finally, the effects of 
the crucible and radial thermal resistance within the charge <>re shown 
analytically to result In a decrease of the effective thermal coupling 
between the furnace and the charge. 

IH IE- THQ-D1MEKS1 Q MAL T HERMAL 

Thermal criteria for satisfactory Bridgman crystal growth Include 
requirements for the shape of the sol Idlf leaf Ion temper? e Isotherm. 
The moving fin model does not permit evaluation of the isotherm shape 
since It neglects radial temperature variations within the charge. For 
this purpose, a two-dimensional model of the Bridgman growth system Is 
developed In the present section. Results from Ihe two-dimensional mode! 
are used In Section 4 In order to establish the causes of radial tempera- 
ture variations within the charge and to study their dependence on the 
thermal parameters of the Bridgman growth system. 


94 


The Concentric Fin Approach 


ORIGINAL RAC . U 
OF POOR QUALITY 


The present two-dimensional thermal model of vertical Bridgman 
growth treats the entire system as a set of coaxial fins. Radial temper- 
ature variations are thus accounted for by differences In the fin temper- 
atures at the same axial location. The model can therefore be considered 
a hybrid between a continuous description In the uxlal direction and a 
finite difference description In the radial direction. This approach, 
rather than a complete finite difference or finite element formulation. 
Is taken for several reasons: 

(1) The primary objective of the two-dimensional model Is to 
develop Insight Into the origin and behavior of radial 
temperature gradients near the growth Interfac for the 
purpose of controlling the shape of this Interface 
through system design. The Increased accuracy achievable 
using a fine grid network In a finite difference or 
finite element model is nor considered essential for this 
objective and, furthermore, does not appear warranted for 
systems for which reliable thermal property data are no' 
yet available. The prosent model Incorporates all the 
critical thermal effects; It Is expected, therefore, to 
be able to predict the radial gradient behavior at least 
to first order approximations. (By utilizing enough 
coaxial fins, the present model should. In principle, be 
as accurate as a finite difference model.) 

(2) The present model Is easily Implemented on th3 computer 
and Is a simple extension of the one-d I mens tonal moving 
f I r. mode I . 

(3) Even with the simplifications Inherent to + he concentric 
fin approach, the added comp I ex It y of a two-dimensional 
formulation necessitates the use ?f a computer to calcu- 
late the two-dimensional temperature distribution of the 
charge. The concentric fin model can be simplified. 
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however, without less of pertinent Information, to a 
point which permits presentation of come Important 
results In analytical form. 


I ha Describing Equatio ns For The Con " antrlc Fin Modal 


The factors of concern for the development of the concentric fin 
mode! of the Bridgman growth system depicted In Fig. 1.1 are Indicated In 
Figs. 3.1 and 3.2. The two-dimensional thermal model Incorporates the 
following assumptions. 

(1/ Radial temperature variations within the charge are 
accounted for by modeling the charge as two radially 
d I s+ 1 net f I ns. The I nner fin s a cy I I ndwr of d I ameter 
D] n and has an axial temperature distribution T| n (Z). 

The outer fin Is an annulus of Inrer and outer diameters, 

D| n and D, and has an axial temperature distribution 
T ou t(Z). The diameter D| n is chosen so that the Inner 
and outer fins have equal cross sectional a rear The 
temperature of each charge fin Is not considered a func- 
tion of radius. The difference In temperature between 
T| n and T ou f Is 1 idlcatlvo o: radial temperature varia- 
tions within the charge. 

(2) The crucible Is modeled as a coaxial annular fin !n 
contact with the outer charge fin. The crucible fin 
temperature, T cr (Z), Is not :j function of radius. 

(3) As opposed to the '-jr responding assumption Ir, the moving 
fin model, +he gradient zone 1 s not sdlobatlc. This 
assumption Is rolexed In the two-dimensional mooel so 
that the effect of radlnl heat transfer between a non- 
perfect I y Insulating gradient zone and the charge cc r bo 
studied. 
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Fig. 3.1. Schematic of the charge, crucible, gradient 
zone annulus and furnace heat pipes in the 
concentric fin model. 
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Cut-away sketch of charge, crucible and 
gradient zone annulus showing the geometry 
of the concentric fins. 
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(4) The portion of the furnace In the gradient zone Is con- 
sidered to be an annulus of conductivity kg and cross- 
sectional area Ag located coaxially with tli3 charge and 
crucible. The temperature of the gradient zone annulus 
Is assumed constant In the radial direction; Its axial 
temperature distribution, Tq(Z), depends on Its thermal 
interaction with the charge and crucible. Temperature 
boundary conditions are used at the ends of the gradient 
zone annulus ana no heat loss Is assumed at Its outer 
d I ameter . 

(5) Hot and cold zone furnace temperatures, Tf •; anl Tf t g 
respectively, are uniform, reflecting the heat pipe 
behavior. The lengths of the heat pipes are assumed 
Infinite. 

(6) The system Is at all times In a quasl-steady state; ..e., 
transients are neglected. 

(7) No heat transfer by natural convection In the melt. 

(8) Heat exchange between the furnace and the outer crucible 
surface Is described by a heat transfer coefficient, h, 
which Is constant within each zone. 

(9) Thermal properties do not depend on temperature except 
for a difference In conductivity between the melt and 
crystal portions of the charge. 


Assumptions (5) through (9) are also used In the moving fin model. 

On the basis of the above assumptions, each of the concentric fins 
Is analogous to the moving thin rod treated by Cars I aw and Jaeger [l]: 


, , <J*T , <4T 



[3.1] 
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where: k ■ thermal conductivity of a fin 

A x ■ cross sectional area of a fin 
T ■ fin temperature 

Z ■ axial coordinate measured from the center of the 
gradient zone toward the cold zone 
V * lowering rate 
f * mass density of a fin 
Cp * specific heat of a fin 
q * heat transfer per unit length to the fin from Its 
surroundings 

In Appendix A, the q term of eq. [3.0 Is approximated resulting In the 
following set of coupled ordinary differential equations for the axial 
temperature distribution of the fins: 


dc ; 4 

^ e loc dC + A, ^ ^ 

[3.2aH 

d*fc 

ics 


[3.2b] 


^«c< jt (■A* A,- Ad 4 A*(A - Ad = o 

[3.2c] 

d*^ 

+ A s (.A, - Ail * o 

[3. 2d] 

where: £ * dimensionless axial coordinate, Z/D 



+ A») 

A ■ - A . ) 

?e„= VD(fC 
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subscript "loc" refers to local charge properties or 
zone heal transfer coefficient 
subscript ,s .'" -->fers to the crucible 

The coefficients A A 3 , A.-, and A 5 depone on the radial thermal 

resistance between the concen+ric t'ns. They are given by: 



8 


rr * 

j o v , \ Z i 


’ ' *oc 



A K k (£‘-^ 

*3- o 


'<* jloc 


[3.3a] 


[3.3b] 


[3.3c] 


[3.3d] 


[3.3e] 


where; g ^ "Dcr/15 

D c< - crucible diameter 
"D = charge diameter 

V’ 

MJV A 

The fin equations tor and have beer, added to and subtracted 

from one ano+her 1 n order +0 produce eos. [3.2sJ and [3.?b] for and ft 
respectively. Si ice the cross sectional aceas of the 1 nner and outer 
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charge fins are equal, <p Is interpreted as the mean cross sectional 
temperature of the charge; ^ Is a measure of the radial temperature 
difference within the charge. Within the hot and cold zones, the furnace 
temperature,^ , Is the temperature of the hot and cold heat pipes, 
respectively, and eq. [3.2d3 for the gradient zone annulus In not needed. 
Within the gradient zone, Is Identified with and eq. [3.2d3 has 

no Pec let term since the gradient zone annulus Is stationary. 

Solution Qf The Concentric Fin Equations 

The set of eqs. [3.23 has constant coefficients If a separate set of 
equations is employed for each axial region of uniform thermal properties 
and furnace boundary conditions. As In the moving fin model, four such 
axial regions are present In the concentric fin model: the hot and cold 

zones and the liquid and solid portions of the gradient zone. (It Is 
assumed that the growth Interface is In the gradient zone.) The hot and 
cold zones are therefore described by eqs. [3.2a3, [3.2b3 and [3.2c3; the 
two regions within the gradient zone Include, In addition, eq. [3.2d3, 
representing the gradient zone annulus. 

The solution of eqs. [3.23 within a particular axial region Is 
described In many texts (e.g., [2,33). A system of p second order 
differential equations can be changed to an equivalent system of 2p first 
order differential equations by the definition: 

<f>' a [3.43 

Performing this transformation and putting the result In matrix notation, 
eqs. [3.23 become: 

- PX « F C3 ' 5] 
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'< 


V 




*cr 


<*>cc 

4m 

i / 

V - 



“ 





and cold zone 

<< 





gradient zone 


(Variables In boldface denote matrices.) The elements of P and F are 
found by direct substitution. In the g’-.-'jient zone, F Is an empty 
matrix. 


The solution of eq. [3.5] Is of the form: 


X-c,E, 




■2.pW2 


e Wl ^ 


[3.6] 


The values of w are commonly called the eigenvalues of the matrix P and 
the values of E are commonly cal led the eigenvectors of P. The eigen- 
values and eigenvectors are easl ly found by standai d computer subrou- 
tines. The coefficients denoted by c In eq. [3.63 are the unknown con- 
stants of Integration. Xp denotes the particular solution of eq. [3.5]. 


Equation [3.6] Is written for each axial region. The hot and cold 
zones each contribute six unknown constants of integration; the liquid 
and solid parts of the gradient zone each contribute eight unknown con- 
stants of Integration. There are, therefore, 28 unknown constants of 
Integration which must be determined by using boundary conditions of (a) 
equality of temperature and (b) continuity of flux between corresponding 
fins of adjacent axial regions as well as (c) heat transfer coefficients 
at the ends of charges of finite length and (d) Imposed temperature 
boundary conditions at the ends of the gradient zone annulus. Applica- 
tion of the boundary conditions results in a set of 28 simultaneous 
linear algebraic equations whose unknowns are the required constants of 
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Integration. The solution of these simultaneous equations Is easily 
accomplished by standard computer subroutines. 

Approximation For Tiia Growth Interface Boundary Condition 

The solution of the fl.i equations requires a trial and error proce- 
dure when the extent of each axial region Is not Initial ly known. For 
example, the moving fin model shows that the Interface location, Qj, and 
the non-dimensional Interface temperature, ©i = <T|-Tf,c>/<Tf,H-T,.c>. 
are functionally related. When 0^ Is chosen as the Independent variable, 
must be found by Iteration. When £; Is taken as the Independent 
variable, the location of the Interface boundary conditions Is known 
a priori resulting In a direct solution for 6; as well as the entire 
axial temperature distribution. 

Figure 3.3 shows that the ax' a I position of the growth Interface In 
the Inner and outer charge fins Is, In general, not the same; t K ey are 
the same only when the Interface Is flat. The correct application of •’he 
concentric fin mode! requires that a separate set of eqs. [3.2] be 
written for the short region between the axial locations of the Inter- 
face. In this case, the Interface temperature, 9^, Is coupled with the 
Interface locations w Ithln both the Inner- and outer charge f Ins^i^n 
and respectively: 

d.7] 

Since this functional relationship Is not known a priori. It Is not 
possible to choose Initial values of and Ci 09 t which satisfy 

eq. [3.7]; l.e., fin boundaries at the Interface are not Initially known 
and therefore the solution requires an iterative procedure In order to 
locate these boundaries. The remainder of this section describes an 
approximate model of the Interface which precludes the use of a trial and 
error solution for the concentric fin model. 

The length of the region separating the Interface In the Inner and 
outer charge fins Is normally small since It Is desired to have an 
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Approximated Interface 
Location 



Location 


Fig. 3.3. Approximation for the interface 
boundary conditions used in the 
concentric fin model. 
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Interface of small curvature. The present model neglects this axial 
region and applies the boundary conditions of the Interface at the same 
axial location for each of the Inner and outer charge fins. The axial 
location of the Interface Is takon to be the location where the mean 
charge temperature, <f> m , Is equal to the Interface temperature, 0i (see 
Fig. 3.3). In this sense, the axial location where ^,-01 denotes an 
"average" Interface position (l.e., between the actual Interface posi- 
tions of the Inner and outer charge fins). The value of at this 
location, It;) , is Indicative of the curvature of the Interface. 

The use of the above approximation for the Interface boundary condi- 
tion permits the a priori choice of the "average" Interface location and, 
therefore, the solution for the temperature field In the charge Is not a 
trial and error procedure. The approximation Introduces no error If the 
thermal conductivities of the melt and crystal are equal and In the 
absence of the generation of latent heat, or If the interface Is flat. 
Otherwise, when the Interface curvature Is small, the error should like- 
wise be smal I and the predictions of the present model concerning the 
radial temperature variations In the charge are expected to be correct to 
a first-order approximation. 

The Interface flux boundary conditions, assuming that the melt/ 
crystal conductivity change and the generation of latent heat occur at 
the same axial location in each of the charge fins, are: 



[3.8a: 


[3.8b] 


[3 . 8cD 
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[3.8d] 


where: R K x k L /k s 

Sm ■= AWc.(Tf,«-*T^ t ) 

= latent heat of solidification 

(The subscripts and "$" refer +o the melt and crystal side of the 
Interface, resnec - 1 1 \ ■? I y.) Equation [3.8a] assumes that the solidifica- 
tion rate !:. equal +c +re 'ower Ing rate, V. 

SAQ1AL— T.E A r J .QHS. w 1 :e .QH£33E 

Radial tempe-a^u. e var iaf ions near the growth interface are shown In 
this section +c -esu t from severs 1 effects usually present In Bridgman 
crystal growth, E/s ferns with an adiabatic gradient zone, first studied 
by Chang aid '••'i.cov addresses only one of these factors. The other 

factors have not deer adequatley considered In the literature: non- 

perfec+!y 1 rssu ! e~ ' no ar^c’ent zone (hg J 0), and thermal conductivity 
charge and/or generz+'on ot latent heat at the growth Interface In the 
presence of a crucible. Computer solutions of the concentric fin model 
are utl i Ized to graph ical ly III ustrate the parametric behavior of the 
radial temper alure variations. Special attention ’ given, however, to 
the development of approximate analytical expressions which serve to 
emphasize the relative Importance of the various growth parameters on 
radial temperature variations. 

The variable defined by the concentric fin model of section 3 Is 
used to Indicate the radial temperature variations within the charge. 
Ti.e results tor <p A to be shown In this section are conveniently approxi- 
mated to the interface shape by the following relation developed In 
Appendix B: 


I L 


[ 4 . 1 ] 
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where N Is the radius of curvature of an Isotherm at a particular axial 
location In number of charge radii. (N Is negative for an Interface 
which Is concave toward the melt as shown In Fig. 1.1.) Accordingly, the 
degree of flatness of Isotherms Is enhanced by large axial temperature 
gradients (d and small radial temperature gradients ($ a ). 

A system defined by the following parameter values Is used as a 
reference system against which changes In Individual parameters can be 
compared: 

1. Equal heat transfer coefficients In the hot and cold 
zones. 

2. Axial convection of heat due to charge motion Is 
small; i.e., Pe Is neglljlble. This condition Is 
shown by the moving fin modei to be satisfied If: 

PeACfti**)''* "1 

where Bi** Is called the effective B ot number. BI** 
is related to the Blot number, BI, additionally 
accounting for the presence of the crucible and 
radial thermal resistance within the charge on the 
axial temperature distribution of the charge. 

3. No generation of latent heat; I.e,, Pe$Rn = 0. 

4. Equal thermal corduc Ivitles In the melt and crystal; 

I.e., R« = 1. With condition (1) above, Bl H = BIq. 

5. Equal charge lengths In e hot and cold zone; I.e., 

= where Xls the dimensionless axial length of 
a furnace zone, L/D. 

6. Temperature boundary conditions on the ends of the 
gradient zone annulus which satisfy ^(-X^/2^ + 

QclKI 2 ) = 1 - 
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A system which has the above parameter values Is referred to as a 
symmetric system. The temperature distributions of a symmetric system 
are of the following form: 


* 4w('0 * I 
<M 0 ♦ 

4>o>(C.) * <£<s(-C) = l 

* M-z) s ° 

Accordingly, It Is seen that It Is only necessary to consider one-half of 
a symmetric system, either C, < 0 or C, > 0. 

Causes Of Radial Temperatur e. Variations A : t hin The G rodlanl Zone 

Subtracting eq. t3.2b]] from eq. [3.2al, and assuming Inat the Pec let 
term Is negligible (one-dlmenslonal results Indicate that the Pec let term 
Is small for typical semiconductor crystal growth conditions), y e ! ds: 


A . dNja 
AC, X 


[4.2] 


The particular solution of eq. £4.23, denoted oy <p?, may be solved by 
using the appropriate Green’s function: 


Using eq. [4.3D» the particular solution of eq. [4.2] Is: 


L4.3] 




\s.z(c V 



£4.4] 
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where the value of A 2 given by eq. [3.3b] has been used, x Is a dummy 
variable of integration, and C> c and Ch ore the axial locations of the 
cold and hot ends of the charge, respectively. (Note that eq. [4.4] is 
not the complete solution for since ^ and are not Independent by 
virtue of the coupled nature of the fin equations. Nevertheless, the 
actual distributions for and must satisfy eq. [4.4].) 

JVu/atf is ncn-^ero when the combined axial heat conduction within 
(he Inner and outer charge fins is not constant with axial position; 
i.e., when there is heat addition to or removal from the charge. Equa- 
tion [4.4], therefore, clearly demonstrates that radial temperature 
gradients within the charge are caused by heat exchange between the 
charge and its surroundings. Such heat exchange occurs: 

(1) if. In the Idealized case of no crucible, there is b 
difference between the temperatures of the outer charge 
fin and the furnace, and the heat transfer coefficient 
between them Is not zero. 

(2) If, when a -.ucible is present, there Is a difference 
between the temperatures of th? crucible and the outer 
charge fin. 

Equation [4.4] expresses ^ at a particular axial location as the 
superposed effects of heat transfer to the change along Its entire 
length. However, the effect of heat transfer to the charge at one loca- 
tion, x, on the radial temperature gradient at another location, C, 
attenuates rapidly with the distance Detween the two locations, being 
proportional to the factor exp(-7.63|^- xj ). For example, when - x[ 
= 1, the va I ue of th i s exponent la I function Is 0.00049. Therefore, the 
radial variation In temperature at the location £ Is primarily affected 
by heat exchange with the charge within the range £±| (i.e., one 

charge diameter to either side of the given location). 
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Analysis of Radial Tem perature Variations In the .Grad 1 ant Zona 

The total solution for <£ A (C) Includes, In addition to eq. [ 4 . 4 ], the 
hoaiogenous part of the solution: 


<MO = C,e ' wc - Cte'" 5 ' - <p\ 

Equation [4.5] Is applied separately to each region of the charge of 
uniform properties; given the assumptions of ths concentric fin model, 
there are two such regions — the I Iquld and the sol Id portions of the 
charge. The coefficients Cj and Co In eq. [4.5] are found by applying 
boundary conditions at the ends of the charge and at the growth Inter- 
face. 


Radial temperature variations near the growth Interface are of 
primary concern. Since the Interface Is ossumed to be In the gradient 
rone, eq. [4.5] Is applied. In the remainder of this section only to the 
gradient zone; l.e., C < . It Is also assumed that there Is 
at least one diameter of charge In both the hot and cold zones: 

U > (x 6 /z - 0 C'- 6I >3 

As noted earlier, heat exchange with the charge at distances greater than 
about one diameter from the ends of the gradient zone contribute little 
to the radial temperature variations In the gradient zone. Equations 
[4.6] thus permit the limits of Integration In eq. [4.5] to be contracted 
If desired, and further, the effects of the boundary conditions at the 
ends of the charge on <t> & can be neglected. 

When the solid and liquid thermal properties of the charge are 
equal, l.e., R« 3 eq. [4.5] need be applied only once to the entire 
length of the charge. The homogenous part of the solution dlsv^pea’ s 
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under the restrictions of eq. [4.63 and the total solution for 0^ (c) 
within the gradient zone becomes: 




[4.73 


where: 


- Klz <C* \U 


When R« ^ 1, eq. [4.53 must be appi led separately to the I Iquld and 
solid portions of the charge using the following boundary conditions at 
the growth Interface: 




[4.8a3 


R 


K 



’ J <A,(0’ 

1 « 1 L 

dC ] 


[4.8b3 
(also eq. [3.8b3) 


Applying these boundary conditions yields: 


ML* LM d , + 1lW C4.n 

where: - X<, 1 2 < C, £ £■ 


, Is w c4.io] 
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The variables 1 1 _ (C) and i${5) represent the respective contributions of 
radial heat exchange with the liquid and solid portions of the charge to 
the value of ^ at any axial location C, within tho gradient zor-.e. They 
are defined as follows: 



[4. Ha] 


[4. 11b] 


The first term In each of eqs. [4.93 and [4.103 represents a contri- 
bution to Ich Is produced at the growth Interface and which 

attenuates with distance from the Interface. This contribution Is apsent 
when R« = 1 . 

The value of at the Interface, A te) , Is determined from either 
eq. [4.9} or eq. [4.103: 

4>t (C.) - ^7 [l,fc,) + Rk.Il («,•)! C4 -’ 2] 

An attract I ve possl b 1 I I tv tor the quant I tave utilization of eqs. 
[h. 93 and [4.103 Is to provide <J from + he results of one- 

dlmenslonal modeling. In this way, a two-dimensional temperature distri- 
bution Is obtained only through these equations and a one-d I mens Iona ! 
model. The error In calculated in this manner would be due to approx- 
imations in (1) the one-dlmenslonal model which estimates 
(2) the concentric fin model which provides eq. C4.21. 
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One-dimenslonal models do provide good approximations for the charge 
mean temperature distribution. Taking two derivatives of such a distri- 
bution In order to obtain a good approximation for Is «o r e 
I Ikely to be unrel table. Results presented later In this section. In 
fact, demonstrate that the presence of a crucible produces a distribution 
of radial heat exchange with the charge which Is not predicted by a one- 
dtmenslonal model; therefore, the use of cal culated from one- 
dimenslonal models Is not appropriate for use In eqs. C4.11D. Neverthe- 
less, even In these cases, the Green's function approach developed in 
this section may be used to obtain a qual Itatlve understanding of the 
behavior of radial temperature variations within the charge. 

Penetration Of Radial Temperature Variations Into The Gradient Zone 
From The Hot And Cold Zones 

Heat transfer to the charge In the hot zone and from the charge In 
the cold zone must occur so that the required axial temperature gradient 
at the growth interface Is obtained. Such heat exchange, according to 
the results of section 4.1, produces radial temperature variations 
with'* the charge which, even within an adiabatic gradient zone. 
Influence the shape of the growth Interface. 

Fu and Wilcox Q6] studied the penetration of radial temperature 
variations into an adiabatic gradient zone for systems without a crucible 
through the use of a two-dimensional finite difference model. Isotherm 
shape Indicate, for various B1 and Pe, that radial temperature varla- 
tk s within the charge In the hot and cold zones attenuate rapidly from 
rne ends of the gradient zone. If the gradient zone Is suff Iclenth’ 
long, there ex'sts a region within which the shape of Isotherms are quite 
flat. 


^sterns.. M-l thout a Crucible 


Figure 4.1 shows the variation of ^ In an 
obtained by the concentric fin model of section 
crucible. (Note that the vertical scale of 


adiabatic gradient zone, 
3 for systems without a 
Fig. 4.1 changes from 
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logarithmic to linear at ilO~* so that both negative and positive values 
of ^ may be plotted In the same figure.) Accordingly, It Is seen that 
attenuates at a rate approximately proportional to exp(-7.6 5*) where 
C*Is measured from the ends of the gradient zone. It Is also seen that 
4^ l n the gradient zone Increases with the Blot number of the hot and 
cold zones reflecting greater amounts of radial heat transfer to the 
charge In these zones. 

The variation of the radius of curvature of Isotherms for the same 
cases considered In Fig. 4.1 Is shown In Fig, 4.2. These curves were 
calculated using eq. [4.1] along with values of and eJ^/jC determined 
from the concentric fin model. Accordingly, there Is a central region of 
the gradient zone where the Isotherm shapes are quite flat. To obtain an 
interface which Is slightly concave toward the crystal, l.e., small but 
positive N"1, requires that the Interface be located slightly to tne hot 
side of The center of the gradient zone. Also seen In Fig. 4.2 Is that 
Isotherms are more curved as B! Increases. Therefore, for the parameter 
values considered. Increasing the Blot number of the hot and cold zones 
has a greater effect on Increasing than upon Increasing the axial 
temperature gradient of the charge In the gradient zone (c.f., eq. 
[4.1]). 

The effect of unequal hot and cold zone Biot numbers In an otherwise 
symmetric system Is shown In Fig. 4.3. As the hot zone Blot number 
Increases relative to the constant cold zone Blot number, 4^ becomes more 
positive In the hot portion of the gradient zone. In this region, there- 
fore, Isotherm shapes become more concave toward the crystal and the 
location of the flat Isotherm shifts slightly toward the cold zone. The 
effect of changing thermal conditions In the hot zone attenuates rapidly 
from the end of the hot zone and effects only slightly the radial 
temperature variations near the cold end of the gradient zone. 

Figure 4.4 shows the variation of In systems which are symmetric 
except R k i 1. In the systems considered, the cold zone Blot number 
remains constant} the hot zone Blot number therefor* changes with 
(B1 h = BIq/R^). As R« Increases, It Is seen that In the hot part of 
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hot end £ in gradient zone cold end 


Fig. 4.1. Variation of <t>^ in an adiabatic gradient zone 
without a crucible. The system is symmetric. 
116 






hot end C, in the gradient cold end 

zone 


Fig. 4.3. The effect of unequal hot and cold zone 
Biot numbers on <{>* in an adiabatic 
gradient zone. Tne system is otherwise 
symmetric. 



Fig. 4.4. Variation of in an adiabatic gradient 
zone without a" crucible. The system is 
symmetric except R« f 1 . 
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the gradient zone decreases whereas <P tn the cold part of the gradient 
zone Is little affected. This behavior Is similar to that of a changing 
value of BI h exhibited In Fig. 4.3. Further, a disturbance appears at 
the Interface which Increases as R K Increases and which dies out with 
distance from the Interface. This disturbance represents the first term 

In eqs. [4.93 and [4.10] where It Is shown that It appears only when 

* '■ 

An analytical expression for the variation of in an adiabatic 
gradient zone may be obtained by using the mean temperature distribution 

determined by the moving fin model In order to describe In eqs. 

[4.11]. In order to distinguish between the mean temperatures predicted 
by the moving fin and concentric fin models, 6U is used to denote the 
former while ^continues to denote the latter. The charge Is assumed 
Infinite in length; In this case, 0 m Is described by a simple exponential 
function In both the hot and cold zones. Further, the Peclet number Is 
assumed small enough to be negligible; In this case the argument for the 
exponential function Is given by 2(B I **) 1 

hot zone: 


cold zone: 


C » A Jz 


U) “ en.(A*/*)exp|> (KU-Z)] 


[4.13b] 


Substituting the second derivatives of eqs. [4.13] Into eqs. [4.1 1] and 
noting that = 0 In the adiabatic gradient zone, yields: 



sCLi-eu(-AtfcV\ 

TFfbtjpV'*] 


exp[-?.fc3(£+A a /2)] 


[4,14a 
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exp[-?.ijk,/2-S)] 


[4.14b] 


•here: - >*/2 f £ * X,/Z 


Equation [4.14a] represents the effect of radial heat transfer between 
the charge and the furnace In the hot zone on the radial temperature 
variations In the gradient zone. The effect of the hot zone dies out 
exponent (ally with distance from the hot end of rhe grad I ent zone at 
approximately the same rate observed In Fig. 4.1. Equation [4.14b] has 
an analagous Interpretation with regard to heat transfer In the cold 
zone. 


Substituting eqs. [4.14] Into eqs. [4.9] and [4.10] produces the 
desired approximation for In the gradient zone due to the generation 
of radial gradients In the furnace zones. Evaluating eqs. [4.14] at 
C = Ci and substituting the result Into eq. [4.12] yields ^ at the 
growth Interface. 

The one-dlmenslonal results required by eqs. [4.14] are the effec- 
tive Blot numbers and the mean charge temperatures at the hot and cold 
ends of the gradient zone. For a system with Big = 0 and which has a 
negligible Peclet number, the axial temperature gradient In the gradient 
zone, G, Is constant except for a step change at the growth Interface. 
The temperatures at the ends of the gradient zone are therefore given by 
the following expressions: 


e v - (.xj2*c,) 

©mt WO* 9l *<*.( AJz-O 
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where: Pe Is negligible 

B'6 ■ 0 

-K/i « C, s X 4 /2 

The non-dimensional Interface temperature, 0£ , G|_, and G s are given by 
the moving fin model. 

The circled points In Figs. 4.1, 4.3, and 4.4 indicate vc’ues of 
calculated through the use of eqs. C4.14], It Is seen t the values 

of obtained by *-he two methods agree very closely lr<dl tg that ths 

moving fin model provides a good approximation for . , for the 

systems considered. The largest difference occurs for the largest value 
of B! reflecting a decrease In accuracy of the effective Blot number as 
Bi increases. 

S ystems With a Crucible 

The effec+ of the presence of a crucible on radial temperature 
variations In an adiabatic gradient zone are shown In Figs. 4.5 and 4.6. 
The systems examined are symmetric and Infinite In length; Fig. 4.5 
cons I ders variable S (w 1 th Kg * 1 ) while Fig. 4.6 cons I ders variable Kg 
(with & = 1.5). The primary effects of the crucible are seen to be a 
reduction In the value of at toe ends of tne gradient zone end a 
reduction In the rate at which attenuates toward the center of the 
gradient zone. 

The reduction In the value of ^ a+ the ends of the gradient zone 
when a crucible Is present Is attributable to the effect of the crucible 

ft ft 

on the effective Blot number, BI . In most cases, the presence of a 
crucible decreases BI compared to a charge without a crucible. A 
reduction In the heat exchange between the charge and furnace Is there- 
fore expected, resulting In decreased radial temperature variations 
within the charge. The Individual curves of Figs. 4.5 and 4.6, labeled 

ftft 

with their respective values of Bi , demonstrate the correlation between 
and BI**. 
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hot end L in gradient zone cold end 


Fit], 4.5. Variation of <j>a in an adiejatic gradient 
zone with a crucible. The system is 
symmetric with variable diameter ratio f 
(K = 1.0). 
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hot end £ in gradient zone cold zone 

Fig. 4.6. Varv* ion of <j>a in an adiabatic gradient 
zone with a crucible. The system is 
symmetric with variable thermal 
conductivity K (6 = 1.5). 
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The reduci ion in the rate of attenuation of ^ when a crucible Is 
present results from radial heat exchange between the crucible and the 
charge near the ends of the gradient zone. Such heat transfer occurs 
when the crucible and the charge temperatures are not equal. A value of 
zero for the local Blot number Is no longer sufficient to eliminate heat 
transfer to the charge In the gradient zone. 

The variation of near ' t ' he h°+ encl +h® gradient zone. Is 
shown in Fig. 4.7 for the case of 5 = 2 from Fig. 4.5. It Is seen that 
the difference In the crucible and charge temperatures created by heat 
transfer In the hot zone furnace extends Into the adiabatic gradient 
zone. The resultant effect on radial heat exchange w I th the charge I s 
Indicated by the distribution of In Fig. 4.8 (lower solid 
curve). It Is seen that Is not zero In the adiabatic gradient 
zone; the effect of the crucible Is to extend radial heat transfer taking 
place In the hot zone Into the gradient zone. Analysis of eqs. [4.11] 
verifies that this extension of radial heat exchange serves to decrease 
the apparent rate of attenuation of ^ at the ends of the gradient zone. 

Assuming that the Pe term Is negligible, the one-dlmenslonal moving 
fin model predicts: 


AC 





[4.161) 


Since Big = 0, the moving fin model therefore predicts that 0 
In an adiabatic gradient zone, whether a crucible Is present or not. The 
dashed curve of Fig. 4.8 Is the distribution of determined from 
the moving fin model. The moving fin model Is unable to account for the 
particular nature of radial heat transfer between the charge and crucible 
caused by the abrupt change In furnace boundary conditions at the ends of 
the gradient zone. Although the error Is of minor consequence In calcu- 
lating the distribution of results from the moving fin model. In 
particular are Inadequate to describe the local variation of 
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Fig. 4.7. The distribution of 4> cr - <j> m near the hot and 
cold end of an adiabatic gradient zone, 
cal elated by the concentric fin model. 



c 


Fig. 4.8. Distributions of <}> m and d^/d^ near 
the hot end of the gradient zone, 
calculated by the moving fin and 
concentric fin models. 
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radial heat transfer near the ends of an adiabatic gradient zone when a 
crucible is present. 

Table 4.1 compares values of 4^ at the end of the hot zone (i.e., 
-0.5) calculated by eq C4.14a3 (I.e., with the use of the moving fin 
model) and by the concentric fin model for the cases shown In Figs. 4.5 
and 4.6. It appears that eq. Q4.14a3 predicts well the quallta+’''e 
trends as K and change as well as oelr.g, to a certain degree, quantita- 
tively accurate. This conclusion provides furr.’<&r evidence the 
correlation between BI** and notices ! n Figs. 4.5 and 4.6. 

Such a favorable comparison may seem, however, surprising In view of the 
limitations In the moving fin model discussed It. the pr sceedlng para- 
graph. This apparent anomoly Is explained by a closer examination of the 
distributions of and shown In Fig. 4.8. 

The abrupt change In predicted by the moving fin model at 

the end of the gradient zone becomes smoothed out In The distribution of 
given by the concentric fin model. The moving fin model over- 
predlcts the heat transfer to the charge In the hot zone while under- 
predicting the heat transfer In the gradient zone. The overpredicted 
heat transfer In the hot zone contributes to while the lack of 

heat transfer In the gradient zone reduces j If the effects of 

the two errors In /<;)£*■ approximately balance, the moving fin model 
provides a good estimation of • The balance of the two errors 

Is Illustrated by the dotted curve in Fig. 4.8 which Is plotted by 
reflecting the section of the curve for labeled "a" about the 

axial location Q - -0.5 and adding to the section of the same 

curve labeled M b". If the distribution of ca I cu I ated by the 

moving fin model were to exactly match the dotted curve, the value of 
hl-Kti) from eq. [_4. 1 4aU would equal that given by the concentric fir. 
model. Figure 4.8 shows that these curves, while not exactly coinciding, 
are nevertheless close. Therefore, the use of the moving fin model 
approximation for can be expected to provide a reasonable esti- 
mate of at the end of an adiabatic gradient Tone. Table 4.1 Indi- 
cates that this conclusions is valid af least over the range of K (oc and £ 
exam I ned . 
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TABLE 4.1 


Comparison of eq. [4.14a] with the concentric fin model In 
predicting^ at the end of the gradient zone. The systems 
are symmetric with an adiabatic gradient zone. 




moving fin model 

. . 


K loc 

$ 

BI** 



eq. [4.14a] 

concentric 
fin model 

^ ...... 

1.0 

1 .001 

0.471 

-0.407 

0.703 

0.00407 

0.00407 

1.0 

1.5 

0.291 

-0.350 

0.675 

0.00284 

0.00263 

1.0 

2.0 

0.205 

-0.312 

0.656 

0.00216 

0.00193 

1.0 

5.0 

0.0619 

-0.199 

0.600 

0.000797 

0.000683 

0.25 

1.5 

0.357 

-0.374 

0.68? 

0.00331 

0.00318 

0.50 

1.5 

0.356 

-0.374 

0.687 

0.0C331 

0.00310 

2.0 

1.5 

0.200 

-0.309 

0.654 

0.00212 

0.00197 

4.0 



1,5 



0.121 

-0.258 

0.629 



0.00141 

0.00132 


Bi * B Iq = B! h = 0.5 
Symmetric system 

V K0 

X tt = X c = CD 
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The axial temperature gradient In the charge, , must undergo 

an abrupt change at the growth Interface when Pe$Rn 4 0 or when R« 4 1. 
(c.f., eq t3.8a]). The axial gradient In the crucible, on the other 
hand, does not change at the growth Interface. Figure 4.9 shows that 
such conditions produce a difference between the charge and crucible 
temperatures at the Interface which diminishes with distance from the 
Interface. The exchange of heat between the charge and the crucible 
which results from this temperature difference has a profound effect on 
the distribution of In the vicinity of the Interface. 

Figure 4.10 shows the variation of In an adiabatic gradient zone 
for systems which are symmetric except R« 4 I. The presence of a 
crucible Is Indicated by a value of 1.5 for the diameter ratio $. As R« 
Increases, It Is seen that the distribution of 4> k , and more Importantly, 
the value of <p k at the Interface, becomes more negative, causing Isotherm 
shapes to be curved in the adverse direction (concave toward the melt). 
Even a small difference In melt and crystal thermal conductivities (e.g., 
R K = 1.1 In Fig. 4.1 0) has a pronounced effect on <£ & (q;) . 

Values of R« greater than unity, being typical of semiconductors, 
hav6 been considered In Fig. 4.10. The variation of <f> k which occurs when 
R« Is less than unity Is similar to that shown In Fig. 4.10 except that 
the curves would be displaced toward more positive values of compared 
to the symmetric case of = 1. 

The effect of the generation of latent heat at the growth Interface 
on 4^ ,n +he gradient zone Is shown In Fig. 4.11. The systems considered 
are symmetric except that PesRn 4 0. It is seen that as Rh i ..creases, 
becomes more negative, slml lar to the effect of R« > 1 shown In 
Fig. 4.10. 

Although the curves of Fig. 4.11 are plotted with Rh as a parameter. 
It must be noted that the effect of the generation of latent heat Is 
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r k ,m 

-or- 

Pe s R H * 0 




Fig. 4.10. Radial temperature variations generated 
at the i .erface when R K ^ 1 . 
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proportional to the product PegR^i* The curves of Fig. 4.1) may also be 
Interpreted as Indicating the effects of changing Pe with a constant 
value of Rh> The value of Pes = 0.05 used in Fig. 4.1 I has been chosen 
so that PegRn would be sufficiently large to have an effect on for 
typical values of Rh» 

The generation of raalal gradients at the Interface Is due to the 
presence of a crucible. Figure 4.12 shows the effect of the crucible 
diameter ratio, £, on for systems which are symmetric except 
R* * 2. Accordingly, very thin crucibles are required to eliminate the 
Interface effect. 

Figure 4.13 shows the effect of changing the location of the Inter- 
face on the distribution of dfc In the adiabatic gradient zone for systems 
which are symmetric except that R« = 2. As the interface approaches the 
hot zone, the value ot <fr t at the Interface becomes more positive due to 
the Influence of radial temperature variations created In the hot zone 
(c.f., section 4.2). For the systems considered In Fig. 4.13, however, 
the dominant factor Is clearly the conductivity change at the Interface. 
In general, the Importance of radial temperature variations generated at 
the interface depends on their relative magnitude compared to radial 
temperature variations created by other effects. 

The Two-Fin. Thermal Mode.1 

In order to obtain from the results of section 4.1 a useful analyti- 
cal relation which expresses the primary effects of the therma 1 parame- 
ters on the radial temperature variations near the growth Interface, a 
simple and accurate analytical expression for Is necessary. 
Radial heat exchange between the charge and crucible In the vicinity of 
the interface Is not accounted for In the one-d i mens Iona 1 moving fin 
model; approximating d^Vn/dt* by d*9,*/d£*, therefore. Is not accurate 
near the Interface when a crucible Is present. The concentric fin model 
Is too complex to yield a simple analytical solution for . This 
section describes a thermal model of the heat transfer near the Interface 
which provides a simple approximation ford^/dC*. 
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hot end 


£ in gradient zone 


Fig. 4.12. The effect of the thickness of the 
crucible on the radial tenperature 
variations generated at the 
interface. 
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The new thermal model Is equivalent to the concentric fin model 
described In section 3 except for the following: 

(1) The charge Is modeled as a single fin as shown In Mg. 

4.14. Since there are only two concentric fins — the 
charge and the crucible — the present thermal model Is 
called the "two-fln" thermal model. The reprs;*»ntatlve 
radial location for the charge temperature Is taken as 
D| n /2, as first defined In Fig. 3.2, since this location 
divides the charge equd'ly with regard to Its cross 
sectional rea. The representative radial location for 
the crucible temperature Is again D 3/2 (c.t., Fig 3.2). 

(2) The Peclet number Is assumed negligible. 

(3) The gradlenh zone Is assumed adiabatic; hence the 
gradient zone amulus fin Is neglected. 

(4) Only the region In the vicinity >~- c ^he Inferface Is 

modeled. This region extends : a Interface, In 

either direction, at least the distance required for th« 
Interface effect to significantly die out. The effect of 
the remainder of the growth system Is Included only 
through Imposing the average charge/cruclblo axial 
temperature gradients In the gradient zene, G|_ and G$, In 
the liquid and solid respectively. These parameters are 
obtained from the one-d Imens Iona I moving fin model. 

« 

Radial temperature gradients generated by other factors 
are neglected. 

In accordance with the assumptions listed above, the fin equations 
for the charge and crucible are: 



T^-Tc 


cr 


'mtrjioc 



[4.17] 
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The thermal resistance, Rmcr,loc* defining the thermal coupling between 
the crucible and the charge, Is given by 1 4 J j 

, . -ML . L lpl c.,,3 

" ,t< > nc 2 . ir k c r 

Oeflne the following non-d I mans Iona I temperatures: 


3 A ky A y &C r 
^\©t A "* key Atr 


[4.20a] 


**■ v A*l 


[ 4. ?0bJ 


Is the average charge/crucible temperature weighted by their respective 
axial 1 her ma I conductances. Using eqs. [.4.19] and L4.20] In egs. [4.17] 
and [4.18], and non-d I mens Iona I I j I ng, results In the following form for 
the two-fln model fin equations: 




G 0 


L4.21] 



[4.22] 


[4.25] 


fcquatlon 14.21 ] Indicates that varies linearly In the gradient 
rone under the present assumptions. The axial gradient of G|_ and G$ 
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In the liquid and solid respectively, are assumed known from the results 
of the moving fin model. 

The general solution to eq. [4*22] is: 

In the I Iquld: Q < 

In the solid: £ > £ 

* Cisexp C-W 5 C) + C'is e *P ( w s6) 

where the coefficients denoted by C are the unknown constants of 
Integration. 

The assumptions of the two-fin model require that as 

C-*®; therefore, Cjl ■ C 2 S = 0. Continuity of temperature at the 
Interface requires u = ; therefore C 2 L = Cjs> Finally, 

continuity of flux in the crucible fin at the Interface requires 
k^r/dOi. = (A 4 > cf » * h,s boundary condition yields: 


[4.24a] 


[4.24b] 




zu 




Ws ^kWu 

I + K’ s (£ i -0 + ** + Ks(S l -0 


[4.25] 


The difference between and G$ Is, in dimensionless form: 


r r -a ~ (Rw. ~0 

UK S U 4 -I) 


[4.26] 


Substituting these results Into eqs. [4.24] produces the solution for 
r m, cr • 
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) L * C exp [- v t ( C-, - C ) ] 

C t^crX 58 C exp [- v s (C - C •, ) J 

where: Pg t R H ~ . 

, Rk^i » + »vsU 1 -^ 

l + K.U 1 -!) R**K s (SM [4.28] 

The constant C defined by eq. [4.28] Includes In the numerator the 
superposed effects of the generation of latent heat and change In charge 
thermal conductivity at the Interface. The axial gradient In the liquid, 
Gj_, Is always negative; If R* > 1, as Is typical of semiconductors, these 
two effects reinforce each other, causing a larger disturbance at the 
Interface than would occur by either effect acting alone. On the other 
hand. If R« < 1, 1+ may be possible to eliminate the Interface 

disturbance by choosing a value for the product Pe$RH that would make the 
numerator of eq. [4.28] equal to zero. This could be accomplished by 
changing either the lowering rate (thus changing Pe$) or the temperature 
difference between the hot and cold zone furnaces (thus changing R^). 
Either possibility, however, has potentially harmful side effects on the 
constitutional supercooling requirement. 

Radial Tempe rature Var ! at Ions Near.- the. Growth I nterface 


[4.27a] 

[4.27b] 


Substituting eqs. [4.27] Into the non-dimensional form of eq. [4.18] 
yields an expression for due to the heat exchange between the 
charge and the crucible near the Interface: 


Wu 

<U“ ' R.*K s (S*-0 


Cexp[-v L U-C;)} 


[4.29a] 



I + Kj(S‘ -0 


C e*p[-w s (C-C;)] 
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Using eqs. [4.29] In eqs [4.11] and [4.12] produces the value of <p. at 
the Interface: 



C Ks(S*-0 

J'.b'b Rv,-+\ 


l R^Ks(S*-0 


w 

I + k 4 (s*-o 


[4.30] 


The value of at the Interface determined by eq. [4.30] Is 
compared to the results shown previously In Figs. 4.10 to 4.13 In Table 
4.2. In general, the comparison Is favorable. Indicating that the two- 
fln model provides a reasonable approximation -or * In the 

vicinity of the Interface. The largest differences occur for those cases 
which exhibit effects not Included In the two-fln model. For example, 
the cases considered In rig. 4.11 use a larger Peclet number, 0.05, 
compared to that used for the other figures, 0.001, In order to 
accentuate the liberation of latent heat. The effect of the larger 
Peclet number, not Included In the two-fln model. Is to create a more 

i 

negative value of <#> A (£,). 

A smal I error In eq. [4.30] compared to results from the concentric 
fin model is also apparent for the cases of small & considered In 
Fig. 4.12. This error Is attributable to the neglect, by the two-fin 
model, of the radial gradients generated In the furnace zones. In order 
to approximate the attenuation effect, the value of determined by 

eqs. [4.9], [4.10], and [4.14] Is simply added tc that determined by eq. 
[4.30]. The results of this approximation, shown In the last column of 
Table 4.2, are seen to Improve the estimation of for small values 

of &. As S Incraases, the contribution to ^(Ci) due to the attenuation 
effect becomes less Important. 

Table 4.2 also demonstrates that eq, [4.30] Is unable to account for 
the effects of the changing Interface iocatlon on . This error Is 
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mLL /L2 

Comparison of the approximate methods of calculating ^(Q;) 
to results of the concentric fin model for the cases 
considered In Figs. 4.10 to 4.13. 




— 




variab'e 

concentric 


1 ■” 

eq. [4.303 plus 
attenuation effect 

Figure 

parameter 

fin model 

eq. [4.303 

from eq. [4.93, 
[4.10], and [4.14] 

4.10 

R k = 1.0 

.000000 

.000000 

.000000 


1.1 

-.000493 

-.000464 

-.000464 


1.25 

-.00110 

-.00114 

-.00114 


1.5 

-.00186 

-.00176 

-.00175 


2.0 

-.00274 

-.00261 

-.00260 


4.0 

-.00343 

-.00330 

-.00329 

4.11 

Ru = 0*0 

-.000023 

.000000 

.000000 


0.25 

-.000219 

-.000179 

-.000179 


0.5 

-.000416 

-.000358 

-.000360 


1.0 

-.000810 

-.000716 

-.000720 


2.0 

-.00160 

-.00143 

-.00144 

4.12 

& = 1 .0001 

-.000021 

-.000002 

-.000021 


1 .001 

-.000040 

-.000021 

-.000039 


1 .01 

-.000214 

-.C00 183 

-.000201 


1.05 

-.000818 

-.000730 

-.000744 


1.1 

-.00136 

-.00122 

-.00123 


1 .25 

-.00227 

-.00207 

-.00207 


1.5 

-.00274 

-.00261 

-.00260 

4.13 

A ; = 0.0 

-.00274 

-.00261 

-.00260 


-0.1 

-.00245 

-.00258 

-.00248 


- 0.2 

-.00210 

-.00254 

-.00231 


-0.3 

-.00166 

-.0025? 

-.00200 


i -0.4 

-.00108 

-.00249 

-.00138 


The values of G|_ and 6i required for use In 
eqs. [4.143 and [4.30] are I isted In Tables 4.3 and 4.4. 
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again attributable to the relative Importance at the Interface of radial 
gradients generated In the hot and cold zones compared to those generated 
at the Interface. As the Interface approaches the end of the gradient 
zone, the Influence of radial gradients generated In the adjacent furnace 
zone becomes Increasingly Important. 

The last column of Table 4.2 provides an approximation for the 
effect of changing Interface location on • The correction provided 

by these values to the results calculaled by eq. [4.303 (I.e., the second 
column of Table 4.2) assumes that the radial gradients generated In the 
hot and cold zones attenuate at a rate proportional to exp(-7.63£ ). The 
actual rate of attenuation Is, however, smaller when a crucible Is 
present (c.f., section 4.2). It Is expected that use of the actual rate 
of attenuation would produce a better correction to eq. [4.30] for the 
effect of changing Interface location. Unfortunately, a simple 
analytical expression for the attenuation of from the ends of the 
gradient zone In the presence of a crucible Is not available. 

Correction to the Axial Temperature Gradient In the Liquid 

At The Growth Interface : 

The moving fin model predicts the value of G|_ = for the 

axial temperature gradient In the liquid at the Interface. The 
constitutional supercooling requirement, however, depends on the value of 
at the I nterf ace. The d I f ference between these temperature 
gradients Is approximated from the results of the two-fin model. 
Differentiating eqs. [4.203 with respect to £, and solving fo r W^/dC) L 
yields: 


<K 


JL 


r + K&(S 2 -0 1 

L JC K 


[4.313 


Differentiating eq. [4.27a3 with respect to £ and substituting Into 
eq. [4.31] yields: 
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[4.32] 



KtCSM 
Rk ■*■ Ks (s l -0 


w L C 


The second term on the right side of eq. [4.32] Is a correction to be 
applied to Gj_ In order to account tor the effect of radial heat transfer 
near the Interface on the axial gradient In the charge at the Interface. 


The axial temperature gradient In the liquid at the Interface 
calculated by the moving fin model, Gl, and by eq. [4.32] are compared to 
results from the concentric fin model In Table 4.3 for the cases 
considered In Figs. 4.10 to 4.13. It Is seen that eq. [4.32] provides an 
excellent approximation to pfuKdAKL even though the correction term 
may be significant (e.g., Arge R« or large R^). The approximation Is 
most In error as the Interface approaches the hot zone (c.f., £; * ~0.3, 
-0.4 In Fig. 4.13) due to the additional Influence of radial heat 
transfer In the hot zone. 


Correction To the Interface Temperature : 


The moving fin model predicts for the Interface temperature the 
average charge/cruc I b I e temperature; I.e., Is Interpreted as 
The difference between ^,(Cj)and &i found by the two-fln model may be 
used as a simple correction to : 

0t - <MO - [<? (C) - 4UO] 

Rearranging eqs. [4.20] yields: 


\$(0- 'UC;)] L => 


Rtf + ^4 (^“O 





1+ K* (6 x -0 ^ c/ J 


[4.34a] 


[4.34b] 


S 1 nee and 4-fa) are equal In the liquid and solid at the 
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TABLE 4.3 

Comparison of the approximate methods of calculating d 
to results of the concentric fin model for the cases 
considered In Figs. 4.50 to 4.13. 



axial temperature gradient In the 
liquid at the growth Interface 




moving fin 



variable 

concentric 

model 


Figure 

parameter 

f I n mode 1 

g l 

eq. L4.32] 

4.10 

R* = 1.0 

-.351 

-.350 

-.350 


1.1 

-.332 

-.340 

-.332 


1.25 

-.308 

-.326 

-.307 


1.5 

-.275 

-.306 

-.275 


2.0 

-.228 

-.272 

-.228 


4.0 

-.140 

-.190 

-.140 

4.11 

Ru = 0*0 

-.351 

-.350 

-.350 


0.25 

-.344 

-.348 

-.344 


0.5 

-.338 

-.345 

-.338 


1.0 

-.325 

-.339 

-.325 


2.0 

-.300 

-.328 

-.300 

4.12 

S = 1.0001 

-.251 

-.252 

-.252 


1.001 

-.251 

-..252 

-.252 


1 .01 

-.251 

-.252 

-.251 


1.05 

-.249 

-.256 

-.249 


1.1 

-.247 

-.260 

-.246 


1 .25 

-.240 

-.267 

-.239 


1.5 

-.228 

-.272 

-.228 

4.13 

C; = 0.0 

-.228 

-.272 

-.228 


-0.' 

-.225 

-.269 

-.225 


-0.2 

-.222 

-.266 

-.223 


-0.3 

-.217 

-.262 

-.220 


-0.4 

-.211 

-.260 

-.217 
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' quality 


Interface, eqs. [4.34] imply that 4 l that Is, the two-fln 
model predicts a step change In £ at the Interface. This Is due to the 
definition of $ which averages the charge and crucible temperatures 
according to their respective axial thermal conductances. 

In order to approximate the effect of radial heat transfer near the 
Interface on the Interface temperature, the cor-ectlon term within 
brackets In eq. [4.33] will be taken simply as the arithmetic average of 
eqs. [4.34a] and [4.34b]. Using eq. [4.27] evaluated at C * Q, yields: 




KsU*-0 


[4.35] 


The Interface temperature determined by eq. [4.35] Is compared to 
that found from the concentric and moving fin models In Table 4.4 for the 
cases considered In Figs. 4.10 to 4.13. It Is seen that the moving fin 
model, provides a reasonable estimation to the Interface temperature, , 
without the correction afforded by eq. [4.35]; this Indicates that 
radial heat transfer at the Interface between the charge and the crucible 
has a relatively small effect on the interface temperature or location. 
Utl tiling the corresponding values of Gl listed In Table 4.3, the moving 
fin model would mlslocate the Interface by at most 0.1 charge diameter 
for the cases considered (c.f., R« x 4 from Fig. 4.10). Results calcu- 
lated from eq. [4.35], however, do compare yet more favorably to those of 
the concentric fin model. 

Dlahfltic- Gradient Zona 

Two additional parameters are required by the concentric fin model 
when heat transfer In the gradient zone Is considered: Big and Rg. The 
Blot number In the gradient zone Indicates the hermal coupling between 
the gradient zone annulus and the charge. Rg Is the ratio of the axial 
heat conducting capacities of the gradient zone annulus and the charge 
(c.f., eqs. [3.3]). A diabetic gradient -.one approaches abladatlc condi- 
tions as either Big or Rg approach zoro. Otherwise, the temperature 
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TABLE ...4^ 

Comparison of the approximate methods of calculating 
to results of the concentric fin model for he cases 
considered In Figs. 4.10 to 4.13. 




temperature of the growth 



1 nterf ace 


variable 

concentr I c 

moving fin 


Figure 

parameter 

fin model 

model 

©L 

eq. [4.35] 

4.10 

R* = 1.0 

.500 

.500 

.500 


1.1 

.505 

.507 

.509 


1.25 

.522 

.518 

.522 


1.5 

.540 

.534 

.540 


2.0 

.570 

.561 

.572 


4.0 

.648 

.634 

.651 

4.11 

HHBf M 

.517 

.500 

.500 


wmSa? 

.522 

.504 

.505 



.527 

.508 

.509 


1.0 

.536 

.516 

.519 


2.0 

.554 

.532 

.537 

4.12 

OO 

N 

• 

O 

o 

o 

.619 

.619 

.619 


1.001 

.619 

.618 

.618 


1.01 

.617 

.617 

.617 


1.05 

.610 

.610 

.611 


1.1 

.603 

.602 

.604 


1.25 

.587 

.584 

.588 


1.5 

.570 

.561 

.572 

4.13 

5 : - 0.0 

.570 

.561 

.572 


-0.1 

.602 

.593 

.604 


- 0.2 

.633 

.624 

.635 


-0.3 

.664 

.655 

.665 


-0.4 

.693 

.685 

.695 


The values of G|_ required for use In eqs. C4.353 are 
1 1 sted I n Tab I es 4.3. 
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d I f ference between the gradient zone annu I us and charge ( ) pro- 

duces heat transfer which affects axial temperature gradient's and also, 
as shown In this section, radial temperature gradients within the charge. 

Results From the Concentr ic Fin Modal 

Figures 4.15a, 4.15b and 4.15c show the variation of In a 
dlabatlc gradient zone as calculated by the concentric fin model. Tha 
systems In these figures are symmetric and do not consider a crucible. 
Since the systems, are symmetric, the curves have been plotted only for 
the hot side of the gradient zone. Each figure considers a constant Blot 
number throughout the furnace while the separate curves In each figure 
consider various values of Rq. The graphs Inserted In each figure are 
axial temperature distributions of ^ and <p m corresponding to the Indi- 
vidual curves of As a reference, the distribution of ^ In an 

adiabatic gradient zone (I.e., BIq a 0) Is also Included In each figure. 

According to Ftgs. 4.15, the value of Rq significantly affects the 
distribution of ^ In the gradient zone. When Rq Is small, heat exchange 
between the charge and the gradient zone annulus occurs only near the 
ends of the gradient zone; toward the center of the gradient zone Is 
nearly equal to v m . In such cases, the distribution of <P t Is only 
slightly altered from that In the corresponding system with an adiabatic 
gradient zone. As Rq Increases, the temperature difference - <p^ Is 
more pronounced Increasing the radial heat transfer to the cha r ge; as a 
consequence, > Increases as well. For sufficiently large values of Rq, 
is approx I mate I y linear; further Increases In Rq do not affect 
either or but serve only to conduct more heat from the hot to 

the cold zone furnaces through the gradient zone annulus. 

Figure 4.16 demonstrates the effect of Increasing BIq for systems 
which have a linear variation of \ as calculated by the concentric fin 
model. A I Inear variation of has been obtained by choosing a large 
value of Rq (I.e., Rq*1000). Increasing BIq has a si Ight effect on the 
axial temperature distribution In the gradient zone, tending to reduce 
the driving force for heat transfer In the gradient zone, . This 
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Bi^» Bi^ *Bi Q * 0.05 

No crucible ORIGINAL PAGE 19 

Symmetric system P00R Q UAL,TY 

V 1.0 

X H’ X C* 00 



hot end center 


l in gradient 
zone 

Fig. 4.15a. <j> m , and <j>g in a dia h ..tic gradient 

zone. Bi * 0.05, no crucible and 
symmetric system. 
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Symmetric system 

V'-o 


X U *X„* 00 



adiabatic 


BU«0) 


C in gradient 
zone 


Hn. 4.15c. 4 > a , 4 >, n , and in a diabatic gradient 

zone. Bi = 5.0, no crucible and 
symmetric system. 
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No crucibla 
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Symmttric tyatam 


linoor : Rq ■ I0 3 


V'° 


Xu* X r ■ 00 



V 10 ' 


BL-0.5 


BU« 1.0 



cantor 


£ In gradiant 
zona 


Bi *5.0 

6 


Fiti. 4.16. 4 > a , 4> m and in a diabatic oradient 

zone for increasinq Bi^. No 
crucible, symmetric system and 
linear ^n. 
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effect, however. Is more than offset by the augmented thermal coupling 
between the charge and gradient zone annulus represented by the larger 
Big. As shown In Fig. 4.16, the result Is that ^ Increases with Big. 
The curve for Big = 5 shows, for example, that ^ Increases, rather than 
attenuates, from the end of the gradient zone when the gradient zone Blot 
number Is greater than those In the hot and cold zones. 

Figure 4.17 shows that the effect on ^ of heat transfer In the 
gradient zone when a crucible Is present Is similar to those systems 
examined above without a crucible. (The system considered In Fig. 4.17 
Is analogous to that of Fig. 4.15b except for the presence of a cru- 
cible.) 

Analytical Expression For the Effect flJLa Dla hatiC 

Gradient Zone On Radial Tempera ture Variations 

The development of an analytical expression for In a dlabatlc 

gradient zone, through the use of the Green's function approach discussed 
In section 4.1, requires a simple analytical expression for 
within the d I abat Ic gradient zone as we 1 1 as w Ith In the furnace zones. 
The present model assumes that the primary effect of heat transfer 
between the gradient zone annulus and the charge Is attributable to 
thermal conditions near the axial location of Interest. This assumption 
Is qualitatively Justified by the high rate of attenuation of the effects 
of radial heat transfer distant from the location of Interest; Indicated 
by eq. [4.43. Accordingly, a linearized description for dV^/di* about 
the location of Interest,^, Is employed: 



- — 


dx* ~ dt 1 

lx*£ <d X 

L 


U-c) 




[4.36] 


Substituting eq. [4.36] Into eq. [4.73 yields: 



b<o.l 




H4 


[4.373 
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Next, eq. [4.16] Is used to approximate JVjdc* i th*s approximation 
assumes that of the moving fin model accurately represents 

d'tfu/JC* of the concentric fin model and that !>£ varies 

linearly near the location of interest. These assumptions are not always 
vcl Id since: 

(1) B!*^ may abruptly change at the ends of the gradient 
zon9 or at the Interface. 

(2) may not vary approximately linearly at the ends 
of The gradient zone, at the Interface, or when Rq Is 
smal I . 

(3) When a crucible Is present, d*"®. at the 
ends of the gradient zone. 

Consequently, the expressions which are obtained through the present 
model are suggested for use under the following constraints: (1) quali- 

tatively, to Indicate the thermal parameters most affecting radial 
temperature gradients caused by a diabetic gradient zone, and (2) quanti- 
tatively, as order of maglnltude estimations. 

Subsltu+ing eq. [4.16] Into [4.37] yields: 


d> lr\~ 

14 . 55 " 


[4.38] 


(Within the gradient zone, becomes 4^.) According to eq. [4.38], the 
value of 4 ^ at a certain location, £ , depends only on the local Blot 
number and the local temperature difference between the furnace and the 
charge. 

Constant Interface Shape in the Gradient Zone 

Equation [4.38] suggests. In accordance with Its limitations 
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discussed above, that ^ may be maintained at a constant value over a 
portion of the gradient zone If the product Bl^t^,- Is constant 
over a sufficiently long region near the center of the gradient zone, 
l.e., not near the ends of the gradient zone. The following set of 
"design rules" Illustrates how the heat transfer In the gradient zone can 
be "tailored" so that these conditions may be approximately achieved for 
systems In which the Interface disturbance Is not present: 


(1) Calculate In the gradient zone, assuming that 

Big = 0, by using results of the moving fin model. 

(2) Allow Rq to be sufficiently large so that ^ Is linear. 

(3) Alter the end temperatures of the gradient zone annulus 

so that djjAt, * . In this way, ^ 

constant In the gradient zone, neglecting the effects of 
the crucible near the ends of the gradient zone. In the 
present development. It Is arbitrarily chosen to maintain 

= 1 while altering In this case: 


<t> a hji) * 1 + 

\ <J 
Ac, 

£4.39] 

(<V <P m ) - i 1 

f 1 + \ 4le\ 
l 1 At , ) 

£4.40] 


(4) The desired value of (^determines the required value of 
Bl§* by substituting eq. £4.40] Into eq. £4.38]: 


n .** _ £3.1 4 >& 


£4.41] 


A non-zero value of Bl§* leads to a non-symmetrlc system 
and non-constant (^- thus contradicting the 
assumptions of steps (1) and (3) above. These effects 
are shown, however, to be of secondary Importance In the 
results which fol low. 
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(5) The desired radius of curvature of the Interface, N, Is 
next used to el Imlnate ^ from eq. C4.41]. Using eq. 

CB.6]: 

Bl** = - — i- 4i42: j 

* "» N l + X 6 (<J8„/i^ 

Three test cases are employed to demonstrate the procedure described 
above. The parameters defining the systems and the results of the calcu- 
lations of steps (1) through (5) are shown In Table 4.5. Figures 4.18, 
4.19, and 4.20 plot the resultant temperature and Isotherm curvature 
distributions as calculated by the concentric fin model: part (a) of 

each figure plots the distribution of 0^ln the gradient zone while part 
(b) shows the distribution of ^and tf>„and the Inverse radius of 
curvature, N“1 of the isotherms In the gradient zone. The desired con- 
stant radius of curvature chosen for use In eq. C4.43J is N = 20; there- 
fore, N _1 * 0.05. 

Figure 4.18 considers a system without a crucible and with a 
gradient zone length 1*0. It is seen that the distribution of 
and N tend to flatten out near the center of the gradient zone, approxi- 
mating the desired values of 0 A = 0 00127 and N“^ * 0.05. It Is expected 
that the distributions of 0 & and N will be yet more flat when the 
Influence from the ends of the gradient zone are even more diminished at 
Its center. This Is accompl Ished In Fig. 4.19 by Increasing the length 
of the gradient zone to 2.0. The curves for ^and N are quite flat 
In a large portion of the center of the gradient zone In this case. Note 
that the required value of = 0.000903 for case 2 (c.f.. Table 4.5) Is 
larger than actual ly present as shown In Fig. 4.19. This Is due to heat 
transfer to the charge In the gradient zone, not considered In eq. 
C4.403, which causes to be slightly smaller than the value 

predicted for a symmetric system. Reduced heat transfer to the charge, 
and therefore a reduced value of ^results. If desired, the value of ^ 
can be adjusted upward by an appropriate alteration In Big* 
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hot end C in the gradient cold end 

zone 


Fig. 4.18a. Demonstration of the procedure used to 

obtain an appoximately constant isotherm 
shape in the gradient zone. No crucible, 
symmetric system and Ag = 1.0. 
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Bin* B»q *0.5 

Big = 0.0625 
No crucible 

Symmetric system excepti 

4> g <V 2) * 0 593 

R r = 10 

v 10 

x H * x c * ® 



0.2 


0.1 

N 

0 

- 0.1 


- 0.2 


hot end £ in the gradient cold end 

zone 


Fig. 4.18b. Demonstration of the procedure used 
to obtain an approximately constant 
isotherm shape in the gradient zone. 
No crucible, symmetric system and 
Xq = 1.0. 
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Bl|_j * B'iq* 0.5 
Bi Q * 0.0625 

No crucible 

Symmetric system except! 
<^> G (X G /2)* 0.422 

V° 

V 20 

X H -X C* °° 


» '> ’> 7 \ 
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hot end £ in the gradient cold end 

■zone 


Fig. 4.19b. Analogous to Fig. 4.18b except the 
gradient zone is longer; i.e., 

Xq 3 2.0. 
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1.0 0 1.0 

hot erd C in gradient zone cold end 


Fig. 4.20a. Analogous to Fig. 4.19a except a 

crucible is present; i.e., 6 = 1.5. 



srssa $$ 


8^* Bi c * 0.5 

Bi *0.0735 
8 *1.5 K * 1.0 

Symmetric system except. 
<J> q (X g /2)* 0.481 

v° 

Xq-2.0 

w® 



zone 


Fig. 4.20b. Analogous to Fig. 4.19b except a 

crucible is present; i.e., 6 = 1.5. 
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TABLE 4.5 

Description of the systems used In 
Figures 4.18, 4.19, and 4.20. 



Case 

1 

2 

3 

Figure 

4.18 

4.19 

4.20 

s 

1.0 

1.0 

L . 

Xc, 

1.0 

1.0 

o 

• 

CM 

B! H - BI C 

0.5 

lillliM Hi 

0.5 

Blfi* = Big* 

0.471 

0.471 

0.291 

moving fin model 

-0.407 

-0.289 

-0.259 

(eq. [4.39]) 

0.593 

0.422 

0.481 

(eq. [4.40] 

0.297 

0.211 

0.241 

Big* 

(eq. [4.43]) 

0.0625 

0.0625 

0.0491 

bi g 

0.0625 

0.0625 

0.0735 


0.00127 

0.000903 

0.000811 


Symmetric system except t 0 

K* X t - <*> 

Desired value of N Is 20 
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Figure 4.20, which considers a system similar to that of Fig. 4.19 
except 1.5, demonstrates that the presence of a crucible Is ade- 
quately accounted for in this procedure. 

The procedure outlined In this section Is not Intended as a gener- 
al I zed design method leading to optimized heat transfer for a typical 
Bridgman growth system. Systems of practical Interest do not behave as 
those In rigs. 4.18 to 4.20 due to the generation of radial temperature 
variations at the growth Interface (c.f., section 4.3). l+s Intent, 
rather. Is to accentuate three features concerning the radial temperature 
distribution In the charges 

(1) Heat transfer In the gradient zone Is a useful tool for 
altering the Interface shape In the gradient zone. 

(2) Whore appropriate, results from one-dlmenslonal models 
can be used to help predict radial temperature variations 
In the charge. In the model developed In this section, 
t . [4.38] Is seen to provide at least order of magnitude 
accuracy. 

(3) Generation of radial temperature variations at the 
interface (R* ¥ 1 and/or Pe$Rn ¥ 0) Is detrimental for 
achieving a region about the Interface which has a small 
and relatively constant Isotherm curvature. 


Radial temperature var' ^lons within the cii^-'ge have been shown In 
section 4 to attenuate r«r-. •. with axial distance from the location at 
which they are generated. This behavior permits the formulation of a 
thermal design criteria, leading to satisfactory Interface shape, for 
systems in which radial temperature variations created at the Interface 
are not present. In summary: 
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The Interface should be located In a gradient zone which 
Is made approximately adiabatic through the use of a 
thermal Insulating material (l.e., low kg) placed between 
the hot and cold zones and radiation shields In the fur- 
nace cavity gap to decrease kg. The gradient zone should 
be sufficiently long so that radial temperature gradients 
generated In the hot and cold zones have largely attenu- 
ated at the Interface location; this length Is approxi- 
mately one to two charge diameters and depends on the 
effect of the crucible In reducing the rate of attenuation 
of these radial gradients (see section 4.2.2). The exact 
location of the Interface within the gradient zone Is 
chosen so that the shape of the Interface Is slightly 
curved toward the solid In order to Inhibit the propoga- 
tlon of defects generated at the boundary crucible sur- 
face. Additionally, small amounts of heat transfer to the 
charge In the gradient zone can be used to alter the 
Interface shape as desired, or to maintain a large axial 
region of nearly constant Isotherm curvature (c.f., sec- 
tion 4.4.3). 

The thermal design Is complete at this stage If the axial temperature 
gradient In the melt at the Interface Is sufficiently large to satisfy 
the const I tut Iona I y supercooling requirement. If the axial gradient Is 
not sufficiently 'arge, however, techniques to Increase Its value can be 
Investigated while attempting to maintain satisfactory Interface shape. 

Radial temperature variations generated at the Interface are, how- 
ever, the dominant factor with regard to the Interface shape. The proce- 
dure outlined above, therefore, does not lead to desirable Interface 
shapes unless ,ne Interface effect can be eliminated. This Issue Is of 
partlc lar Importance for the growth of semiconductor materials which 
have values of R* greater than unity; for such materials, the Interface 
effect tends + o produce an Interface shape wh'ch Is curved In the adverse 
direction, l.e., toward the melt. 
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Counteract 1 ng the Effects of Radial Temperature Variations 
Generated at the Growth Interface 

The analysts of section 4.1 suggests that heat transfer from the 
furnace to the charge tends to produce Isotherm shapes which are concave 
toward the solid. Such heat transfer In the gradient zone thus serves to 
make the Interface shape more favorable If It Is originally concave 
toward the melt. The gradient zone. In such a case, resembles an exten- 
sion of the hot zone furnace. The potential for completely eliminating 
an adverse Interface curvature may be estimated by comparing the value of 
<b. produced at the Interface due to the Interface effect (eq. C4.30]) to 

A 

a representative value of ^ caused by heat exchange with the furnace at 
the preval I Ing Blot numbers. The value of ^ at the end of an adiabatic 
gradient zone (see section 4.2) may be used to estimate this latter value 
of ^ . If radial temperature gradients caused by heat transfer to the 
charge are of comparable magnitude, or larger, than those generated due 
to the Interface effect, such heat transfer In the gradient zone may 
prove useful In reversing an adverse Interface curvature. If, on the 
other hand, radial temperature gradients caused by heat transfer to the 
charge are less than those generated at the interface, this technique can 
not be expected to change the direction of the shape of the Interface. 

A decrease In the value of £ likewise decreases the magnitude of the 
radial temperature gradients generated at the Interface, as shown In 
Fig. 4.12. The value of Scan be decreased by either an Increase In the 
charge diameter or by a decrease In the thickness of the crucible. A 
decrease In S, while not eliminating the Interface effect, may be useful 
In combination with heat transfer to the charge as described In the 
previous paragraph. 

The Ideal method of dealing with the Interface effect Is to elimi- 
nate It entirely. If this Is possible, the design procedure for obtain- 
ing the desired Interface shape, described In section 5, once again 
becomes useful. The cause of the Interface effect Is the change In the 
axial temperature gradient within the charge at the growth Interface 
which Is not present In the crucible. Efforts to eliminate the Interface 
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effect therefore must either (1) remove the change In axial gradient 
wtthln the charge at the Interface or (2) produce an equivalent change In 
axial gradient within the cruclblo. 

Equation [4.26] suggests that the slope change at the Interface 
disappears for an appropriate value of the lowering rate, V, when < 1. 
The efficacy of this technique depends, however, on the Interplay between 
the va I ue of V requ I red to e 1 1 m I nate the s I ope change at the I nterf ace 
and the effect of V on the constitutional supercooling requirement. 
Further, this technique Is not applicable when R« > 1. 


The electron potentials are not equal in the solid and liquid por- 
tions of the charge. The passage of an electric current across the 
growth Interface produces or absorbs heat (depending on the direction of 
the current) by virtue of the Peltier effect. Such heat can be used to 
eliminate the change In axial temperature gradient at the Interface 
caused by the change In thermal conductivity and the liberation of latent 
heat. Further, this technique also counteracts the decrease In the axial 
temperature gradient In the I (quid at the Interface which occurs when 
Rk > 1. Rewriting the heat balance at the Interface In order to Include 
Peltier heating/cooling yields: 




S e T, i." D 

-\c) 


[5.1] 


where: G[ * Gj In order to ellmlnai the slope change 

at the Interface 
Se * Seebeck coefficient 
I" * current per unit cross sectional area 
SeT|l" * rate of generation of Peltier heat per unit 
cross sectional area 

T] * absolute solidification temperature 


Using germanium as an example (Se * 70 pV/K, T| ■ 1211 K), eq. [5.1] 
uhows that the required current density Is In excess of 200 amps/cm? for 
typical operating conditions. Such large current densities may be diffi- 
cult to attain experimentally and may produce significant Joule heating 
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effects. Although each case should be examined Individually, It appears 
that Peltier heating Is not a vlab'e method for el Imlnatlng the slope 
change at the Interface* 

Peltier heating/cooling has, however, greater potential for modi- 
fying the Interface shape In small amounts. For example, a flat Inter- 
face shape may be made sllghlty concave toward the solid by the applica- 
tion of a small amount of Peltier cooling. This technique yields a 
similar result as that of transferring small amounts of heat to the 
charge In the gradient zone as described in section 4.4.3. 

The Interface effect may also be eliminated by producing a slope 
change In the crucible at the Interface which Is equal to that In the 
charge. This may be done by adding the appropriate amount of heat to the 
crucible In a very narrow region near the axial location of the Inter- 
face. The remainder of the gradient zone may be adiabatic If desired. A 
thin crucible, being to a greater extent thermally one-dimensional than a 
thicker crucible. Is more desirable In this respect. 

Recommendat Lon For Further Research 

This work considers furnace boundary conditions as an Independent 
variable which. In combination with the other system parameters, deter- 
mines the Interface shape. It would be useful, however, to be able to 
solve the Inverse problem. I.e., what must the furnace boundary condi- 
tions be In order to produce a desired Interface shape? There may be 
multiple solutions or no solution at all. In the first case, such 
results would provide a direct Indication of the thermal design required 
for satisfactory interface shape; In the second case, the knowledge that 
no set of furnace boundary conditions exists which can produce the 
desired Interface shape Ind'cates that other parameters of the growth 
system require alteration. The ability to make such conclusions requires 
i computer solution of the appropriately modeled Inverse heat transfer 
problem. 


« 
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The harmful consequences of the Interface effect for semiconductor 
crystal growth warrant an Independent experimental Investigation. Growth 
experiments utilizing variable thickness and thermal conductivity of the 
crucible and thermal conductivity ratios of the charge (metals for RK < 1 
and semiconductors for R* > 1) will provide results for comparison with 
the analysis of section 4.3. 

The technique presently considered most suitable for counteracting 
the Interface effect Is that of local ly heating the crucible near the 
Interface In combination with a thin crucible. Analysis Is first 
required In order to verify the potential of this method. The relation- 
ship between the thickness of the crucible and the amount and distribu- 
tion of heat transfer to the surface of the crucible to the radla! 
gradients generated In the charge will require the so'utlon ot an appro- 
priate two-dimensional thermal model. Thin crucibles are expected to 
require less heat transfer compared to thicker crucibles In order to 
eliminate the Interface effect. Tnln crucibles of sufficient structural 
strength may be obtained through the use of a metallic crucible; a 
uniform and chemically Inert Interior coating, able to withstand stresses 
created by differential thermal expansion, remains problematic. 
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APPEND IX A; DEVELOPMENT OF THE FIN EQUATIONS FOR THE TWO-DIMENSIONAL 

THERMAL HQPEL 


With the assumptions of Section 3.2, each of the concentrically 
located fins (l.e., Inner charge, outer charge, crucible, gradient zone 
annulus) can be considered as a moving fin. According to Carslaw and 
Jaegar [1], the axial temperature distribution for the fins In the hot 
and cold zones Is given by the following equations: 


A _ A \ i ~ <n~in . _ n 


[A. 1 a] 


- A ~^^f * T 0 - o 


[A. lb] 




kr A, 


fT c 


cr 


cr J i 1 


W W - " 0 


[A.IcD 


where: A * cross sectional area of a fin 

R = thermal resistance per unit length between fins 

subscripts: 

loc * local charge phase, liquid or solid 
In * Inner charge fin 
out * outer charge fin 
cr * crucible 

f * hot or cold zone furnace 
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When considering the fins within the gradient zone, Tf and R cr ,f In 
eq. [A.lc] change to Tq and R C r,G» +he subscript "G" denoting the 
gradient zone annulus. Further, there Is an additional equation represen- 
ting the gradient zone annulus: 



[A. Id] 


The radial geometry of the concentric fins shown In Fig. 3.2 
specifies the cross sectional areas and thermal resistances In eqs. 
CA.1]. The cross section areas Aj n and Ao U f are chosen equal; therefore, 
Dj n * D//2. The radial locations for the representative temperatures of 
the fins are D;/2, D2/2 and D3/2 for the Inner and outer charge fins and 
crucl’jle fin respectively; they are chosen so that there Is an equal area 
within the fin to either side; 


D, = ~-iz 

tA.2a] 

T> 

CA.2b] 


Dj,= 


CA.2c] 


The thermal resistances can then be expressed as (Rohsenow and Choi C4]): 
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[A, 3a] 
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^cr>F Rc^ 
where: S ■ D cr /D 

hioc = heat transfer coefficient from furnace or 

gradient zone annulus to the surface of the 
crucible, based on the outer diameter of the 
crucible 

Substituting eqs. [A. 2] Into eqs. Into eqs. [A.l] and non 
dimensional Izlng the resultant expression yields eqs. [3.2] and [3.3]. 
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This appendix employs simple approximations In order to develop a 
useful relation between Isotherm curvature and the axial and radial 
temperature gradients obtained from one- and two-dimensional thermal 
mode I s . 

An Isotherm within the charge Is assumed to be spherical; Its shape 
can then be represented by a single number — Its radius of curvature. 
The axial distance between the locations of an Isotherm at the center and 
the surface of the charge Is denoted by if. The assumed geometry, shown 
In Fig. B.1, yields the following relation between if and the Isotherm 
curvatures 



KN (NMV /2 ~ | 

£ 4 k; 


far !w|*5) 


[B.1] 


where N Is the radius of curvature of the Isotherm In number of charge 
radii. Equation CB.l] yields a positive value of N when the Isotherm Is 
concave toward the solid as shown in Fig. B.1. If the axial gradient 
does not vary greatly over the cross sections 


r ^ ‘Ac 

* d^/dC 

,her0! P«. a 7>(f- y 0- P(f*o) 

Substituting eq. [B.23 Into [B.1]: 

. , - I d&w/d C 
N 4 ~K~ 


Cb.2] 


[B.3] 


For the purpose of relating <P SC to , the heat conduction equation 
within the charge Is used: 
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CB.4] 


If the Pe term Is neglected and If It Is assumed that Is not a 
function of -f, the resulting temperature distribution In the radial 
direction Is parabolic. In this case. It Is easily shown that: 


Using eq. CB.5J In eq. CB.3]: 

M ~ L CB.6D 

N lb <K 


Equation [B.6] Indicates that Isotherm curvature Is small (I.e., large N) 
when the axial temperature gradient I; large as well as when tiie radial 
temperature gradients are small. 
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Fig. B.l. Assumed geometry used to calculate 
an approximate isotherm curvature. 
N is positive for the isotherm 
curvature shown above. 
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OF InSb IN CONTACT WITH ITS MELT 

The developed theoretical framework presented In the last report has 
been used In conjunction with experimental results (current- I nduced 
growth la/er thickness data) to calculate the Peltier coefficient of 
doped InSb In contact w Ith Its melt. Estimates were made of the error 
resulting from time dependent geometry (A i n(r»/4)2, D = f(t)), and 
composition. The Peltier coeff Iclent thu; obtained was compared with 
those calculated from published Seebeck coefficients (thermoelectric 
power); agreement Is found to be excel leni. 

PROCEDURE AND RESULTS 

The Peltier coefficient Is calculated using equation COs 


J 


co 


and values of Rj and R 2 obtained from the layer thickness and the tem- 
poral characteristics of the current pulses: 



cm 


where: Xj = layer thickness (current-on or current-off) 
tj * duration (current-on or currenf-off) 

The Peltier coefficient. If determined by eq. CO Is directly depen- 
dent on the heat of fusion, H. The search for pr^-;i.ie H data Ir the 
literature revealed values rangl.ig from from ?,'■*( j/g to 213.4 J/g. 
Consequently, the heat of fuslo.i of InSb was experimentally determined 
using a Perkln-Elmer DTA 1700 ‘erentlal thermal analyior *-51h, <1 
Perkln-Elmer System 7/4 thermal analysis ccitroller. Aluminum 
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(H = 398.2 j/g) and antimony (H = 164.8 j/g) were used as calibration 
standards for the system. For a heating rate of 5°C/mln, the data from 
the standards were reproducible and se I f cons I stent to 12 %, The heat of 
fusion for undoped semiconductor grade InSb was measured to be 238.1 
i6 j/g (13.48 kcal/mol). Tne crystal diameter data needed to determine 
the current denity were obtained using a precision micrometer stage 
mounted on a Nomarski Interference contrast microscope (with Integral 
cross-hair... The measurements were made on polished and etched samples 
normal to the growth direction. 

It was found that crucible design and the location of the control 
thermocoup'e have profound Influences on the characteristics of growth 
layers produced oy passing pulses of dlrecT current across the crystal/ 
melt Interface dur- owth. The effect of Joule heating In the hot zone 
on growth layer characteristics varied from extensive to almost non- 
existent (manuscript in preparation). 

The Peltier coefficient of InSb was calculaled from the dimensions 
of contiguous current- Induced growth layers. Representative growth lay' 
hlckr.oss data, illustrating the effect of various crucible designs ativ, 
temperature control schemes, are shown In Figs. 1, 2 and 3. Also shown 
are the calculated Peltier coefficients. Defining the error associated 
with each set of data as tne standard deviation divided by the arithmetic 
average expressed as a percentage. It is seen (Fig. 1) that even though 
J ne measured layer thickness varies by more than 60J&, the calculated 
Peltier coefficient (-0.082 V) has a precision of better than tl2j For 
conditions of current puising with temperature control from the base of 
the crucible (Fig. 2), where the layer thickness varied by more than 11 %, 
the overage calculated Peltier coefficient Is found to be -0.084 V 
The Peltier coefficient calculated from growth layer thick- 
nesses obtained using a special BN/composite crucible and current pulses 
of 19.8 A/cm^ Is -0.097 V (13$) (see . ig. 3). This last result Is from a 
series of curre pulsing experiments with different current densities 
but otherwise identical configurations. Using 9.8 A/cm^ current pulses, 
the average calculated Peltier coefficient Is -0.097 V (tlOf). The 
larger error associated with these data can be traced to an inability to 


180 



ORIGINAL PAG* sy 
OF POOR QUALITY 



CMQU3XN) 333WM3IMJ. U3A.V1 


CA> -dd3Q3 M3I1HW 


181 


Figure 1. Growth with current pulses of short durations graphite 

crucible, control thermocouple at heat pipe. 

a) growth layer thickness 

b) Peltier coefficient calculated from growth 
layer thickness 
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Figure ?. Growth with current pulses of short duration: graphite 

crucible, control thermocouple at base of the crucible. 

a) growth layer thickness 

b) Peltier coefficient calculated from growth 
layer thicknesses 
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Figure 3. Growth with current pulses of short duration: 

Decomposite crucible, control thermocouple at heat pipe. 

a) growth layer thickness 

b) Peltier coefficient calculated from growth 
layer thicknesses 
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duration (20s) this assumption Is Inaccurate and the overall layer thick- 
ness cannot be used as data to calculate k. However, Interface demarca- 
tion superimposed on pulses of long duration permits the measurement of 
the Interface velocity In the critical transition region (current-off to 
current-on) with a femporal resolution that Is consistent with the con- 
straints placed on the use of eq. CO and such experiments are therefore 
expected to yield viable data for the determination of Peltier coeffi- 
cients. The growth velocity changes that occur on the appl Icatlon of 
direct current for extended time periods wero, therefore, also used to 
calculate the Peltier coefficient. The average Peltier coefficient cal- 
culated from three 9.5 A/cm^ long duration current pulses was -0.1064 V 
(+6$). Two long growth layers produced by 17.3 A/cm^ current pulses 
resulted In an average Peltier coefficient of -0.096 V (±2f). Long 
duration growth layers produced by Peltier heating also produced viable 
data for calculation of the Peltier coefficient. The average Peltier 
coefficient calculated from three -12.2 A/cm^ (20s) layers was -0.113 V 
(i3Jli. The results are tabulated In Table 1. 

FRROR ANALY SIS FOR THE EXPERIMENTALLY OETFRMINED PELTIER COFFF I HI PMT 

It should be pointed out that the area of the crystal /me It Interface 
(even for 1 00% central facet) Is not exactly related, through the physi- 
cal constant n/4, to the width of the sample segment. Crystals grown In 
the [ill] crystallographic direction exhibit threefold symmetry and 
related to It a varying degree of external faceting which prevents the 
sysiem from assuming a cylindrical shape. The extent of external facet- 
ing can vary from near zero to 1 00% depending on the radial temperature 
gradient, growth rate acceleration and growth kinetics of the given 
system. For the case of minimal external faceting the relationship A = 
(ji/ 4)D2 (for the area of the crystal /me. it Interface) Is a good approxima- 
tion If the crystals are cut so as to contain the axis of rotation In the 
surface analyzed (as In this study). This relationship Increasingly 
underestimates the area of the crystal cross section as the extent of 
external faceting approaches ICOf. (For this exteme case, the external 
shape of the crystal would be triangular.) In effect, at lOOJt external 
faceting conditions, n assumes the value of 9/4(3) ^ * 3,8971 which Is 
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accurately measure the thickness of the growth layers produced by current 
pulses of such relatively low density. As a result of the small current- 
induced segregation change between growth layers, the differential chemi- 
cal etching, used to reveal the growth layers, was Insufficient to 
clearly delineate the current-on to current-off transition. For the 
reason cited above: greater precision In layer thickness measurements, 
growth layers produced from current pulses of 29.6 A/cm^ led to a calcu- 
lated Peltier coefficient of -0.099 V with a precision of t3jl. In an 
experiment where current pulses of 39.7 A/cm^ were used, the Peltier 
coefficient was found to be -0.087 V with a precision of t6%. An 
explanation for this relatively low value of the Peltier coefficient and 
the low precision (when compared with the other experiments In the 
series) Is given below. 

The Peltier coefficient Is a materials property and as such expected 
to remain constant over the entire growth process and for various growth 
configurations. Conspicuous for these results Is that the same Peltier 
coefficient Is obtained although pulse duration, current density, dopant 
level and dopant type were varied. Moreover, the various experiments 
used as a basis for the calculations were conducted In different thermal 
configurations (l.e. with different crucible and temperature control 
configurations). This finding ts taken as a confirmation of the validity 
and applicability of the theoretical framework. Furthermore, It Indi- 
cates that, for the growth of (short time duration) current Induced 
layers, the effect of bulk heat flow perturbations (Joule heating) can be 
separated from transient Interface heat balance considerations (Peltier 
effect). Thus, fhe long term variation In layer thickness Is controlled 
by the changing temperature gradients In the liquid and solid caused 
principally by transient bulk heat flows. On the other hand, the funda- 
mental relationship between the dimensions of contiguous layers (layers 
generated by a pulse sequence) and the Peltier effect (l.e. eq. CO) 
remains unchanged. 

The basic premise for using eq. CO to calculate the Peltier coeffi- 
cient Is that the Interfaclal temperature gradients remain virtually 
unchanged over a given pulse sequence. For d! r ect current pulses of long 
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TABLE 1 


EXPERIMENTALLY DETERMINED VALUES OF THE PELTIER 
COEFFICIENT OF InSb 


on/ + off 

(sec) 

Dopant 

Current Density 
(A/cm^) 

Peltier Coeff. 
(V) 

Error 

<±*> 

2/6 

Te 

16.8 

-0.082 

12 

2/6 

Te 

13.0 

-0.084 

5 

1/3 

Ge 

9.8 

-0.097 

10 

1/3 

Ge 

19.8 

-0.097 

3 

1/3 

Ge 

29.6 

-0.099 

3 

1/3 

Ge 

39.7 

-0.087 

6 

20/100 

Ge 

9.5 

-0.106 

5 

20/100 

Ge 

17.3 

-0.096 

2 

20/100 

Ge 

-12.2 

-0.113 

3 
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an error In excess of 24?. InSb crystals have been grown exhibiting the 
full range of external faceting though those that reached 100)1 were 
highly susceptible to twinning. Those crystals from which current- 
induced growth data were taken were Intermediate In external faceting. 
It Is estimated that the areas obtained using the relationship for circu- 
lar cross section are In error by less than 10?, An underestimation In 
area translates directly to an error in calculated current density and 
thus leads to a calcualted Peltier coefficient which Is lower than the 
actual value. 


Due to changes In meniscus height and/or perturbations In the bulk 
heat flow, the diameter of a growing crystal Is time dependent In nature. 
In addition, the external morphology will change as the diameter assumes 
a new dimension. The concomittant changes In current density also appear 
as variations In the calculated Peltier coefficient. This Is evident In 
the Increased error associated with the data obtained from layers pro- 
duced by 39.7 A/cm^ pulses (Table I). In the related experiment the 
diameter Is 1.12 cm at the beginning of the pulsing sequence. It Is 
1.19 cm at the end, an area Increase of 11?. The calculated Peltier 
coefficient, when corrected for this change In diameter. Is -0.097 V 
which Is closer to the values (-0.097 V, -0.097 V, -0.099 V) calculated 
for the other current densities (9.8 A/cm^, 19.8 A/cm^, 29.6 A/cm^ 
respectively) In that sequence of experiments. 


The effect of meniscus height change on crystal diameter was Inves- 
tigated analytically using the expression of Bardsley et al (1 *: 


h - 




ill 


_ ft 

4-r 


[3] 


where: h <* meniscus height 

r ■ c-ystal radius 

QT “ surrac* enslon 

*,5 

g « gravitational acceleration 
0 L = angle made to the vertical by the meniscus at 
the llquld/solld Interface 
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The derivative of eq. C3] with respect to r was evaluated to determine 
the change of meniscus height with respect to crystal radius: 


jk. r ; 


$Cei8t 

4r z 


C4] 


For the representative values: 

£ * 0.21 (<r = 600 erg/cm 2 ) 

% - 20 ° 
r = 0.50 cm 

the derivative gives dh/dr = -0.17 cm/cm. During the growth of 20 A/cm 2 
layers (Fig. 3), the crystal radius was found to Increase by 117 pm. 
Thus, according to the above expression (eq. CO, the meniscus should 
decrease by 20 pm. The net Interface displacement calculated from the 
Inll lal layer thickness transient Is 24 urn. Considering the Inherent 
error ’ n calculating the steady state pul I rate (Including melt height 
drop) and the uncertainty In the values of p and 6 L the agreement between 
the calculated and experimental ly determined values of the Interface 
displacement Is remarkable but may be fortuitous. 

C OMPARISON OF THE EXPERIMENTALLY DETERMINED PELTIER CQEE-F-LC I ENT_W 1 TH 
THOSE CALCULATED FROM THE PUBLISHED SEEBECK COEFFICIENTS: 

DEPENDENCE ON DOPANT TYPE AND DOPANT LEVEL 

The Peltier coefficient can be directiy calculated from the measured 
thermoelectric power (Seebeck effect) using the Kelvin relations^ 2 *. The 
Peltier coefficient Is related to the Seebeck coefficient by the absolute 
temperature: 

7T(j'l = «(t)T C5] 

where a If fhe Seebeck coefficient. The Seebeck coefficient, as measured 
by various Investigators, Is tabulated In Table 2. Though not all 
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measurements were made at the melting point It has been shown both theo- 
retically and experimentally that the thermoelectric power of InSb Is a 
very weak function of temperature and Is Independent of doping level and 
type at high temperatures*^'**. The value of r calculated Is the differ- 
ence In the Peltier coefficients cf the liquid and solid at the melting 
point: 


For all calculations the Seebeck coefficient for liquid InSb was taken to 
be -40 nV/°K*^* (except Blum and Ryabtsova*®* who give the difference, 
o g - oj). The experimentally determined values of *\/ s (Table 1) are In 
excellent agreement with those listed In Table 2. 


It has been shown that the Peltier coefficient of semiconductor/melt 
systems can be calculated directly from ih. dimensions of contiguous 
pairs of current-induced growth layers. This technique Is applicable to 
systems which (1) exhibit high resolution differential chemical etching 
behavior or (2) permit measurement of the thickness of current- Induced 
layers (eg. high spatial resolution measurement of changes In dopant 
distribution caused by the current pulses). It should be noted that the 
InSb/melt system Is optimized with respect to the boundary conditions and 
assumptions of this analysis, namely, the heat absorbed at the crystal/ 
melt Interface due to the Peltier effect Is of the same order as the heat 
released as a result of solidification (l.e. heat of fusion). In this 
way, the effect of the current pulses was to produce a substantial change 
In the microscopic rate of growth for that portion of the duty cycle when 
the current was on. This permitted the accurate determination of the 
Peltier coefficient since It was directly >ro?ortlonal to the difference 
In the rate of growth current-on/current-of f. For other systems. It will 
be necessary to choose an Impressed pulling rate and current density that 
are commensurate with this condition. The assumption of constant rates 
of growth during both portions of the pulsing sequence should, be verified 
through superimposed Interface demarcation. 
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TABLE 2 

PELTIER COEFFICIENT OF InSb CALCULATED FROM PUBLISHED 
VALUES OF THE SEEBECK COEFFICIENT 


Investigator 

Seebeck Coeff. 
(|iV/°K> 

BaJUfir Coaff. 

(V) 

Rose et al 

-130 (773°K) 

-0.0715 

Glazov <s)< 5 > 

-150 (798°K) 

-0.0878 

Glazov (l) (5) 

-40 (803°K) 


Tauc< 3) 

-150 (750°K) 

-0.0878 

Ehrenrelch (calc.)* 4 ^ 

-160 (750°K) 

-0.0958 

Blum and Ryabtsova (s/i)^ 

-235 (798°K) 

-0.1875 
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III. COMPUTER OPERATED BRIDGMAN GROWTH SYSTEM 

WITH GRADIENT CONTROL ZONE 
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A Br Idgman-Stockbarger crystal growth system which allows for 
controlled heat transfer during crystal growth has been designed and 
constructed. It consists of two vertically aligned Na-fllled heat pipes 
separated by a gradient control region which Is variable In length as 
well as composition. The furnace operating temperature range Is 500°C to 
600°C In vacuum and Inert gas atmospheres, excluding hydrogen. The 
furnace can accommodate crucible sizes of up to 1.9 cm diameter and 46 cm 
In length. 

Two dc servo motors control led by a standard dost .-loop servo 
system provide for crucible lowering/raising and rotation functions. An 
optical encoder on i he motor shaft al lows the crucible position to be 
monitored at all times. Crucible lowering rates can be varied from 
0.036 cm/hr to 180 cm/hr using a gear reduction from 171 to 171,000; 
crucible rotation can be varied from 0.05 rpm to 30 rpm. All motor 
functions controlling lowering/raising are digitally Interfaced. 


DIGITAL TEMPERATURE CONTROLLER 

The temperature control function of the two Independently operated 
heal pipes Is performed digitally. Although there Is some degradation In 
the monitoring signal (tempera + ure In mill Ivolts) because of Inherent 
samp 1 1 ng effects, there are significant advantages to digital control 
which Include the greater range of compensators and time varying control 
functions. 

The first step In the design of the compensator Ir.olved the formu- 
lation of a model of the open- loop system performance, l.e., relating the 
measurable outputs to the variable li.puts. T e temperature of the heat 
pipe is related to the power Input by a first-order differential equa- 
tion: 


A 
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Ht = C + RC 

where: T * temperature 

T a * ambient •. imperatcre 

Q = powr (in Watts) 

R » overall furnace resistance (°C/Watt) 
c = overall furnace capacitance (Watt*sec/°C) 
t = time 

The control ler could then be designed according to the following 
control specifications: 

(1) No steady state error to step Inputs of power; 

(2) No overshoot (undershoot); 

(3) Temperature settles In <500 secs. 

A P~l (proportional and Integral) compensator was designed In the 
continuous time domain: 

K p (s + K i ) 

PRC - t. 


(K p R + l) 2 

< 

4R 2 K 0 C 

t s Is + ho desired settling rime. 

In this manner the control variables, K p und Kj. could be calculated 
based on the furnace parameters and control requirements. 

Character l zaf Ion of the furr.ace parameters P nd C was done by 
measuring the static and dynamic open-loop sys+em resjjonse of temperature 
for step Inputs of power. 


K(s) 


and 


wh^ro 
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Experiments revealed that (a) the heating efficiency Is greatest for 
vacuum and deteriorates as the Inert atmosphere changes from argon to 
helium because heat losses Increase by conduction and (b) the system time 
constant decreases as the atmosphere changes from vacuum to argon to 
het turn. 

Sire > he compensator Is implemented digitally, it must be discrete 
a e well as the temperature Input to the controller. A 16 bit A/D 
converter with a theoretical temperature resolution of t0.04°C for a 
chrome I /a i ume I thermocouple is used to sample the furnace temperature. 
Discrete representation of the control ler Is achieved by using the bi- 
linear transformation: 

M(nT) ■= M[(n-1)T] + E(nT)[K p *^] t E [(n-l)T][-K p + i2|iI] 

where: M = control response, i.e., power Input 

E = error = temperature set point temperature 
T = sampl Ing time 
n = present time step 
(n-1) = one time step prevtous 

The output from the computer Is discrete and is reconstructed to give a 
continuous output to the furnace using a D/A converter. The functional 
block diagram of the closed- loop system is: 



During the control loop, five thermocouples can be monitored, the point 
reference Is timed digitally and power Is measured at each time step. 
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CHARACTER I ZAT I ON OF THFRMAL BEHAVIOR 


Furnace characterization for determination of the temperature 
distribution In the heat pipe regions and gradient relgon has been 
Initiated. The temperature was measured along the central axis of a 
graphite rod completely filling the heat pipe cavity. Measurements were 
made for vacuum (400-600 >im), argon and he 1 1 urn w Ith T H = Tq * 600°C and a 
gradient region length of 2D (D = diameter of graphite * 0 . 75 "). The 
material used for the gradient region was an Insulating, low density 
AI 2 O 3 cylinder. The measurements show that within the top heat pipe, as 
anticipated, the graphite rod Is about 65 % Isothermal. Near the top, the 
temperature decreases because of heat losses to the surroundings; the 
femperature decrease Is greatest for helium and smallest for vacuum. In 
the lower heat pipe there Is a significant drop In temperature toward Its 
lower end. For helium only about 45J of the cavity length Is Isothermal; 
for argon and vacuum, the temperature decrease extends from top to 
bottom. The non- Isothermal behavior orglnates In (a) the heat losses at 
the enci of the graphite and In (b) the reduced thermal coupling of the 
heat pipe to the graphite In lower conducting atmospheres. In the gra- 
dient region the temperature drops off about 6 °, Indicating radial heat 
losses to a non- Idea I Insulator. 

Characterization of the thermal behavior of the furnace cavity Is 
being continued Iteratively making use of the two-dimensional heat trans- 
fer model so as to arrive at an optimized configuration for the gradient 
control region. 
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IV. MAGNETIC FIELD ASSISTED MELT STABILIZATION 
IN VERTICAL BRIDGMAN GROWTH 


SUMMARY 

Making use of the prototype ASTP multipurpose furnace, the effects 
of transverse magnetic fields of up to 3.8 kGauss on macro- and micro- 
segregation was Investigated during growth of Ga-doped germanium. This 
particular system was selected for the study since It permitted a compar- 
ative analysis with results obtained during growth In a re' <~ed gravity 
environment. 

It was found that the macrosegregation behavior for the same c ling 
rate as used In the space experiment was significantly different ar was. 
In fact, similar to that observed for regular ground based growth experi- 
ments. The results Indicate clearly that while transverse magnetic 
fields will suppress turbulent convection, they will not result In diffu- 
sion controlled segregation during vertical Bridgman growth. 

The microsegregation behavior In the presence of a magnetic field is 
of Interest Insofar as current Induced Interface demarcation generates 
distinct magnetic field Induced segregation anomalies which are normally 
not encountered. While no comprehensive analysis of this effect has thus 
far been conducted, there are strong Indications that the electric field 
Induced perturbations of segregation can provide details concerning the 
solute boundary layer characteristics for magnetic field assisted growth. 

Evidence was found for a reduction In the degree of radial segrega- 
tion asymmetry. This finding Is of Interest since the change In radial 
segregation was not associated with a measurable change In growth Inter- 
face morphology. It suggests that the convective melt flow pattern, 
rather than the growth Interface morphology, controls radial segregation, 
and that the melt flow pattern changed with the application of the 
external magnetic field. 
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Fluid Flow in Crystal Growth: 
Analysis of tha Floating Zona Procass 

Principal Investigator: Robert A. Brown 


SUMMARY 


This research program is aimed at a fundamental understanding of 
the interactions of heat, mass and momentum transport in crystal growth 
from the melt. Emphasis has been on studies of the small-scale floating 
zone process and on a prototype of the vertical Bridgman growth system. 

In both systems detailed numerical calculations have been used to disect 
the interplay between fluid convection and dopant segregation. Calcula- 
tions are based on newly developed finite-element techniques that make 
feasible the complete solution of solidification problems that include 
convection in the melt, and on computer-implemented perturbation analysis 
of solution multiplicity and stability. These numerical methods have also 
been successfully extended to include the interactions between a melt/gas 
meniscus with heat transfer in meniscus-defined crystal growth systems and 
to model. the transition to a cellular morphology in a solidification inter- 
face. 


1 . Scope of Research Program 

Our research is composed of a set of theoretical and computational 
studies directed at a fundamental understanding of the interactions of 
fluid mechanics and heat and mass transport in crystal growth from the 
melt, especially by the vertical Bridgman and floating zone process. 
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The techniques for analysis and the quantitative insights developed in 
this research have broad application in a variety of melt growth systems 
(see list of Publications at the end of this report) and have made it 
possible to model Bridgman-Stockbarger and floating zone systems that 
are being used in ground-based experimental studies. Combining our 
detailed models with careful experiments hopefully will lead to quanti- 
tative pictures of the transport phenomena in each growth system and to 
well controlled space experiments. 



Figure 1 : 

Schematic of Floating Zone Process. 


Modelling of a melt crystal growth system is a formidable task. Most 
of the complications are seen by considering the floating zone process pre- 
sented schematically in Fig. 1. Here a cylindrical rod of a multi -component 
and polycyrstalline material is melted and re-solidified into a single 
crystal by using a short circumferential heater (either resistance or RF) 
that translates slowly along the axis of the rod. The molten zone forms 
between the crystal and feed rods and is held in place by surface tension. 
Its size is governed by heat transfer and limited by instabilities that 
originate at the melt/gas meniscus. In practice on earth, the radius of 
the rods and the height of the zone are nearly equal and range roughly be- 
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tween 0.5 and 3 centimeters. A quantitative understanding of the physics 
of a small-scale floating zone process requires analysis of the couplings 
between the convection in the melt, the heat transfer in melt, crystal 
and surroundings, the shapes of the melt/solid and melt/gas interfaces 
and the transport of solutes In the melt. 

Our first approach toward understanding these interactions has been 
to decouple them along lines dictated by the physics of typical semicon- 
ductor growth systems; the separations used here are shown on Figure 2 
where the numbers in parenthesis refer to publications listed on page 
3. In the first three years of NASA support, research projects have fo- 
cused on analysis of heat transfer in growth systems with melt/solid and 
melt/gas interfaces, on the convection in melt growth driven by crystal 
rotation, buoyancy and surface tension differences and on solute segrega- 
tion caused by these flows. 

In each research project, a combination of analytical modelling and 
computer-aided calculations is being used to develop both a qualitative 
understanding of the relevant physics and quantitative models of the ver- 
tical Bridgman and floating zone processes. To do this, state-of-the-art 
numerical methods have been developed for handling solidification problems 
and for studying the multiplicity of flow fields in nonlinear convection 
problems. These methods give us the capability of accurate calculations 
of the flow fielu and mass transfer near the solidification interface. 

This ability has been demonstrated through analyses of the effects of melt/ 
solid interface shape and natural convection on radial segregation in a 
prototype of the vertical Bridgman system [7-fi] and by a study of the evol- 
ution of a microscopic cellular morphology at a solidification interface 
[12-13]. Both studies are highlighted below. 

The Bridgman calculations will be compared directly to experimental 
results gathered in the new vertical Bridgman system under construction in 
A.F. Witt's research group at M.I.T. This collaboration promises the first 
detailed comparison between the modelling of convection in melt growth and 
experimental segregation results. Similar calculations are underway for 
small-scale floating zone system and we hope to compare these to experiments 
that are being carried out at Westec Corporation in Pheonix. 
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To understand the complicated interaction between solute segregation 
and natural convection in melt growth geometries we have analysed the de- 
tails of the flow structure and transitions in natural convection in rr.*--' *• 
confined within cylindrical and rectangular containers. These stud’ev 
have lead to systematic methods for tracking changes in steady state flows 
and for detecting the onset of time-periodic convection that leads to axial 
dopant stria'. ions in melt growth (Hurle et al . 1974 and Carruthers 1976), 
These calculations and the insight they supply for simulation of melt crys- 
tal growth are described in [5, 7, 9]. 


The four research projects described here and the personnel involved 

are: 

1. Coupling between rotationally-dri ven , surface-tension-driven, and 
electromagnetically-driven convection and solute transfer in a 
small-scale, gravity-free zone (6.M. Harriott and J. Duranceau). 

2. Fundamental studies of flow transitions in natural convection (Y. 

"k ★ 

Yamaguchi , C.J. Chang , P. Sackinger). 

3. Effect of buoyancy-driven convection on melt/solid interface shape 
and radial segregation in vertical Bridgman growth (C.J. Chang*). 

4. Microscopic modelling of melt/solid interface dynamics (L.H. Ungar). 


RESEARCH SUMMARY 


Fluid Mechanics of a Small-Scale Floating Zone 

Analytical and numerical methods have been applied to study the flows 
driven by rotation and surface-tension-differences is a prototype for the 
small-scale floating zone system proposed for space experiments. Both mech- 
★ 

Graduated during period '0/1/81 - 9/31/82. 
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anisms for convection are important for these experiments. Marangoni flows 
driven by temperature (and possibly composition) differences may be un- 
avoidable in space, while convection induced by deliberate crystal rotation 
causes well characterized radial segregation patterns that can be used to 
calibrate experimental results to theory. 

G.M. Harriott has developed analytical methods for calculating the 
cellular flows driven by these mechanisms and the deflections of the melt/ 
gas meniscus caused by these flows when the zone is almost a circular cy- 
linder. The studies of rotating flows are based on expansions for small 
Reynolds number Re h nRH/v , where R is the radius of the rods, H is the 
height of the zone, and v is the kinematic viscosity. The details of this 
aralysis are contained in the manuscript [6]. One major result is the cell 
ular structure of the flows produced by rotating the too rod at the rate 
sn and the bottom surface at rate n . Depending on the rotation ratio s 
and the aspect ratio of the zone ^ * R/H , there may by either one or two 
toroidal secondary cells in the zone. For a range of these parameters the 
second cell does not fill the cross-section of the zone and a stagnation 
line may intersect the surface of t.he growi g crystal. This stagnate’' 
surface can have a drastic effect on solute segregation as discus 
[6]. The flow structures calculated for small Reynolds number a.c . jwn 
on Figure 3. 

Harriott's analytical results have been compared to detailed finite 
element calculations that are possible for a wide range of Reynolds number 
[10]. The asymptotic results are within ten percent of the numerical solu- 
tions up to Re = 50 and hence are very useful estimates of the fluid flow 
in the microzone. These flows have also been used in radial segregation 
studies [16]. 

Finite element calculations for flows with exact counter-rotation (s = 
-1) at higher Reynolds numbers have uncovered multiple steady flows beyond 
a critical value that depends on aspect ratio, as depicted in Figure 4. 

The new flow fields do not have the reflective symmetry between the two 
toroidal cells normally expected for s = -1 and e/olve toward increasing 
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Reynolds number. Some flows in the symmetric and asymmetric flow families 
at the same value of Re are shown as Figure 5. Mathematical arguments in- 
dicate that the stability of the fluid flow changes at the bifurcation 
point so that only the asymmetric flows would be observed experimentally. 
Also, as explained in [10], changing from exact counter-rotation or add- 
ing gravity to deform the cylindrical meniscus ruptures the bifurcation 
point into two separate families of flows with stable flows existing on 
each family. These are new and exciting fluid mechanical predictions and 
experimental verification is underway using the micro-zone apparatus con- 
structed by A.D. Little for NASA. 

The mathematical concepts and numerical techniques for detecting, 
tracking and calculating the stability of multiple steady states was first 
laid out by L.H. Ungar n our research group in a related study [3] of 
the effects of melt volume, gravitational acceleration, and rotation rates 
on the stability of a "igidly rotating liquid zone held together by surface 
tension. Harriott has also applied these ideas to analyze the effect of 
the rotational ly-driven flows on the maximum length of a cylindrical zone. 

Our approach to the study of Marangoni convection in the floating 
zone geometry has raralled the research on rotating flows described above. 
First, analytical forms, similar to those derived in [6], have been calcu- 
lated for flow fields that are valid in zones only slightly distorted from 
a cylincor at Ijw Marangoni and Prandtl numbers. The results are cast in 
terms of a temperature distribution (see Figure 1 for a qualitative pic- 
ture) that is incident on the meniscus and the structure of the cellular 
flows Inat result from the differences in surface tension. The constraints 
of the melting temperature at the crystal faces limit realistic temperature 
profiles to cubic polynomials with the maximum temperature between the two 
ends. This distribution always results in two toroidal cells that flow 
inward at the crystal faces and so are opposite from the flow induced by 
counter-rotating the seed and crystal. The combined effects of Maragoni 
and rotationally-induced flows are presently being studied to determine 
the strategy for crystal rotation that leads to the smallest amount of rad- 
ial segregation of dopant. 
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Real floating zones have neither flat melt/solid interfaces nor 
cylindrical menisci. J.L. Duranceau has adapted finite -element solu- 
tion schemes to solve for the rotating flows In deformed floating 
zones. Meniscus shapes are calculated by solution of the Young-Laplace 
equation (see [3]) which neglects the coupling of viscous normal stress 
to the deformation of the zone. The analysis in [6] and computational 
checks in [14] show this to be an extremely good assumption for the low 
capillary numbers (C^ = yftR/a) characteristic of silicon floating zones. 
Figure 6 demonstrates a typical change in the cellular flow structure 
caused by the meniscus deformation for a one centimeter diameter silicon 
zone. Changing the shape of the melt/ gas interface has caused an order- 
of -magnitude decrease in the intensity of the top cell; this change will 

have a significant effect of radial dopant segregation. 

To relate these results to actual floating zones requires heat 
transfer analysis for determining the proper length and volume of the 
melt as well as the shapes of the freezing and melting phase boundaries. 
The framework for this analysis has been laid in related studies by Et- 
touney ([1-2], also see Ettouney et al . (1982)) of Edge-Defined Film-Fed 
Growth of thin silicon sheets and the analysis of the floating zone is 
underway. 


Fundamental Studies of Flow Transitions in Natural Convection : 

The complicated patterns for buoyancy-driven convection in melt 
crystal growth systems anc! the importance of the location in terms of 
temperature gradients of the transition between laminar and time-periodic 
convection have emphasized the need for the development of numerical 
techniques for probing structure and transitions in natural convection. 

We have undertaken such an investigation while focusing on the simply 
defined heat transfer systems of a vertical cylinder heated from below 
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and a rectangular cavity heated from the side; these two systems are des- 
cribed next. Finally, the analysis has been extended to Include the In- 
teraction between natural convection and a deformable melt/sol Id phase 
boundary In a study of flow structure In a two-phase analog of the verti- 
cal cylinder problem. 


The Vertical Cylinder Heated From Below 

The vertical cylinder configuration was chosen because of geometrical 
similarities to the vertical Brlagman system and because of the availa- 
bility of previous theoretical (Charlson and Sani 1970, 1971, 1975) and 
experimental (Liang et al. 1969, Olson and Rosenberger 1980) results for 
checks of our results. We have combined finite-element techniques with 
the computer-aided methods developed In [3] for detecting and tracking mul- 
tiple solutions to map out the steady axl symmetric flows that satisfy the 
Boussinesq equations for natural convection. The results of this study 
as discussed in [5] and in more detail in [9]. A typical solution structure 
is shown as Figure 7 for a cylinder with height L twice its radius R and a 
r nducting sidewall; here the Nusselt number for the top surface Nu* of the 
cylinder Is plotted against Rayleigh number Ra ■ eg(T H - T c )l /a) for each 

family of flows identified in our calculations. The static state corresponds 

t 4 

to Nu «1 and is the stable state up to the first critical value Ra * 1.2 x 10 

where two families of single-celled flows are stable and equally probable 
In experiments. The flows In these families differ only in the direction 
nf the circulation, either up or down along the center of the ampoule, and 
are otherwise identfcal . Our predicted location for the critical point Is 
within three percent of the value calculated by Charlson and Sani 0970) 
using linear stability analysis. 

For a Prandtl number of one, the single-celled flows increase in in- 
tensity with increasing Ra until both families turn toward decreasing values 

4 

of Rayleigh number at a limit point near Ra * Ra^ = 6 x 10 . This limit 

point marks the loss-of- existence of single-celled flows for increasing Ra 
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rigure 7. Bifurcation diagram for a vertical 
cylinder with conducting sidewall; 
Pr = 1 and A = 0.5. 
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Figure 8. Loci of three lowest critical Rayleigh nwnbers as a 

function of aspect ratio for a cylinder heated from below 
with rigid sidewalls. 
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and may well mark the upper bound for steady state axisymnetrlc flows. 

As depicted In Figure 7, the two families of single-celled flows join 

families of flows with two cells stacked axially In the cylinder. These 

two cell flows branch from the static state at a second critical Rayleigh 

a 

number Ra » 3.1 x 10 

Extensive calculations reported In [9] have shown that the limit 
point In the primary family of cellular flows (see Figure 5) is generic 
for Prandtl numbers between 100 and 1, for either insulated or conducting 
rigid sidewalls, and a wide range of aspect ratio. This finding differs 
from results obtained in our research group and by Jones et al. (1976) 
for a cylinder with shear-free walls, where the first family of cellular 
flows that bifurcate from the static state seem to evolve continuously 
toward higher values of Ra. We believe that the differences in the flow 
structures for cylinders with slippery and rigid walls is linked to the 
breaking of second-order critical points that exist for the problem with 
shear-free sidewalls. The loci of the lowest three critical Rayleigh 
number for a cylinder with rigid, insulated sidewall is shown in Figur® 8-, 
each curve denotes the location of simple bifurcation points and inter- 
sections correspond to second-order bifurcation or double-points. Note 
that the curve for the mode occurring at the lowest value of Rayleigh 
number does not Intersect the curves for higher modes and that the form 
of the cellular flow that evolves from each critical value changes con- 
tinuously along the lowest mode. This behavior is distinctly different 
from the result for the cylinder with shear-free sidewalls (Liang et al . 
1969 and Rosenblat 1982) where the form of the convection changes abruptly 
at double points. We feel that the behavior of the critical curves and 
the existence of the limit points for rigid walls can be explained through 
the transition between slippery and no-slip boundaries acting as a higher- 
order imperfection that destroys double points. Asymptotic analysis is 
underway that will hopefully confirm this hypothesis. 
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The Onset of Time-Periodic Flows In a Rectangular Cavity Heat From 
the SW . 

Although the Importance of the transition between laminar and time- 
periodic convection has been well established (see the review paper by 
Carruthers (1976)) few theoretical or experimental studies have been 
directed at quantitatively understanding this transition In systems of 
Interest In crystal growth of semiconductors. The most notable experi- 
ments are those of Hurle et al. (1974) on the onset of oscillatory con- 
vection In a rectangular boat of gallium heated from the side. Crochet- 
et al. (1982) have used time-dependent finite difference methods to track 
the evolution of a GaAs melt In a rectangular boat with slightly different 
thermal boundary conditions from those considered by Hurle and have found 
values of Rayleigh number where the flow Is oscillatory. 

Our approach to detecting the onset of these oscillations has focused 
on the approach of bifurcation theory (Iooss and Joseph 1980) that Iden- 
tifies the onset of stable time-periodic solutions as a Hopf bifurcation 
point between the family of time-periodic states and one of laminar flows. 

As outline In the Ph.D. thesis of C.J. Chang and In the manuscript [11], 
these bifurcation points occur where a pair of complex conjugate eigen- 
values of the linearized and discrete version of the Bousslnesq equations 
cross the Imaginary axis. Chang has calculated these eigenvalues for the 
horizontal boat geometries considered by Hurle and Crochet. Sample Iso- 
therms and streamlines are shown In Figure 9 for Pr ■ 0.015 and an aspect 

ratio (length to height) of two; the stability of each flow Is denoted 

2 

along with the eigenvalues (scaled with v/H ) for the unstable modes. Com- 
plex eigen-values correspond to the onset of oscillatory convection. 

One of the most Important results obtained thus far Is the sensitivity 
of the structure of the laminar flow and the critical Rayleigh number for 
the onset of the convection to the each value of the Prandtl number. This 
point Is made by the results for Pr « 0.03 shown on Fig. 10. Here the first 
transition detected from the base laminar flow Is to a second family of 
steady flows, not to an oscillatory state. The structure of the new laminar 
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Figure 9. Linear stability analysis of Hurle's experiment for Pr = 0.15, 
L/H = £ , insulated top and bottom and rigid walls. 
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flow for Pr » 0.03 and further transitions tc oscillatory states will 
be discussed in [11]. 


Effect of Melt/Solid Interface on Flow Structure and Transitions 

An Important consideration 1r. the analysis of convection and mass 
transpor* In melt crystal growth is the role of the melt/solid in in- 
fluencing the structure of laminar flows and the location of transi- 
tions to time-periodic convection. One goal of our research has been 
to extend the methods for computer-aided analysis of nonlinear structure 
to solidification systems that combine both buoyancy-driven convection 
In the melt and the complications of the unknown location of the melt/ 
solid phase boundary. We have successfully reached this goal through 
a sequence of calculations using progressively more advanced numerical 
methods . 

Finite element methods [2] were developed first for steady solidi- 
fication models that neglect the fluid flow in the melt. The reasoning 
behind the use of these simple models Is two-fold. First, they are good 
descriptions of heat transport for many crystal growth geometries, es- 
pecially in a low-gravity environment; our notion of such an "operating 
model" for a meniscus-defined growth technique is displayed in the work 
of Ettouney et al. (1982). Secondly, novel numerical methods are more 
effectively developed for these less complicated systems and then extended 
to account for bulk convection in the melt. This strategy has been success- 
ful in the development of finite-element techniques by Chang [4] for analy- 
zing the interaction of melt/solid Interface shape and buoyancy-driven con- 
vection In a model system where melt and solid are held in a cylindrical 
ampoule placed in a vertically destabilized temperature gradient. 

When the cylinder Is heated from below (Figure 11a) the meltAolid 
system Is a two-phase analog of the natural convection problem considered 
above and in [5,9], The same techniques developed for analyzing multiple 
steady flows in the single-phase problem have been extended to the 
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(b) PROTOTYPE OF VERTICAL BRIDGMAN- 
STOCK BARGER SYSTEM 


Figure 11. Model two-phase systems: (a) Raylelgh-Benard and 
vertical Bridgman geometries. 


m/mm 




two-phase geometry to understand the role of the melt/solid Interface 
on the structure of laminar flows. As is seen by comparing the bi- 
furcation diagrams Figures 7 and 12, the oresence of the melt/ sol id 
Interface has very little effect on the structure of the weakly non- 
linear cellular flows that originate along the static state, but 
changes qualitatively the structure of the flow families far from the 
bifurcation points. Most noteably, the 1U and ID flows result In 
different shapes of the melt/solid Interfaces as shown by Figure 13. 

The loss of synmetry between the 1U and U flow fields signals 
an imperfection that ruptures the secondary bifurcation po'nt that ex- 
ists along the 2U family In the single-phase problem and results in 
the continuous transition In the two-phase system between flows in the 
ID and 2U families. 

Although the boundary conditions for the two-phase Rayleigh-Benard 
problem are extremely idealized, the flow structure for this problem 
gives much detail about the more complicated vertically destabilized 
Bridgman system (see Fig. lib) of Interest In crystal growth. As is 
discussed in [7], the Rayleigh-Benard and Bridgman systems car be con- 
nected mathematically through the thermal boudnary conditions and the 
flow structure can then be thought nf In terms of a continuous transi- 
tion between the two systems. The fruitfulness of this approach is 
demonstrated on Figure 14 where a sample bifurcation diagram and flow 
fields for the Bridgman system are displayed. Notice that the patterns 
of the streamlines and the existence of the limit point Ra^ are quail ta 
tively unchanged between the Rayleigh-Benard and Bridgman systems. 


Effect of Buoyancy-Driven Convection on Melt/Solid Interface Shape and 
jjadiaT Segregation in Vertical Bridgnan Growth . 

Melt/solid Interface shape and segregation In the crystal are the 
most easily measured indicators of the role of convection in melt growth 
systems (Kim et al . 1978, Wang and Witt 1980). Previous calculations of 
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Figure 14. Bifurcation diagram for vertically destabilzed 
Bridgman system. 
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convection in melt growth do not account for the shape of the melt/ sol id 
interface or for realistic heat transfer conditions and thus cannot sim- 
ulate experimental results. The numerical methods developed by Charg [A] 
make possible these calculations In a wide variety of small-scale growth 
systems . 

Chang's algorithm has been used at M.I.T. to study convection and 
mass transfer of a dilute component In thermally stabilized and destabilized 
vertical Bridgman systems; the stabilized system is the same as shown in 
Figure lib with the ampoule and furnace flipped over. These caluclations 
are aimed at modelling the Gallium-doped Germanic experiments of Chris 
Wang In the Materials Science Department at M.I.T. and will supply a basis 
for evaluating convection In other earth-bound experiments related to NASA's 
program on directional solidification. The results discussed below are 
taken from the more complete report [8]. Sample isotherms and streamlines 
calculated in the stabilized geometry are shown as Figure 15 for Pr * 0.01 
and a dimensionless growth rate Pe * 0.01 which corresponds to moving the 
ampoule at roughly 10 ym/sec. The melt/solid Interfaces In these calcula- 
tions are not flat, but slightly curved because of radial temperature gra- 
dients at the edges of the gradient region and because of convective heat 
transfer. Interfaces are shown In Figure 16 that correspond to the flow 
fields given in Figure 15. These curved interfaces together with the cor- 
vection in the melt cause radial segregation of the dopant in the crystal. 

Chang has calculated the solute distribution In a steady state analog 
of Bridgman growth using the velocity fields and interface shapes predicted 
from his convection calculations. Figure 17 shows the distribution of solute 

4 

along the solidification Interface for Rayleigh numbers up to 10 for the 
gallium-doped germanium system of Wang; the Schmidt number Is 10 and the 
segregation coefficient Is 0.1 . Almost no radial segregation was present 
in the absence of natural convection (Ra * 0) because the melt/solid Inter- 
face was sensibly flat ; see Figure 16. As the convection level was Increased, 

the delivery of solute to the Interface became uneven radially until at Ra » 

3 4 

1 x 10 the radial segregation of gallium was maximized. Ra * 1 x 10 

radial segregation is still appreciable and the much used boundary-layer 
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approximations for calculating dopant distributions are not applicable. 

The cencentratlon fields at higher convection levels tend toward a uni- 
form core, with regions of rapid change In dopant level near each boundary, 
this structure marks the start of the formation of boundary layers on 
these surfaces. The amount of radial segregation Is plotted on Figure 
18 as a function of Rayleigh number and segregation coefficient for Sc ■ 

10 . 


The sensitivity of melt/solid Interface shape and dopant segregation 
to such important parameters as the length of the gradient zone and the 
ratio of conductivities between melt and solid are established In [8]. 


Microscopic Modelling of Helt/Solld Interface Dynamics 

Morphological stability theory describes the tendency of a micro- 
scopically smooth solidification front separating a binary melt from its 
solid to develop undulations that lead to cellular and dendritic inter- 
faces. Although the onset of these Instabilities has been studied using 
linear stability analysis (Mullins & Sekerka 1964) and by weakly nonlinear 
analysis (Wollkind and Segal 1970), metallurgists still lack a theoretical 
understanding of the microscale behavior of highly deformed Interfaces. 

Our research has concentrated on efficient calculation of large amplitude 
deformations of a melt/solid Interface and on the effect of various system 
properties (temperature gradient, growth rate, material properties, and 
grain boundaries) on the structure of the interface. 

By adding the variation of the melting point temperature with the mean 
curvature of the solidification interface, Ungar and Ettounay [12] have 
developed a numerical method for calculating tne shapes of solidification 
Interfaces on the microscale. These numerical methods have been used to 
study a simplified model for directional solidification where the temperature 
gradient Is assumed to be linear throughout both melt and crystal and un- 
effected by the morphology of the melt/solid Interface. The assumptions 
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leading to this model are often valid and have been made by other re- 
searchers (Hunt 1979, Langer 1980). The calculations shown here are 
for parameter values close to those appropriate for the Pb-Sb system 
studied theoretically in Coriel 1 and Sekerka (1980) and experimentally 
by Morris and Winegard (1965). When the interface shape was taken to 
be periodic with a given wavelength the linearized results of Mullins 
and Sekerka (1964) were recovered and steady state interfaces with 
finite-amplitude cellular morphologies were calculated. 

The evolution of these interfaces with decreasing temperature gra- 
dient G (increased constitutional instability) is depicted on Figure 
19a for a dimensionless growth rate Pe = 0.8 by the maximum deflection 
of the solidification from; a sample shape is shown as Figure 19b. The 
shape family first evolves subcriticaHy (to large temperature gradients) 
and is unstable, ht a critical value of G the shape family turns toward 
smaller values of |G| and the cellular forms regain stability and should 
be observable. The characteristic "molar" shape of the melt/solid inter- 
face was reported by Holmes and Gatos (1981). Increasing the growth rate 
changes the structure of the shape family from sub- to super-critical and 
results in an "s-shaped" shape family with two regions. of stable inter- 
faces; see [12] for details. 

Introducing a grain boundary in the melt/solid interface changes the 
structure of the constitutional instability. Analytical and numerical 
results for the one-sided model discussed above [13] show that the grain 
boundary ruotures the bifurcation between the planar and cellular inter- 
faces. When the bifurcation point is subcritical , adding the grain boundary 
leads to loss of stability at a limit point, as shown on Figure 19. When 
the original transition is super-critical (the case for high growth rates), 
adding the grain boundary leads to a smooth transition to a highly curved 
interface with no change in stability. These results are contrary to pre- 
vious linear analyses (Coriell and Sekerka 1973a, 1975b) for grain angles 
near 90° which predict that the pressure of the grain does not affect the 
point of exchange of stability. Ungar has shown that, for temperature gra- 
dients near the critical value for the cellular instability, the steady 
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interface shapes are singular functions of the grain angle a and a regular 
perturbation analysis in this parameter, as performed by Coriell and Sekerka, 
fails. The correct singluar analysis predicts the critical temperature 
gradient for the existence of the limit point mentioned above varies as 
(o)2/3. The details of this analysis will be presented in [13]. 
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Solutal Convection and its Effects 
on C r ystal Growth in Binary and Pseudo-Binary Systems 
with Large Liquidus-Solidus Separation 

Principal Investigator a: Robert A. Brown 

and 

August F. Witt 


SUMMARY 

This report details the first year of a theoretical and experimental 
study into solutal convection and its effects on crystal growth and seg- 
regation in binary and pseudo-binary systems with large liquidus-solidus 
separation. Significant progress has been made toward the qoals of this 
research. Theoretical analyses of melting and equilibration on melt and 
solid have yielded new insights into the first stages of crystal growth 
of binary compounds; a simple, one-dimensional model of directional so- 
lidification has been used to examine transient behavior in non-dilute 
systems; and finite-element methods have been developed for simultane- 
ously calculating thermal- and solutal-driven convection, temperature, 
and concentration fields and the melt/solid interface shape in vertical 
Bridgman growth. Experiments at MIT have focused on the melt growth of 
PbSnTe which has been carried out successfully in both a liquid-encap- 
sulated Czochralski and vertical Bridgman configurations. Research in a 
parallel collaborative study in the Metallurgy Division of the Centre 
d' Etudes Uucleaires de Grenoble has centered on the growth of GeSi in a 
Bridgman furnace. 


SCOPE OF PROGRAM 


During solidification, conventional doped semiconductor systems ex- 
hibit growth and segregation behavior which is in quantitative agreement 



with established theoretical concepts for mass transfer and fluid flow. 
There exists mounting evidence, however, that systems with large liqui- 
dus-solidus separation exhibit segregation and morphological character- 
istics that cannot be reconciled with the conventional pictures for 
convection in the melt which rely on thermal gradients to create den- 
sity gradients for driving motions of the melt. In non-dilute systems 
density differences caused by the rejection of solute during solidifi- 
cation give rise to buoyancy-driven flows that are not accounted for in 
segregation theory. The goal of our research and of the parallel re- 
search programs led by Or. J. J. Favier of the Metallurgy Division of 
the Centre d'Etudes Nucleaires at Grenoble is to understand the role of 
solutal -driven convection in crystal growth of materials with large sep- 
aration of the liquidus and solidus curves. Exoerimental and theoretical 
programs are underway at both institutions and the collaboration has been 
fueled by trips by E. D. Pourret and R. A. Brown to Grenoble and by J. J. 
Favier to M.I.T. 

The experiments focus on separating out the effects of solutal- 
driven convection in vertical Bridgman growth of binary and pseudo-bi- 
nary cyrstals. To do this, two systems are being considered that 
should exhibit qualitatively different flow behavior. The French group 
is studying (Rouzaud 1982) Bridgman growth of GeSi (Ge - 5 % Si alloy) 
with the crystal on the bottom. In this configuration, the mean 
(cross-sectional average) axial temperature increases with distance 
into the melt and is stabilizing with respect to buoyancy-driven flows. 
During growth, the heavier component (silicon) is rejected at the so- 
lidification front and thus the mean axial composition profile is also 
stabilizing. Hence, the only convective motions expected in this sys- 
tem are those caused by radial gradients of temperature and composition 
induced by curvature of the melt/solid interface. The furnace recently 
designed in Grenoble for these experiments (Rouzaud 1982) will have 
small radial gradients (theoretically zero) over a nine centimeter 
Tength of the furnace while maintaining a 50°C/cm axial gradient. 

Experiments at MIT have centered on the growth of PbSnTe, a 
pseudo-binary compound, which when grown in the thermally stable Bridg- 
man configuration rejects the less dense component (SnTe) at the inter- 
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face and so leads to a density profile in the melt that is unstable. 
Much more rampant convection is expected here than in the GeSi system. 
Progress toward these experiments is described below. 

The two experimental programs are coupled together by theoretical 
calculations boing developed at MIT for the steady convection in pro- 
totypes of these experiments. The calculations are based on the fi- 
nite element algorithm developed by Chang (Chang and Brown 1982) for 
computing simultaneously the temperature, velocity, and pressure 
fields in the melt, the temperature field in the crystal , and the 
shape of the melt/solid phase boundary for a single-component melt. A 
great deal of thermo-physical data is necessary to successfully model 
eithet of these experimental systems. Only in vhe GeSi system is the 
data available (Dismuke and Kim 1974). A thorough literature survey 
of the PbSnTe system (Adornato 1981) has exposed large gaps in the 
knowledge of this system, especially in the melt. For example, such 
routine data as the thermal conductivity of the melt has not been mea- 
sured. Without this data, the matching of calculations to experiments 
will be extremely tenuous. 


RESEARCH SUMMARY 

Basis Factors Control ing Seeded Melt Growth of Concentrated Alloys . 

experiment design and execution: 

Basic studies of semiconductor crystal growth from the melt have 
in the past been mainly concerned with the behavior of dilute (doped) 
alloy systems. Recently, however, focus is being placed on more, com- 
plex pseudobir.ary systems (HgCdTe, PbSnTe, GeSi, GalnSb, etc.) in 
which the 1 iquidus-sol idus separation becomes larger. A large li- 
quidus-sol idus gap may lead to major complications for experiment de- 
sign and execution since the growth interface temperature is a pro- 
nounced function of composition and thus subject to change by several 
tens of degrees during the course of a growth experiment. Complica- 
tions in solidification of such crystals arise primarily from uncer- 
tainties concerning the growth interface temperature and the location 
and shape of the crystal-melt interface within a given thermal gradi- 
ent region. A seeded melt-growth experiment in a Bridgman-type con- 
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figuration can be conducted either by establishing melt-seed contact or 
by partial melting of a single crystal ingot and subsequent regrowth. 

This last procedure was adopted, for example, in growth experiments 
under reduced gravity conditions (Witt et al. 1978). In both approaches 
the actual location of the seed-melt interface within the thermal gradi- 
ent region is controlled by the prevailing compositions of the contact- 
ing phases. The concentrated alloys, the growth interface temperature 
at the end of the back-melting (seeding operation) as well as during 
equilibration assumes a value which is a priori unknown but is some- 
where between the solidus and liquidus temperatures, which in such sys- 
tems may differ by as much as 100°C. The consequences of this uncer- 
tainty in the value of the interface temperature are apparent in Fig. 1 
which provides a schematic axial temperature distribution in a Bridgman 
growth configuration. Accordingly, if the growth interface temperature 
is close to the liquidus temperature, the indicated location of the am- 
poule is adequate; if, however, the interface temperature approaches 
the solidus temperature, the depicted growth configuration will lead to 
failure of the experiment. 

A further experimental complication arises from the requirement of 
growth interface "planarity" which is considered essential if radial 
compositional homogeneity in the absence of convection is to be achieved. 
While it is feasible in Bridgman configurations to control heat flow for 
the establishment of growth interface planarity at a given temperature 
(Fig. 2), it is beyond the state-of-the-art to achieve interface plan- 
arity at a given temperature (parallel heat flow) over a temperature 
range of about 100 °C. A solution, in principle, can be found if, for 
transient melting and growth, the time (composition) dependence of the 
growth interface temperature is theoretically or experimentally deter- 
mined and heat flow conditions are continuously adjusted so as to main- 
tain interface planarity under transient growth conditions. 

The complications associated with seeded crystal growth of systems 
with large 1 iquidus-sol idus separation are related to diffusive mass 
transfer within and between the liquid and solid phases involved. Since 
mass transfer in the liquid phase is strongly affected by convection, 
any meaningful growth experiment must be carried out either in the ab- 
sence of convection or under conditions of controlled convection. 
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Figure 1. Uncertainty of interface location in systems 
with large 1 iquidus-sol idus separation. 
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Figure 2. Schematic of Bridgman growth of non- 
dilute binary crystal. 
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Melt-Seed Contact and Equilibration . 


In a binary system, such as the one schematically indicated in Fig. 

2, there is only one set of compositions, and C^, at a given temper- 
ature for which the liquid phase and the contacting solid phase are in 
thermodynamic equilibrium; for example, and at 

T l(l) . This particular configuration must be considered as hypothetical 
since it is not normally realized in seeded crystal growth experiments. 

A whole range of concentrations can be encountered for both the solid 
and the liquid phase at the beginning of the equilibration period, i.e., 
when the melt contacts the seed. For example, the seed composition is 
often identical to that of the melt. Any configuration departing from 
equil ibrium yields instabilities, and diffusional transfer takes place 
through the liquid-solid interface. In the subsequent analysis, the 
interface temperature at the time of arrest is taken to be the liquidus 
temperature at the concentration of the melt so that there are no addi- 
tional phenomena due to a supercooled melt. The range of composition 
chosen for the seed should be limited to regions I (Fig. 2) and II since 
in region III the seed would melt to an extent determined by the thermal 
gradient in the furnace and ^he resulting increase in the melt concen- 
tration would decrease its liquidus temperature. 

The free- energy curves drawn at T. \ for the liquid and the solid 

^ l U / / Q \ 

phases are shown in Figs 3 and 4. A solic at C^' ' in region I (Fig. 

3) cannot form a liquid at temperature T^jj. trar) sition from C^(o) 

to a value of lower free energy can occur by diffusion of the solute 
through the interface and yield an equilibrium configuration at the in- 
terface. If diffusion in the liquid is much faster than in the solid 
(D l >> D s ) or in the presence of intense convective mixing, no boundary 
layer build-up in the liouid phase occurs (assuming infinite length to 
eliminate end effects) and the transformation is essentially isothermal 
(path <1> in Fig. 3). As the ratio a = D^/D^ increases, a solute boundary 
layer will form ahead of the solid-liquid interface. The interface comp- 
osition in the liquid phase decreases, which induces a simultaneous increase 
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in the interface temperature according to the shape of the liquidus curve 
(path 2 in Fig. 3). Subsequently, the interface moves toward the liquid 
and solidification occurs durinq equilibration. 

In region II of the phase diagram (Fig. 4) the solid can form a 
liquid, which lowers its free energy toward an equilibrium configura- 
tion. The extent to which the seed melts is limited by the mass trans- 
fer through the interface, the thermal gradient in the furnace and th;. 
dilution effect resulting from melting of the seed. If >> 0^, as 
for the case of the seed composition in region I, the process is iso- 
thermal. As a = D<-/D l increases, the interface temperature drops due 
to the solute diffusion layer in the liquid phase and equilibration 
occurs through melting of the seed. Melting of the seed is expected 
to be a rapid, diffusionless process; the resulting dilution of the 
melt can then resul t In partial resolidification. In region III, an in- 
crease in the melt composition due to melting of the seed will induce 
further backmelting until equilibrium is reached. 

The above thermodynamic considerations are useful to predict if 
the equilibration step is accompanied by solidification (C s ‘ 0/ in region 
I) or by melting (C^ 0 ^ in region II); it, however, fails to give quan- 
titative. data for the evolution of the concentration and temperature at 
the liquid-solid interface and in both bulk phases. Ghez end Small ( 

1981; also see Small and Ghez 1979) have solved the problem of mass trans- 
fer involved in the equilibration process for the formulation of hetero- 
structures by liquid phase epitaxy. An analysis of equilibration during 
seeding in melt-growth differs from their treatment for the thermal field 
governing the process and for the simultaneous change of the bulk melt 
concentration as mentioned earlier. It then Involves a solution for the 
coupled heat and mass transfers in both phases, the interface conditions 
being tied to thephase diagram. Such a Stephan problem has not yet been 
solved analytically. However, the expected results can be qualitatively 
assessed. In Fig. 5, regions I and II are referred to as "stable" and 
"unstable" respectively with regard to the morphological stability uf 
the interface in each configuration. For convenience, one introduces 
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Figure 5. Compositional inhomogeneities due to 
equilibration in the crystal once re- 
growth is initiated. 
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the ratio of the initial bulk seed composition to the initial bulk melt 
composition, p * (C^/C^ 0 ^). From the previously stfted assumption, 
the interface temperature is the liquidus temperature ”, ^ and the 
system Is maintained In a positive stabilizing thermal gradient. The 
solid then is in a stable configuration In region I, ifp<k, k <1 or 
p > k, p > I, since the temperature of the interface is below the 
solidus temperature for the range of concentration between C<^°^ and 
(C s (eq) is ti-e interface composition in the solid one* equili- 
brium is reached.). In the liquid phase, the temperature is at the 
liquidus temperature at the interface and above it further from the 
interface. According to the consitutional supercooling criteria (Rut- 
ters et al . 1953), ' ie system is in a stable configuration. 

In region II, the system is subject tc morphological instabilites aid 
steep temperature gradients are required to maintain stability f p > k, 
k < 1 or p < k, k > 1. Initial conditions such that 

p < k, k < 1 or p > k, k > 1 , corresponding to the bulk seed composi- 
tion in region I are then favorable since it prevents regions of super- 
cooled liquid or superheated solid. The concentration profiles ex- 
pected in the crystal after equilibration and are presented in Figure 5 
for each initial condition, discussed above. If the seed is initially 
in region I, a transition region of uncontrolled grow;h is expected in 
which the ccncentration increases rapidly, due to uncontrolled solid- 
ification. If the seed Is initially in region II, a decrease in the 
concentration is expected in the region of uncontrolled growth, due to 
melcinc, and eventual partial resolidification, as discussed earlier. 

It is important to note that the exact concentration profiles in the 
region of uncontrolled growth are strongly dependent on the level of 
convective mixing in the melt. 


Crystal Growth and Characterizations of PbSnTe . 

Experimental research has focused on techniques for the melt 
growth of single crystal PbSnTe and for the characterization of its 
segregation behavior. Both the vertical Bridgman and Czechralski 
techniques have been used successfully to grow PbSnTe. In the [’ridgman 
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technique, a precast charge of PbSnTe is loaded into a quartz ampoule 
that is tapered on the bottom to seed the crystal . The ampoule is po- 
sitioned in the vertical furnace described by Witt (1981) and melted. 
Directional solidification is driven by lowering the ampoule from the 
bottom of the furnace and the onset of solidification is by hetero- 
geneous nucleation caused by the temperature at the bottom of the am- 
poule dropping below the solidus temperature of PbSnTe melt. When 
the included angle of the tapered ampoule exceeded 60°, the grains in 
the poly-crystal grew toward the side of the ampoule and left a single 
crystal in the cylindrical portion of the ampoule. Heterogeneous nu- 
cleation was suppressed completely by seeding the melt with a single 
crystal of PbSnTe with < 1 00 > orientation. Results of experiments 
with seeded growth are described in [3] and highlighted below. 

The PbSnTe crystals grown by either the Czochralski or vertical 
Bridgman methods were soft and cleaved easily alcng the [100] direc- 
tion. Special care was needed to minimize mechanical stresses in the 
crystal induced during preparation of the crystal; for example, both 
ultrasonic cleaning and high-speed mechanical polishing had to be 
avoided. Good quality polishing of slices cut with a low-speed wire 
saw was obtained by sequentially: 

1. Polishing with 5 pm Al^O^ powder on a glass plate. 

2. Polishing with 0.3 yin AljO^ on a felt pad. 

3. Polishing with 0.05 ym Al^O^ on a felt p*d. 

4. Polishing with a Clorox-Syton mixture (one-to-one by volume) 

on a felt pad. 

The quality of the crystal surface was conti oiled by Normarski inter- 
ference microscopy. This procedure yielded perfectly polished surfaces 
that were the source material for studies of differential etching. 

Striations in the crystals grown by the Czochralski method were 
revealeo using the Norr electrolytic etch (6 volts for 4 minutes). 
Compositional traces made by the electron microprobe confirmed that 
this etching differentiates between PbTe and SnTe rich regions. This 
etchant has failed to repeal the concentration variations caused by 
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Peltier pulsing in the vertical Bridgman configuration. T't.e cause of 
this failure is not completely understood, although one possible ex- 
planation is the lack of sensitivity of the etchant to small concen- 
tration fluctuations. This point is being investigated further. 

Since the large 1 iquidus-sol idus separation in the PbSnTe system 
caused a pronounced composition dependence of the growth on the inter- 
face temperature, it was of interest to analyze both the growth and 
segregation behavior in the region of the seed-melt interface. We 
have done this using composition analyses (wt. % Sn) based on electron 
microprobe scans for crystals grown in the vertical Bridgman system; 
complete results are reported in [3], The axial composition profiles 
are shown in Figure 6 for the two crystals (A and B) whose growth 
parameters are summarized in Table I and indicate the e/istance of a 
"transition region" between the seed and the crystal grown at constant 
lowering rate; in this region the Sn concentration changes from that 
of the seed to the melt composition resulting from dilution of the 
charge by the partial melting of the seed. Using microprobe scans 
with 11 ym spacings this transition region is found to have a thickness 
of 400 and 150 ym for the crystals A and B, respectively. 

The origin of the region of uncontrolled growth was deduced from 
the length of seed meltback (7 mm in both crystals) and the composition 
of the initial region of growth. Upon meltback of the seed, with the 
crystal-melt interface at the liquidus temperature, a melt region is 
generated in which the Sn concentration is significantly less than in 
the bulk melt. Since, during the ensuing equilibration period the Sn 
concentration in this region increases because of the diffusion and 
convection of Sn from the bulk melt, the liquidus temperature will in- 
crease and cause crystal growth at the melt/solid interface without 
pulling of the crystal. The magnitude of this uncontrolled growth is 
a function of the compositions of the seed and bulk melt, the extent 
of meltback, the length of melt, and the degree of compositional equil- 
ibration of the melt. In terms of these parameters, the difference in 
the widths of the uncontrolled growth regions in crystals A and P is 
attributed to the lower Sn concentration in the seed of A which results 
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in a lower interface temperature during meltback and more extensive 
uncontrolled growth during equilibration of the melt, than for crystal 
B. 


It is significant that, within the error limits of the analysis, 
the composition of the initial portion of controlled growth is the 
same for both crystals. The short equilibration period of only 40 
minutes for crystal A points to the fact that convective mixing was 
the prevalent mechanism for obtaining compositional uniformity of the 
meit. As is seen in Figure 6, the axial composition profiles for the 
first 3 mm of growth in crystals A and B are identical within error 
bounds for the microprobe. Tin also was found to be distributed 
uniformly across the radius of each crystal as shown in Figure 7. 
Strictly, this result is inconclusive since a uniform Sn concentration 
can result from either a flat melt/crystal interface in the absence of 
convection or from a curved interface coupled with so rampant convec- 
tion that a thin and uniform boundary layer forms along the interface 
and isolates tne crystal from radial segregation. When taken together, 
the radial and axial segregation behavior point to convection in the 
melt, especially adjacent to the interface. 


Modelling of Growth and Convection in Non-Pi lute Binary Systems . 

We have pursued two approaches in our efforts to model the di- 
rectional solidification of binary crystals. In the first, the tran- 
sients in growth rate and composition caused by changes in the ampoule 
pull rate or the temperature of the surroundings are being studied 
using a simple, time-dependent one-dimensional model for solidification. 

In the second set of calculations, the more detailed finite-element cal- 
culations are being used to determine the fluid flow and melt/solid 

interface shapes in a steady-state prototype of the Bridgman system. j 

The motivation and approach to each of these analyses is laid out j 

below. i 

i 

i 
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Figure 7. Radial solute segregation in Pb Q ,Sn 0 ^Te 
crystals A and B grown at V = 0:5um7Sec. 



Figure 8. Distribution of furnace temperature and 
heat transfer coefficient for the static 
ampoule. 
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One-Dimensional Modelling of Transients in Directional Solidifi- 
cation . 

When the melting point of the crystal is strongly dependent on 
the composition of the melt at the melt/solid interface the rate of 
equilibrium of the growth velocity to changes in the pull rate depends 
on the rate of equilibration of the composition field in the melt as 
well as on the time for equilibration of the thermal field. In fact, 
for typical binary and pseudo-binary melt growth systems, e.g., 

HgCdTe, PbSnTe, and GeSi, the time scale for mass diffusion is much 
larger than the one for heat and it controls the rates of any transi- 
ents in the growth system. 

We have developed a one-dimensional (in the axial direction) 
model of heat and mass transport in a Bridgman system that contains 
all the gross features of the growth when mass transfer is controlled 
by diffusion. Even in this case the resulting equations describe a 
moving boundary problem for the location with time of the melt/solid 
interface, the concentration profile in the melt and the temperature 
profile in both phases. The boundary immobilization method first 
developed by Landau (see Ettouney and Brown 1982) is used to discre- 
tize this equation set on a finite difference mesh fixed in space and 
the system of nonlinear ordinary differential equations that results 
are solved using a variable step-size Runge-Kutta integration scheme. 
Applications of this scheme >.o the growth of PbSnTe are discussed below. 

When the thermal and m---ss diffusivities are much different, i.-. : 
when the Lewis number Le s VMj is not close to unity, the equation set 
is numerically stiff and the Runge-Kutta scheme is forced to march in 
very small time increments in order to maintain a prescribed level of 
accuracy. The levels of accuracy and efficiency in the integration 
scheme necessary for repeated detailed simulations of real growth systems 
have been reached by solving the original set of 60-100 coupled differen- 
tial equations by an implicit (Gear 1976) integration scheme that is numer- 
ically stable for larger time steps. This technique requires the solution 
of a set of nonlinear algebraic equations at each time step, but allows 
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Table 1. Conditions for Growth of PbSnTe. 
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much larger time steps for numerical stability. Comparisons between the 
Runge-Kutta and implicit integration schemes show the implicit methods 
to be most efficient. 

Calculations are reported in [3] for the growth of Pbg gSng 2 Te in a 
vertical Bridgman system under convectionless conditions. The numerical 
simulations are carried out for an ampoule with 0.5 cm radius being pulled 
initially at a growth rate of 0.5 ym/sec. The temperature profile along 
the wall of the furnace is shown in Figure 8 along with the dimensionless 
heat transfer coefficient between the furnace and the ampoule. The tran- 
sient that appears at the beginning of pulling was simulated by starting 
the numerical simulation with the temperature profile for a stationary and 
equilibrated ampoule. As the ampoule moves the temperature changes toward 
the solidus temperature as illustrated by the five temperatures shown on 
Figure 9. These curves correspond to the following steps ir growth: 

<0> the equilibrated system (t = 0 sec). 

<1> displacement of the ampoule by 1 cm: 

V = 0.5 ym/sec and t = 2 x l(r sec. 

<2> displacement of the ampoule by 2 cm: 

V = 0.5 ym/sec and t - 4 x 10^ sec. 

<EQ> equilibration of ampoule aftergrowth 

rate arrest at t = 4 x 10* sec: V * 0 

and t = 6 x l(r sec. 

<3> displacement of the ampoule oy 3 cm 
after resumption of pulling at t = 

6 x l r 4 sec; V *C 5um/sec and t = 

8 x 10 4 sec. 

The 3.5°C rise in interface temperature during the period of equilibration 
in the middle of the run corresponds to a decrease in the SnTe concentration at 
the interface as it's gradient was flattened by diffusion toward the bulk of 
the melt. This point is illustrated on Figure 10 where the solute distributions 
corresponding to each of these points in the experiment are shown. During the 
initial period of pulling, the solute diffusion layer grew at the interface. 

Upon stopping the ampoule, solute diffusion away from the interface dominated 
over the small amount of SnTe rejection caused by the uncontrolled growth and 
the solute composition decreased the interface. After the resumption of pulling 
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Figure 9. Variation of axial temperature profile 
during growth c*f Pb Q g Sn 0 „Te crystal at 
V * 0.5 ym/sec . CQfvesxorrespond to 
specific times during the transition to 
steady growth and equilibrium in the melt. 


Table II. 

Interface Temperatures During Transient 
Growth of PbSnTe . 

Point 

Time (10^ sec) 

Temperature T(h(t))(°C) 

<0> 

0 

900.8 

<1> 

2 

895.3 

<2> 

4 

891 .5 

<EQ> 

6 

895.0 

<3> 

8 

890.0 
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Figure 10. Concentration of SnTe in the crystal during 
initial transient for V = 0.5 ym/sec, a 
short equilibration period, and a resumption 
in pulling. Concentration profiles in the 
melt are also shown for several interface 
positions. 




Figure 11. Interface position as a function of time for 
the sequence of pulling rate changes described 
for Fig. 10. 
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the solute profile started to develop a new gradient, but followed a 
different transient than at the onset of the run, because a non-planar 
solute profile existed at t * 6 x 10* sec. 

The time history of the ratio of growth rate to pull rate R(t)/V 
is shown on Figure 11, with markers to show the different steps of the 
simulated run. Steady state is not reached during the first two centi- 
meters of ampoule translation, as shown by the failure of R(t)/V to reach 
unity. End effects in the concentration profile make it impossible for 
the growth rate to ever reach the pull rate. The coupling between the 
growth rate and the transient transfer of solute is the major difference 
between growth of dilute and non-dilute alloys. To make clear the im- 
portance of this coupling, we calculated the initial transient for the 
growth of a dilute alloy with thermal properties identical to those used 
for PbSnTe, but with the melting point fixed at 901°C. This case is 
plotted on Figure 12 as the dashed curve and shows the much more rapid 
thermal equilibrium that is possible in the dilute system. 

The one-dimensional model for directional solidification described 
in [3] is also being used to model transient behavior in other binary 
(GeSi) and pseudo-binary (HgCdTe) systems with large 1 iquidus-sol idus 
separation; these results will be reported in [4]. 


Finite Elem e nt Calculations of iolutal Convection 

Melt/solid interface shape and segregation in the crystal are the most 
easily measured indicators of the role of convection in melt growth systems. 

New numerical methods developed by Chang and Brown (1982) make possible the 
solution of the full free-boundary problem governing small-scale crystal 
growth systems. Chang's algorithm has been used at M.I.T. to study convection 
and mass transfer of a dilute component in thermally stabilized and destabilized 
vertical Bridgman systems; results of these calculations are described in Chang 
and Brown (1982). 



Peter Adornato has modified this algorithm to include the coupling 
between the solute, temperature and velocity fields caused by solutal- 
driven convection and by the phase diagram. These modifications are straight 
forward for calculations at low solutal Rayleigh numbers and for completely 
stabilized systems like GeSi . In these calculations, the convection driven 
density differences is weak, solutal boundary-layers are broad, and the mod- 
erate separation of liquidus and solidus leads to small displacement of the 
interface. Our calculations to date have focused on the GeSi system being 
studied in Grenoble. 

For systems with a large gap between liquidus and solidus (HgCdTe for 
example) the deformation of the interface caused by radial s']regation may 
be severe. In systems that are unstable with respect to either the axial 
solute or temperature profile convection will be severe and thin solute 
boundary layers will form along the ampoule and the interface. In the next 
year, the multi-grid method developed by Chang (1982) will be adapted to 
attempt to calculate these cases. 
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Analysis of Polymer Melt Processing Operations 

Principal Investigator: David K. Roylance 


RESEARCH ABSTRACT 


This report outlines the equations which govern the 
nonisothermal flow of reactive fluids, and describes the 
means by which finite element analysis can be used to solve 
these equations for the sort of arbitrary boundary condi- 
tions encountered in industrial practice. The performance 
of the computer code is illustrated by several trial prob- 
lems, selected more for their value in providing insight to 
polymer processing flows than as practical production prob- 
lems. Although n good deal remains to be learned as to the 
performance and proper use of this numerical technique, it 
is undeniably useful in providing better understanding of 
today's complicated polymer processing problems. 
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Introduction. 


Finite element analysis offers great promise for reducing 
the empiricism often found in polymer processing design, 
since it is well suited for modeling the complicated bounda- 
ry conditions and material properties encountered in indus- 
trial practice. Although the method is now well accepted in 
structural stress analysis, its use in fluid transport situ- 
ations is less widespread. We have sought to exploit some 
of the many advantages of the finite element method in 
polymer fluid processing analysis, and this paper describes 
some of our work in chemorheology . 

Reactive flows have several advantages in polymer proc- 
essing in comparison with more traditional melt-flow tech- 
niques. Perhaps the most significant of these is the energy 
savings which result from the elimination of the several 
melting stages necessary in such technologies as extrusion 
or injection molding. However, these advantages are offset 
to some degree by the complexity of the process, which ren- 
ders the empirical approach to process development difficult 
in the extreme. 

The flow velocities in such processes are governed by the 
fluid viscosity, which is a strong function of temperature 
and molecular weight. The temperature in turn is c 'ected 
by the viscous dissipation and the heat released or consumed 



by the reaction, and the reaction rate is also a strong 
function of temperature. All of these variables interact 
and change in such a way as to make an intuitive grasp of 
the process almost impossible, and there is obviously an 
advantage to being able to provide some sort of mathematical 
or numerical simulation of the process. 

The finite element scheme to be described below is very 
useful for obtaining numerical solutions for reactive flows 
with arbitrary boundaries, and such a technique is well 
suited for detailed analysis of real industrial processes. 
However, we argue that the greatest value of the method may 
not be in these detailed calculations, but in the degree to 
which the method can enhance the process designer's intui- 
tion. For this latter purpose, it is often sufficient to 
run only very small ar.d inexpensive trial Problems; these 
■ran elucidate the manner in which the various >rcblem param- 
eters interact, so that the designer develops a much 
improved "feel" for the problem. It is likely that most 
real advances in processing technology will continue to be 
made more by Edisonian innovation than by detailed math- 
ematical calculations. But since today's processes have 
become so complicated, such a technique as the finite ele- 
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ment method can be a powerful adjunct to intelligent 
intuition . 

Theoretical background. 


The equations which govern the nonisothermal flow of a 
reactive fluid are derived in several texts on transport 
phenomena and polymer processing (e.g. refs. 1,2). Regard- 
ing velocity, temperature, and concentration of unreacted 
species as the fundamental variables, the governing 
equations can be written as: 

/>(au/at + u« 7 u] ■ - v p + i 7 *u 

/>c[aT/at + u* 7 T] * Q + k 7, T 

[aC/at + u* 7 C] « R + D 7 *C 

(A list of nomenclature is attached.) The similarity of 
these equations is clear. In all cades, the time rate of 
change of the transported variable (velocity, temperature, 
or concentration) is balanced by the convective or flow 
transport terms (e.g. u* 7 C) , the diffusive transport (e.g. 
D VJ C), and a generation term (e.g. R) . 
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The analyst seeks expressions for the space- and time- 
dependent velocities, temperatures and concentrations which 
satisfy these equations and also the problem's boundary con- 
ditions. It is generally the boundary conditions which make 
real problems intractable: even if one were able to 

describe the boundaries mathematically, the resulting 

expressions would not likely be amenable to closed-form sol- 
ution. In addition, many of the "constants" in the above 
equations are often nonlinear functions of the problem vari- 
ables. In reactive polymer processing, one might encounter 
such expressions as the folloving: 

i * ., 0 * n - ’exp^/RgTjexp^p) (mw 3 4 ) 

Q « (,/2)(*:i) + RUH) 

R - -k m exp(E 2 /R g T)C m 

It is cleer that all of these expressions, taken 
together, constitute a mathematical situation which must be 
approached with caution. Even though it is not overly dif- 
ficult to incorporate them into a numerical scheme, which we 
have done, it is important to proceed slowly enough to 
develop the proper experience in the computer code's behav- 
ior before tackling full-blown problems. 
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Comouter Model. 


We have sought to develop a finite element code which is 
able to predict polymer fluid velocities, stresses, temper- 
atures, and degrees of chemical conversion in a variety of 
flow geometries and for a variety of fluid material proper- 
ties. Space limitations prohibit our listing here the full 
derivation of the finite element equations from the above 
difterential equations, but there exist several well devel- 
oped means by which this may be done. The reader is 
directed to standard texts (3,4) for a more complete 
description, and we will just state here that we employ the 
Galerkin method of weighted residuals, together with 
isoparametric mapping and interpolation, to replace the dif- 
ferentials by integrals which can be evaluated numerically 
over small subregions ("elements"). The results of these 
numerical integrations are then assembled into a set of 
simultaneous algebraic equations which can be solved numer- 
ically. 

The salient features of our code can be listed briefly as 
follows: (1) Velocity, temperature, and chemical conversion 
are taken as nodal unknowns, so that coupled incompressible 
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viscous flow and diffusive-convective heat and mass trans- 
port may be modeled. (2) Incompressibility is enforced by a 
"penalty" formulation employing selective reduced inte- 
gration. This approach requires the use of double precision 
arithmetic with a concommitant reduction in the amount of 
available core, but it has several programming advantages 
which usually outweigh this drawback. (3) The code is 
developed primarily for plane and axisymmetric flows. We 
have coded a three-dimensional capability, but generally 
feel that the expense of running three-dimensional problems 
is not justified for most of our modeling research. (4) The 
code includes several models for the effect of shear rate, 
temperature, pressure, and chemical conversion on the fluid 
viscosity. These nonlinear models have not yet been 
researched extensively, however, and the exploration of the 
interative schemes needed for their proper use constitutes a 
major goal for future research. We have also coded a capa- 
bility for viscoelastic fluid effects (ref. 5), but 
currently feel that this difficult aspect of polymer flow is 
being researched satisfactorily by other workers. 
(5) Convective transport of heat or chemical species can be 
handled either by a conventional Galerkin treatment or by 
the convenient but still controversial "optimal upwinding" 
approach. (6) The code can treat transient problems by 
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means a two-point "theta-method" time-stepping algorithm. 
The dynamic algorithm is also useful in nonlinear problems, 
in which the final fluid state may be approached dynamically 
from an estimated initial state. (7) The code is capable of 
a variety of iterative treatments of nonlinear problems, 
including Newton-Raphson iteration and incremental load 
methods. 

Some additional discussion is warranted concerning the 
treatment of convective effects beyond what has been men- 
tioned in item (5) above. Momentum convection (/»u« v u) is 
generall negligible in comparison with the viscous terms 
due to the high viscosities of polymer fluids, but the 
convective terms tend to dominate the energy and mass trans- 
port equations due to the generally low thermal and mass 
dif f usivities. The programming of the convective terms pre- 
sents no special problems in the Galerkin approach beyond 
the need to stor-e and solve unsymmetric matrices, but it is 
well known that the presence of strong convective terms 
tends to create oscillations in the final solutions which 
can be large enough to destroy their value. This instabili- 
ty is related to the tendency of convection to produce large 
downstream gradients which the finite element grid cannot 
resolve. A largely ad hoc procedure known as "upwinding" 
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has been used in both finite element and finite difference 
work which seems to alleviate this problem by providing a 
greater numerical weight to the upstream portion of the ele- 
ment. Hughes (6) has published a very convenient means of 
upwinding, in which the sampling points in the numerical 
integration scheme are simply moved upstream an appropriate 
distance. We have made extensive use of the Hughes upwind- 
ing technique, but the reader is cautioned that this method 
is regarded as controversial by many workers. A provocative 
paper by Gresho (7) details many of the possible pitfalls in 
upwinding, and states a strong preference for grid refine- 
ment as the appropriate cure for convection-induced 
instabilities. 

Selected Numerical Results. 


Nonreactive entry flow. Figure 1 shows the streamlines 
for a 4:1 entry flow. Here a grid of 100 four-node linear 
elements was used to model the upper symmetric half of a 
plane capillary, and a fully-developed Poiseuille velocity 
was imposed on the reservoir entry as a boundary condition. 
The streamlines are identical with published experimental 
and numerical results, although the grid used here was not 
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intended to be fine enough to capture the weak recirculation 
which develops in the stagnant corner of the reservoir. 

The temperature contours for convectionless flow are 
shown in figure 2, which shows a hot region at the entrance 
of the capillary due to the combination of high viscous 
energy dissipation there and its distance from cool bounda- 
ries to which heat may be conducted. These isotherms are 
normalized on the maximum centerline temperature expected 
for Pois^uille flow in the capillary. 

The importance of thermal convection relative to con- 
duction is given approximately by the Peclet number 
Pe = UL/>c/k, where U and L are a characteristic velocity and 
length. Figure 3 plots the variation of temperature along 
the centerline for various values of the Peclet number, and 
it can be seen that the effect of increased thermal con- 
vection is to sweep the cooler upstream flow particles into 
the capillary, with a resulting lowering of the temperatu.es 
overall and a shift downstream of the hot spot near the 
throat. The relatively coarse grid used in this problem 
produced unstable Galerkin results for Peclet numbers higher 
than approximately ten, and the hlghe; degrees of thermal 
convection were computed using the upwinding formulation. 
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Further tests with refined grids should be completed to 
assess the accuracy of the upwinded solutions, although the 
plots in figure 3 appear reasonable. 

One-dimensional reactive flow. As a preliminary trial 
problem in our computations of reactive flow, we have 
studied a simple situation in which a fluid obeying 
first-order chemical kinetics moves at constant velocity and 
temperature in the positive x direction. Here only the 
mass-transport equation is operative, and it takes the sim- 
ple form: 

u(dC/dx) = -KC + Did’C/dx 1 ) 

This equation is solved easily, and for nonzero values of 
the diffusion coefficient D two boundary values for C must 
be specified. One of these is the initial concentration at 
the inlet, and the other requires a consideration of the 
outlet conditions. Here several possibilities exist, and we 
have studied the case in which the concentration of the cut- 
let reservoir is allowed to rise to meet that supplied by 
the flow; this is equivalent to specifying a zero concen- 
tration gradient at the outlet. For the case of negligible 
diffusion (D-0), the second-order term vanishes from the 
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above equation and the downstream boundary condition cannot 
be specified. The solution is then a simple exponential, in 
which the reactive species vanish according to first-order 
kinetics as they are carried downstream at constant 
velocity. 

Figure 4 shows the computed and exact predictions for 
reactive group concentration as a function of distance along 
the channel. The Galeikin values are nearly exact, but it 
is clear that upwinding leads to erroneous results in the 
small-diffusion case. The upwinding has introduced an arti- 
ficially high diffusivity and a zero concentration gradient 
at the outlet, and such artifacts must be considered as pos- 
sibilities when upwinding is used. 

Two-dimensional nonisothermal reactive flow. Figure 5 
shows the contours of constant conversion for a 
two-dimensional analog of the flow discussed in the previous 
section. Here again, a simple uncoupled problem is treated 
in which the material parameters are taken to be. independent 
of the solution variables, and in which the velocity condi- 
tions are prescribed. The concentration is taken to havt a 
fixed value at the inlet and a zero gradient at the outlet, 
as before. The two-dimensionality of the problem is con- 
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tained in two features: the velocity is taken to be 
parabolic, ranging frcm a maximum at the centerline to zero 
at the walls (a Poiseuille flow); and now diffusive heat and 
mass transport can occur in both the x and v direction. For 
the low dif fusivities shown, mass diffusi •- in the x direc- 
tion is negligible, as was demonstrated in the previous 
section. However, the concentration gradients in the y 
direction are substantial, so that diffusive transport in 
that direction is appreciaole even at D = 0.01. At 
D * 0.001, even the ; diffusion is negligible, so the con- 
centration contours simply represent a fluid which moves in 
the x-direction while reacting by first-order kinet : cs. The 
concentrations along the centerline are then identical with 
the D = 0.01 curve of figure 4. 


Figure 6 shows the contours of constant temperature which 
result from this flow (with D * .001), where the temperature 
boundary conditions were set to zero at the entry and along 
the top and bottom surfaces. The tempers ,ure gradient at 
the outlet was allowed to become zero, similar to the con- 
centration gradient. The results obtained fcr tne 
temperature field are of course dependent on the vaiim, cho- 
sen for fluid properties. To avoid using space here for a 
detailed discussion of the dimen, nal analyses used for 
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selecting these parameters, we will state simply that in 
figure 6 the viscous dissipation and reaction heat make 
approximately equal contributions to the internal heating 
(Brinkman and Damkoehler numbers both equal to three). 

Conclusions. 


The nuT'"ical model described above has a significant 
present ability to simulate a wide range of problems in 
polymer processing. \z the s?r.3 t:.ne, it is small enough to 
permit easy implementation in even rather small processing 
facilities, and for quick familiarization by process ^'sign- 
ers. We feel such a capability will lead to a significant 
advance in industry capability for process development and 
optimization. 


Nomenclature . 


c specific heat 

C concentration of reactive species 
D diffusion coefficient 
E, activation energy for viscosity 
E 2 activation energy for reaction 
k coefficient of thermal conduction 
k m reaction rate preexponential factor 
K o.- - i_ 1 1 reaction rate 
m reaction order 
mw molecular weight 
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n power-law exponent for viscosity 
p pressure 

Q internal heat generation 

R internal species generation 

R g gas constant 
T temperature 

u velocity or velocity vector 
p factor for pressure dependency 
% shear rate 
>i viscosity 
i/ C viscosity coefficient 
p density 
AH heat of reaction 
7 gradient operator 
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Figure Captions 

Figure 1 streamlines for nonreactive 4:1 entry flow, 
Newtonian fluid with imposed Poiseuille flow at inlet. 

Figure 2 - Contours of constant temperature for 

convect ionlesu entry flow, with heat generation by viscous 
dissipation only. 

Figure 3 - Entry flow centerline temperatures at various 
Peclet numbers. 


278 



Figure 4 - Degree of conversion along channel in 
one-dimensional reactive flow. 

Figure 5 - I soconversicn contours in two-dimensional 
reactive f^ow at two different mass dif fusivities, Galerkin 
calculations. 

Figure o - Isothermal contours in two-dimensional reac- 
tive flow, Galerkin calculation. Heat generation by viscous 
dissipation and reaction heating. 
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Figure 1 - Streamlines for nonreactive 4 : 1 entry flow, Newtonian 
fluid with imposed Poiseuille flow at inlet. 
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Figure 3 - Entry flow centerline temperatures at various Peclet 

numbers . 
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Studies of Metals Electroprocessing in Molten Salts 
Principal Investigator: Donald R. Sadoway 


RESEARCH ABSTRACT 


Molten salts are very important solvent systems in which the electro- 
deposition of a wide variety of metals may be conducted. However, solid 
electrodeposits from molten salts are typically incoherent, powdery and/or 
dendritic. Fluid flow patterns in the electrolyte are being observed by a 
laser Schlieren imaging technique in order to determine how mass transport 
affects the morphology of the metal deposit and thus improve our understand- 
ing of the electrodeposition process. Studies are beinq conducted on tne 
same metal, zinc, both in aqueous electrolytes in which coherent solid 
electrodeposits are produced, and in transparent molten salt electrolytes. 
Process variables such as current density and composition of the electro- 
lyte are adjusted to change the morphology of the electrodeposit and 
thus to permit the study of the nature of electrolyte flow in relation to 
the quality of the electrodeposit. The results of this work will be help- 
ful in electrochemical cell design and will serve as a data base for sub- 
sequent mathematical modeling of fluid flow in such systems. 

RESEARCH REPORT 


The present study is unique in that it seeks to explain the poor 
quality of solid electrodeposits in molten salts through a consideration 
of the effects of fluid flow of the electrolyte. The same metal is being 
electrodeposited in separate experiments from aqueous solutions and molten 
salt electrolytes. Transparent cells allow observation of electrolyte 
circulation patterns by a laser Schlieren optical technique. 
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Zinc has been chosen as the metal to study. It can be deposited 
both from aqueous chloride and molten chloride electrolytes; furthermore, 
both types of electrolyte systems are transparent to visible light, an 
important property which permits the use of noninvasive optical techniques 
of flow measurement. 

As this project is the first of its kind at the Materials Processing 
Center, much effort was previously spent in building up an electroproces- 
sing laboratory. Facilities such as those for flow visualization by laser 
Schlieren imaging were acquired. These include a pneumatically supported 
research grade optical bench measuring 4* x 6’. The present optical ar- 
rangement is capable of recording the laser images in any of these formats: 
35 mm photographic film, 8 mm movie film, or 1/2-inch videocassette. Of 
course, the necessary electrochemical instrumentation such as a potentio- 
stat, waveform programmer, coulometer, and X-Y and strip chart recorders 
have also been acquired. 

Phase I of the zinc plating study was devoted to the electrorefining 
of zinc in aqueous acid chloride solutions. Working at room temperature 
without the need for furnaces and controlled atmosphere cells allowed 
concentration on perfecting the optics of the system. Following a survey 
of the literature on acid zinc chloride solutions, a study of the chemistry 
cf these solutions was conducted with the purpose of identifying the domi- 
nant ionic species present as a function of ZnCl 2 concentration, alkali 
chloride concentration, and pH. The purpose of this study was to identify 
the conditions under which zinc is present as a single ionic species. The 
literature survey showed that other studies of zinc electrodeposition have 
been confounded by ill-defined chemistries. In this project, two specific 
studies were performed. First, a single thermodynamic model of the aqueous 

acid zinc chloride system was constructed to calculate the distributions of 
2 + + 2 - 

Zn , ZnCl , ZnCl^* and ZnCl^ as functions of concentration and pH. Ther- 
modynamic data were those reported in the literature. 

Figure 1 shows that as long as one is restricted * ^he binary electro- 
lyte system, ZnC^-HgO, there is no concentration for which a single z.nc 
species dominates. As a result, a third component was added to the electro- 
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lyte in order to displace the complexation equilibria through the presence 

of excess chloride ion in the amount of 3 M NaCl . Under these conditions 

the concentration of ZnCl, varied from 0.01 M to 1 M, with zinc present 

c 2 - 

only as fully complexed ZnCl^ . Secondly, titrations were performed to 
determine the operating range of pH such that the solutions were free of 
oxychlorides and hydrc.nde complexes. It was concluded that the optimum 
pH was in the vicinity of 2-3. 

Phase II of the work was the study of the electrochemistry of zinc 
reduction. Many contradictions in electrochemical behavior were discov- 
ered in the literature. These, in many cases, were attributable to in- 
complete control of bath chemistry. Measurements of critical overpoten- 
tial were made potentiostatical ly. Figure 2 shows a typical current vs. 
time curve from such an experiment. Figure 3 shows r .erived quantity, 
the induction time, plotted as a function of applied potential for two 
different, concentrations of ZnCl,,. The critical overpotential, n c * of 
65 mV is defined by the intersection of the two curves on the plot. Con- 
firmation of this quantity was given by similar plots of other time con- 
stants of the system. 

The third phase of the work was the study of convective r.uw during 
electrodeposition. Electrorefining experiments were conducted in cells 
of three different electrode configurations, shown in Figure 4. In the 
TF (top free) arrangement, the anode ’• located below the cathode and is 
facing downward so that its dissolution does not disturb the electrolyte 
flow pattern neat the cathode. The AT (anode on top) configuration com- 
bines the two convective flows, one due to cathodic deposition and the 
other due to anodic dissolution. Finally, in the G (for grid) configura- 
tion, the cathode is no longer a physical barrier to the falling electro- 
lyte. Deposition was performed potentiostatically. Current was recorded 
for a wide range of overpotentials and compositions. Flow patterns in the 
electrolyte were observed by laser Schlieren imaoing. During the experi- 
ments, photographs were made in synchronization with an event marker on 
the chart paper to permit correlations between the electrochemistry and 
fluid dynamics. 
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Figure 2. Current-time plot for the d .ermination of the critical 
overpetential , n c * in Solution I. Applied overvoltage 
was 150 mV. 
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Figure 5. Electrolyte flow patterns during potentiostatic deposition of zinc 
from solution II, n = 250 mV, AT ei r '•ode configuration. Anode 
surface is at the top of the picto e ame; cathode surface is at 
the bottom. 


Figure 5 shows the compilation of the current-time plot and the 
laser Schlieren photographs for an experiment in which the AT electrode 
configuration was used. The salient features are as follows. First, 
dissolution at the anode occurs before deposition at the cathode. The 
plane of intersection of ascending and descending plumes is closer to 
the cathode than to the anode. The anodic plumes, which by this time 
have already started necking, push their rising counterpart outward and 
soon a toroidal convection cell is established. The electrolyte flows 
downward In the center and upward at the edge. Quantitative results of 
such quantities as steady-state current and various critical time constants 
were also obtained below and at the limiting current. Correlations of 
these quantities with concentration agree well with ,'csults in the litera- 
ture. They indicated that the Rayleigh numbers involved were high enough 
(In excess of 5000) tint convection is turbulent and independent of elec- 
trode dimensions. 

These encouraging results are now being extended to molten salt cells. 
This will allow the comparison between electrolyte systems. Furthermore, 
through liaison with NASA staff at the Marshall Space Flight Center, thought 
is being given to studying low gravity electrodeposition. This would fur- 
ther reveal the role of turbulence and convective forces in electrodeposi- 
tion operations. 

LIST OF PUBLICATIONS 

1. A. A. Abdelmassih, "Effect of Fluid Flow on Zinc Electrodeposits 

from Acid Chloride Electrolytes", S. M. Thesis, M. I. T., February 1982. 

2. P. T. Rogers, "Fluid Flow Studies of Molten Salt Electrolytic Cells by 
Laser Schlieren Imaging", S. 8. Thesis, M. I. T. , June 1982. 

3. A. A. Abdelmassih and D. R. Sadoway, "Laser Schlieren Studies of 
Aqueous Zinc Chloride Electrolysis", in Chloride Electrometallurgy , 

P. D. Parker and P. Duby, editors, TMS-AIME, Warrendale PA, 1983. 


294 



Adaptive Materials Processing 


Adaptive Control of Welding Processes 
Principal Investigator: Thomas W. Eagar 


RES EARCH ABSTRACT 


The objective of this research is to use modern digital signal analy- 
sis techniques to determine natural oscillations occurring In the welding 
arc voltage and to attempt correlation of these signals with defects In the 
welds. Work was Initiated In the spring of 1981 by purchase of a dedicated 
minicomputer. During the previous reporting period, the computer was 
brought on-line and a number of test welds were made which proved that 
disturbances to the welding process could be detected In either the time 
domain or In the frequency domain. Unfortunately the noise levels from 
the commercial welding power supply were of the same order of magnitude as 
the measured disturbances. In the past year, a more sophisticated welding 
power supply has been designed and built, computer software and hardware 
have been upgraded, and signals from additional welds have been analyzed. 
The results continue to be encouraging, but the optimum signal analysis 
technique has not yet been determined. 


RESEARCH REPORT 


The two greatest problems Identified In the previous report period 
were the extraneous noise generated by the commercial welding power supply 
end the distortion caused by recording the welding voltage on magnetic tape. 
Both of these problems have been overcome In the current report period. 

ne first problem was overcome by acquiring a 600 amp DC (rms) 
transistorized current regulator manufactured by Alexander Kusko, Inc. 
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This was delivered In January 1932; and after a period of initial testing, 
the regulator has proven extremely useful. The regulator can maintain a 
welding current constant to within 1%, thus permitting a high signal-to-noise 
ratio. The regulator is linear to more than 20 kHz which permits the use 
of precisely pulsed current when necessary. Figure 1 compares the current 
stability of a commercial constant current welding power supply with the 
regulator obtained with the Kusko supply. These tests were made on a large 
carbon block resistor. The coarse steps of +5 amperes are due to the 
sensitivity of the A/D converter used to record the signal. The actual 
regulation of the Kusko supply is better than +0.5 ampere at this current 
level . 

The second problem of distortion of the signal by the AM magnetic 
tape has been corrected by software modifications to the computer which 
allow a direct 10 kHz sampling rate. Figure 2 shows a square wave as 
measured by the new direct input technique compared to the same square wave 
as recorded on the magnetic tape. The direct input technique is clearly 
superior. 

In addition to these software improvements for direct input of the 
welding process data, a complete signal analysis package was acquired from 
Signal Technology, Inc. during the summer. This package has saved consider- 
able time during the past few months by permitting emphasis on study of 
welds rather than on software development. 

The initial welds were produced with constant welding current. 

Although this technique did produce measurable signals during the previous 
report period, recent tests have shown that this technique is not as 
sensitive as an impulse or square wave current signal. Figure 3 shows 
the frequency response of the welding voltage on a stainless steel plate 
and an aluminum plate at constant current. Although there are minor dif- 
ferences In the signal, the spectrum is generally only white nosie. 

Figure 4 shows the spectrum on the same materials using a square wave 
current generated by the Kusko transistorized power supply. The harmonics 
of the square wave input are clearly observed, but it will be noted that 
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the envelope surrounding the peaks of these harmonics is different for 
steel and for aluminum. A simple square wave response generates harmonic 
peaks which decay exponentially. This exponential decay is distorted for 
both aluminum and steel but in a different pattern for each. These pattern" 
are very reproducible. 

It was also noted that the DC arc voltage level varies between steel 
and aluminum. This is shown In Figure 5 which represents the shift in arc 
voltage when welding over a steel -aluminum laminate. These tests show, 
at least in princl-^e, that welding over regions with chemical variations 
results in changes In the arc voltage. Studies have begun to test changes 
in voltage caused by welding over slag inclusions, but the results are 
inconclusive at present. 

Figure 6 shows the change in arc voltage caused by minor disturbances 
in the shielding gas. These disturbances are more readily observed In 
the time domain (Figure 6) than in the frequency domain (Figure 7). 


CONCLUSIONS 


Progress has been made In both the welding equipment and in the signal 
analysis techniques during the past year. The problems have been solved 
sufficiently that the project can concentrate on analysis o actual weldments 
during the current year. This will permit analysis of the frequency response 
of specific defects, as well as determi nation of the ultimate sensitivity 
of this control technique. In addition, welding signals will be analyzed 
by Nyquist plots in order to determine the transfer function for the process. 
These results will permit an accurate evaluation of the ultimate usefulness 
of this technique of weld quality monitoring. 
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